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Abstract 



By taking advantage of the laws of physics it is possible to revolutionize the way 
we communicate (transmit), process or even store infomiation. It is now known 
that quantum computers, or computers built from quantum mechanical elements, 
provide new resources to solve certain problems and perform certain tasks more 
efficiently than today's conventional computers. However, on the road to a com- 
plete understanding of the power of quantum computers there are intermediate 
steps that need to be addressed. The primary focus of this thesis is the under- 
standing of the possibilities and limitations of the quantum-physical world in the 
areas of quantum computation and quantum Information processing. 

First I investigate the simulation of quantum systems on a quantum computer 
(i.e., a quantum simulation) constructed of two-level quantum elements or qubits. 
For this purpose, I present algebraic mappings that allow one to efficiently ob- 
tain physical properties and compute correlation functions of fermionic, anyonic, 
and bosonic systems with such a computer. By studying the amount of resources 
required for a quantum simulation, I show that the complexity of preparing a quan- 
tum State which contains the desired Information is crucial at the time of evalu- 
ating the advantages of having a quantum computer over a conventional one. As 
a small-scale demonstration of the validity of these results, I show the simulation 
of a fermionic system using a liquid-state nuclear magnètic resonance (NMR) de- 
vice. 

Remarkably, the conclusions obtained in the area of quantum simulations can 
be extended to general quantum computations by means of the notion of gen- 
eralized entanglement. This is a generalization based on the idea that quantum 
entanglement (i.e., the existence of non-classical correlations) is a concept that 
depends on the accesible Information, that is, relative to the observer. Then I 
present a wide class of quantum computations that can be efficiently simulated 
on a conventional computer and where quantum computers cannot be claimed to 
be more powerful. The idea is that a quantum algorithm, performed by applying 



a restricted set of gates which do not create generalized entangled states relative 
to small (polynomially-large) sets of observables, can be imitated using a similar 
amount of resources with a conventional computer. However, a similar statement 
cannot be obtained when generalized entangled states (relative to these sets) are 
involved, because this purely quantum phenomena cannot be easily reproduced 
by classical-information methods. 

Finally, I show how these concepts developed from an information-theory 
point of view can be used to study other important problems in many-body physics. 
To begin with, I exploit the notion of Lie-algebraic purity to identify and charac- 
terize the quantum phase transitions present in the Lipkin-Meshkov-Glick model 
and the spin-1/2 anisotropic XY model in a transverse magnètic field. The results 
obtained show how generalized entanglement leads to useful tools for distinguish- 
ing between ordered and disordered phases in quantum systems. Moreover, I dis- 
cuss how the concept of general mean field hamiltonians naturally emerges from 
these considerations and show that these can be exactly diagonalized by using a 
conventional computer. 

In brief, in this thesis I apply several tòpics developed in the context of quan- 
tum Information theory to study the complexity of obtaining relevant physical 
properties of quantum systems with a quantum computer, and to study different 
physical processes in quantum many-body systems. 
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Chapter 1 
Introduction 



...there is plenty of room to make computers smaller...nothing that I 
can see in the physical laws... 

R. P . Feynman, Caltech (1959). 

During the last few dècades, the theory of Quantum Information Processing 
(QIP) has acquired great importance because it has been shown that information 
based on quantum mechanics provides new resources that go beyond the tradi- 
tional 'classical information'. It is now known that certain quantum mechanical 
Systems, named quantum computers (QCs), can be used to easily solve certain 
problems which are difficult to solve using today's conventional or classical com- 
puters CCs. Having a QC would allow one to communicate in secret [IBB84II 
(quantum cryptography), perform a variety of search algorithms IIGro97L factor 
large numbers IISho94l . or simulate efficiently some physical systems HOGKOIL 
ISOG02I . Additionally, it would allow us to break security codes used, for ex- 
ample, to secure internet Communications, optimize a large variety of scheduling 
problems, etc, which make of quantum information an exciting and relevant sub- 
ject. Consequently, the science of quantum information is mainly focused on bet- 
ter understanding the foundations of quantum mechanics (which are different of 
classical mechanics) and the physical realization of quantum controUable physical 
devices. While the first allows clever and not so obvious ways of taking advantage 
of the quantum world, the latter will let us achieve our most important goal: the 
building of a QC. 

When one looks for the word information in the dictionary one finds many def- 
initions: i) a message received and understood, ii) knowledge acquired through 



2 



Introduction 



study or experience, iii) propagated signal through a given channel, iv) broad- 
casted news, and more. Information is then the basis of all human knowledge and 
we usually base our behavior on it. It always requires a physical representation 
to be able to use it, propagate it, or store it, such as a telephone, a computer disk, 
etc. Depending on the physical representation, Information can be classical or 
quantum. 

We define as classical information the one that is manipulated and stored by 
today's CCs. In classical information theory the bàsic unit is the bit. A bit's state 
can be in one of two states represented by the numbers and 1 , which constitute 
the logical basis. A possible physical representation of a bit is given by a system in 
which the state is determined by the distribution of, for example, electrical charge. 
The idea is then to process information through the manipulation of the state of 
a set of bits (i.e., bit sequence) by performing elementary gates. These gates are 
different processes that depend on the particular physical realization of the CC. 

Examples of one-bit gates are the not and reset gates, and of two-bit gates 
is the nand gate. Their action over logical initial states are shown in Fig. II. 1[ 
They suffice for implementing arbitrary state transformations. That is, any classi- 
cal algorithm can be implemented through a circuit that consists of applying these 
elementary gates to a bit sequence. In fact, this method is used by today's comput- 
ers, where the pro gram sets up a particular order for performing the elementary 
gates and the chips implement them physically. Finally, one reads the final state 
where the required information is supposed to be encoded (e.g., the solution to a 
problem). 

The idea of quantum information processing is similar to that of classical in- 
formation but under the laws of the quantum world. One defines quantum infor- 
mation as the one which is stored and manipulated by physical devices obeying 
the laws of quantum physics; that is, satisfying the Schròdinger evolution equation 

thj^\^P) = H\i;), (1.1) 

where H is the Hamiltonian describing the interactions that one manipulates to 
perform the desired evolution, and \^) is some pure state (i.e., wave function) of 
the system. Such devices constitute quantum computers. 

In the conventional model of quantum information theory (QIT) the bàsic unit 
is the quantum bit or qubit. A qubit's pure state can be in any superposition of the 
logical states and is expressed as a|0) + ò| 1), where the complex numbers a and h 
are the probability amplitudes of being in the states |0) and respectively. They 
are normalized to the unity: |ap + |òp = l. A possible physical representation 
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Figure 1.1: Logical action of the single bit gates not and reset, and the two-bit 
gate nand. Here is and fs denote the initial and final bit states, respecti vely. 

of a qubit is given by any two-level quantum system (Fig. 11.21) such as a spin- 
1/2, where the state represented by |0) (|1)) corresponds to the state with the spin 
pointing up (down), or a single atom. Bohm's rule HBohSlll telis one that such 
state corresponds to having probabilities |ap and |òp of being in the state with the 
spin pointing up and down, respectively. 

Due to the superposition principle of quantum physics, a pure state of a set of 
qubits (register) is expressed as ao|0 ■ ■ ■ 00)+ai|0 ■ • • 01) + - ■ -+02^^-111 ■ ■ • 11). 
Again, the complex coefficients are the corresponding probability amplitudes, 
and J2i = 1- The idea is then to perform computation by executing a quan- 
tum algorithm that consists of performing a set of elementary gates in a given 
order (i.e., a quantum circuit). The action of these quantum gates is rigorously 
discussed in Chap. El and requires some previous knowledge in linear àlgebra. 
As in classical Information, these gates involve single qubit and two-qubit opera- 
tions. In order to preserve the features of quantum physics, these operations must 
be reversible (i.e., unitary operations), and are usually performed by making the 
register interact with extemal oscillating electromagnètic fields. 

The great advantages of having a QC are two-fold: First, working at the 
quantum level allows one to make these computers extremely small and, if scal- 
able' IIDiv95L it would allow one to process a large number of qubits at the same 



A Computer is said to be scalable if the number of resources needed scale almost linearly with 
the problem size. 
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X 



Figure 1.2: Different physical representations of the single qubit state \ip) = 
a\0) + b\l). The red dot on the surface of the sphere (left) represents such a 
linear combination of states. Here \g) and |e) denote the ground and an excited 
State of a certain atom. 

time. Second, a computer ruled by the laws of quantum physics should contain 
certain features that go beyond those of classical Information, since the latter can 
be considered a limit of the first one^. For example, one immediately observes that 
the superposition principle gives one more freedom when manipulating quantum 
Information, in the sense that many different logical states can be carried simulta- 
neously (parallelism). 

To analyze the computational complexity to solve a certain problem one needs 
to determine the total amount of physical resources required, such as bits or qubits, 
the number of operations performed or number of elementary gates, the number 
of times that the algorithm is executed, etc. While nobody knows yet the power 
of quantum computation, certain algorithms IISho94[ IGro97l suggest that QCs 
are more powerful than their classical anàlogues. All these algorithms share the 
feature that they not only make use of the superposition principle (which is not 
sufficient to claim that a QC is more efficient), but also of the non-classical corre- 
lations between different quantum elements in the QC. (Interference phenomena 
also plays an important role in the efficiency of quantum algorithms.) Such cor- 
relations are inherent to quantum systems IISch35[ IEPR35II and do not exist in 
classical systems. They are usually referred as quantum entanglement (QE), an 
emerging field of QIT. 



^Any classical algorithm can be simulated efficiently with a QC INCOOI 
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In order to access the quantum information, one needs to perform a measure- 
ment. This is defined as the extraction of some classical information from the 
quantum register. Due to the features of quantum physics, after a measurement 
process the state of the register is collapsed into the logical state corresponding to 
the outcome, with statistics given by Bohm's rule (Fig. 11.31) . This process could 
destroy the efficiency of the computation. For example, if after the execution of 
the quantum algorithm the state of two qubits is ao|00) + ai|01) + 02] 10) +03! 11), 
a measurement process in the logical basis has the effect of coUapsing the state to 
|00) with probability |aop, to |01) with probability |aip, to |10) with probability 
|a2p, and to |11) with probability \a-i\^. Therefore, a single measurement does not 
give the whole information about the state of the register and usually one needs to 
run the quantum algorithm repeatedly many times to obtain more accurate statis- 
tics to recover the state of the register. This is a main difference with the case of 
classical information, where such measurement or convertion is not necessary. 
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Preparation 
of Initial 
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Evolutíon 








4^ 








Unitary 








operators 






Figure 1.3: Circuit representation of a quantum algorithm. A pure state \ip) is 
evolved by applying elementary gates. The evolution obeys Schròdinger's equa- 
tion (Eq. ll.ll) . After the evolution, a measurement coUapses the evolved state with 
statistics given by Bohm's rule. 

The simplest case of an entangled state is the pure two-qubit state (Fig. II. 41) 

|V;) = Í=[|10) + |01)], (1.2) 
or similar states obtained by flipping or by changing the phase of a single qubit. 
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Equation il .21 states that if one qubit is measured and projected in the logical ba- 
sis, the other qubit is automatically projected in the same basis. Moreover, if 
the outcome of a measurement performed in one qubit is (1), the outcome of a 
post-measurement performed on the other qubit will be 1 (0). Remarkably, sim- 
ilar results are obtained for such state if the measurements are performed in a 
basis other than the logical one. These correlations between outcomes cannot be 
explained by a classical theory. 




Figure 1.4: Maximally entangled two-qubit state. The quantum correlations can- 
not be represented by any classical state. 

During the last few years, several authors IIEHK()4[IGKÍIS()K()3llW"()21IVid()3l 

have started to study the relation between different definitions and measures of QE 
and quantum complexity. Naturally, they mostly agreed that whenever the QE of 
the evolved state in a quantum computation is small enough, such algorithms can 
be simulated with the same efficiency on a CC. However, the lack of a unique 
computable measure of entanglement that could be applied to any quantum state 
and quantum system is the main reason why the power of QCs is still not fuUy 
understood. 

One of the purposes of my thesis is to show the computational complexity 
(i.e.., the number of resources and operations needed) to solve certain problems 
with QCs and to compare it with the corresponding classical complexity. In par- 
ticular, I will first focus on the study of the simulation of physical systems by 
quantum networks or quantum simulations (QSs) IIOGK()lllSOÜ02llSÜK03l . As 
noticed by R. R Feynman [|Fey82|| and Y. Manin, the obvious difficulty with de- 



terministically solving a quantum many-body problem (e.g., computing some cor- 



relation functions) on a CC is the exponentially large basis set needed for its sim- 
ulation. Known exact diagonalization approaches like the Lanczos method suffer 
from this exponential catastrophe. For this reason, it is expected that using a 
Computer constructed of distinct quantum mechanical elements (i.e., a QC) that 
'imitates^' the physical system to be simulated (i.e., simulates the interactions) 
would overcome this difficulty. 

The results obtained by studying the complexity of QSs can be extended to un- 
derstand the complexity of solving other problems. For example, different quan- 
tum search algorithms IIGro97l admit a Hamiltonian representation and can be 
equivalently considered as a particular QS. But most importantly, I will show how 
these simulations lead to the definition of a general measure of quantum (pure 
State) entanglement, so called generalized entanglement (GE), which can be ap- 
plied to any quantum state regardless of its nature. Remarkably, this measure is 
crucial when analyzing the efficiency-related advantages of having a QC. 

This report arises from the studies and novel results obtained together with my 
coUeagues at Los Alamos National Laboratory (USA) and Instituto Balseiro (Ar- 
gentina) during my PhD studies. For a better understanding, the main results are 
presented in chronological order. In Chap.|21 1 analyze the problem of simulating 
different finite physical systems on a QC using deterministic quantum algorithms. 
The corresponding computational complexity is also studied. As a proof of prin- 
cipies, I present the experimental simulation of a particular fermionic many-body 
system on a liquid-state nuclear-magnetic-resonance (NMR) QC. 

In Chap. |3l I introduce the concept of quantum generalized entanglement, a 
notion that goes beyond the traditional quantum entanglement concept and makes 
no reference to a particular subsystem decomposition. I show that important re- 
sults are obtained whenever a Lie algebraic setting exists behind the problem 
under consideration. In particular, I apply this novel approach to the study of 
quantum correlations in different quantum systems, regardless of their nature or 
particle statistics, including different spin and fermionic systems. 

In Chap. mi compare the effort of simulating certain quantum systems with a 
QC or a CC. In particular, I show that the concept of generalized entanglement is 
crucial to the efficiency of a quantum algorithm and can be used as a resource in 
quantum computation. Moreover, generalized entanglement allows one to make a 
connection between QIT and many-body physics by studying different problems 

^In general, the QC used to perform a QS is built of quantum elements that are different in 
nature of those that compose the system to be simulated. However, this is not a drawback in a 
simulation because usually one can perform a one-to-one association between the quantum states 
of the QC and the quantum states of the physical system to be simulated. 
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in quantum mechanics, such as the characterization of quantum phase transitions 
in matter or the study of integrable quantum systems. These results are presented 
in Chap.|5l 

Finally, in Chap. |6l I present the conclusions, open qüestions, and future di- 
rections related to this subject. 



Chapter 2 



Símulatíons of Physics with 
Quantum Computers 

Since Richard P. Feynman conjectured that an arbitrary discrete quantum system 
may be simulatedhy another one | |Fey82| , the simulation of quantum phenomena 
became a fundamental problem that a quantum computer (QC), i.e., a system of 
universally controUed distinct quantum elements, may potentially solve in a more 
efficient way than a classical computer (CC). The main problem with the simula- 
tion of a quantum system on a CC is that the dimension of the associated Hilbert 
space grows exponentially with the volume of the system to be simulated. For 
example, the classical simulation of a system composed of N qubits requires, in 
general, an amount of computational operations (additions and products of com- 
plex numbers) that is proportional io D = 2^, where D is the dimension of the 
Hilbert space given by the number of different logical states \i1i2 ■ ■ - Ín), with 
ij = {0, 1}. Nevertheless, a QC allows one to imitate the evolution of the corre- 
sponding quantum system by cleverly controUing and manipulating its elements. 
This process is called a quantum simulation (QS). It is expected then that the num- 
ber of resources required for the QS increases linearly (or at most, polynomially) 
with the volume of the system to be simulated IIAL97II . If this is the case, we say 
that the QS can be performed efficiently. 

To be able to perform a QS, it is necessary to make a connection between the 
operator àlgebra associated to the system and the operator àlgebra which defines 
the model of quantum computation [OGKOl |. The existence of one-to-one map- 
pings between different algebras of operators and one-to-one mappings between 
different Hilbert spaces fB QÜTllSOKOSI . is a necessary requirement to simulate 
a physical system using a QC built on the basis of another system (Fig. l2.1l ). For 
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example, one can simulate a fermionic system on a liquid-state nuclear magnètic 
resonance (NMR) QC by making use of the Jordan-Wigner transformation IIJW28I 
that maps femiionic operators onto the Pauli (spin-1/2) operators. Although these 
mappings can usually be performed efficiently, this is not sufficient to establish 
that any system can be simulated efficiently on a QC. It is then necessary to prové 
that all steps involved in the QS, including the initialization, evolution, and mea- 
surement, can be performed efficiently IIS0GÜ2II . 




Figure 2.1: Relationship between different models of computation (with their as- 
sociated operator algebras) and different physical systems. Question marks refer 
to the present lack of a quantum computer device using the corresponding elemen- 
tary physical components indicated in the box. Diamond-shaped arrows represent 
the natural connection between physical system and operator language, while ar- 
rows on the circle indicate the existence of isomorphisms of *-algebras, therefore, 
the corresponding simulation of one physical system by another. 

This chapter will explore the theoretical and experimental issues associated 
with the simulations of physical phenomena on QCs. In Sec. 12.11 I start by de- 
scribing different models of quantum computation. In particular, I rigorously in- 
troduce the conventional model by means of the Pauli operators, where a natural 
set of elementary gates (i.e., set of universal operations) is obtained. This model, 
roughly described in Chap.[Tl is the one generally needed for the practical imple- 
mentation of a QS. In Sec. 12.21 1 present a class of quantum algorithms (QAs) in 
the language of the conventional model, for the computation of relevant physical 
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properties of quantum systems, such as correlation functions, energy spectra, etc. 
In Sec. 12.31 1 explain how the QS of quantum physical systems obeying fermionic, 
anyonic, and bosonic particle statistics, can be performed on a QC described by 
the conventional model, presenting some mappings between the different opera- 
tor algebras. As an application, in Sec. 12.41 1 show the QS (imitated by a classical 
Computer) of a particular fermionic system: The two-dimensional fermionic Hub- 
bard model. It is expected that such simulation gives an insight into the limitations 
of quantum computation, showing that certain issues remain to be solved to as- 
sure that a QC is more powerful than a CC (Sec. 12.51) . In Sec. 12.71 I describe 
the experimental implementation on an NMR QC of the QS of another fermionic 
system: The Fano-Anderson model. For this purpose, an elementary introduction 
to the physical processes on an NMR setting is described in Sec. 12.61 Finally, I 
summarize in Sec. 12.81 

2.1 Models of Quantum Computation 

When performing a quantum computation, the quantum elements which constitute 
the QC can be universally controUed and manipulated by modulating and chang- 
ing their interactions. This quantum control model assumes then the existence of 
a control Hamiltonian Hp, which describes these interactions. The control possi- 
bilities are used to implement specific quantum gates, allowing one, for example, 
to represent the time evolution of the physical system to be simulated HOGKOll . 

In order to define a model of quantum computation it is necessary to give a 
physical setting together with its initial state, an àlgebra of operators associated to 
the system, a set of controUable Hamiltonians necessary to define a set of elemen- 
tary gates, and a set of measurable operators (i.e., observables). In this way, many 
different models of quantum computation can be described, but for històric rea- 
sons and practical purposes I will focus mostly on the conventional model UNCOOI . 

2.1.1 The Conventional Model of Quantum Computation 

As mentioned in Chap.Ul in the conventional model of quantum computation, the 
fundamental unit of Information is the quantum bit or qubit. A qubit's pure state 
|a) = a|0) + ò|l) (with a, ò e C and |ap + |6p = 1), is a linear superposition 
of the logical states |0) and and can be represented by the state of a two- 
level quantum system such as a spin-1/2. Assigned to each qubit are the identity 
operator 11 (i.e., the no-action operator) and the Pauli operators a^, ay, and a^. In 
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the logical single-qubit basis Í3 = {|0), |1)}, these are 



10\ /01\ (0 MO 

1 ' = 1 ' = U ' = -1 



Because of its action over the logical states, the operator is usually referred as 
theflip operator: 

(22) 

For practical purposes in this thesis, it is also useful to define the raising (+) and 
lowering (-) Pauli operators a± = |((T^. ± iay), and the eigenstates of the flip 
operator |+) = ^[|0) + |1)] and |-) = ^[|0) - |1)], satisfying 

a,|±) = ±|±). (2.3) 

The Pauli operators form the su(2) Lie àlgebra and satisfy the commutation 
relations (/x, z/, A = {x, y, z}) 



i2ef,uxcr\, (2.4) 



where [A, B] = AB — BA and e^^x is the total anti-symmetric Levi-Civita symbol. 
They constitute a complete set of local observables, that is, a basis for the 2x2 
dimensional Hermitian matrices with = (a^)^. The symbol f denotes the 
corresponding complex conjugate transpose. 

Any qubit's pure state can be represented as a point on the surface of the 
unit sphere (Bloch-sphere representation) by parametrizing the state as |a) = 
a|0)+ò|l) = cos(^/2)|0)+e^'^sin(^/2)|l) (Fig.ESl). In order to process a single 
qubit, a complete set of single-qubit gates has to be given. These operations con- 
stitute then, any rotation in the Bloch-sphere representation, which are given by 
the operators R^{^) = e"*^'^/^^''" = cos('i9/2)]l — z sin(í?/2)(T^; that is, a rotation 
by an angle 'í9 along the yu axis. These rotations are unitary (reversible) operations, 
satisfying i?^(í9)[i?^('í9)]t = 1 (i.e., no-action), where = i?^(-ï9). This 

reversibility property allows one to perform these gates with no thermodynami- 
cal cost. In Fig. 12.31 1 present these elementary single-qubit gates in their circuit 
representation. 

Similarly, a pure state of a A^-qubit register (quantum register) is represented 
as the ket = Yln=o^ CLn\n), where |n) is a product of states of each qubit in the 
logical (or other) basis, e.g., its binary representation (|0) = IO1O2 ■ ■ ■ Oat), |1) = 
IO1O2 ■ ■ ■ On-Ún), |2) = IO1O2 ■ ■ ■ 1n-iOn), etc), and J^Ho' = 1 («n G Q. 
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Figure2.2: Bloch-Sphere representation of a one qubit state parametrized as |a) = 
cos(6'/2)|0) + e**^ sin(6'/2)|l). The curved arrows denote rotations i?^^ along the 
corresponding axis. The (arrow) color convention is: |0) ^blue; |1) ^ red; other 
linear combinations magenta. 

Assigned to the jth qubit of the quantum register are, together with the identity 
operator F, the local Pauli operators (with n = x,y, or z); that is 

n factors 

, ^ ^ 

factor 

where ® represents a Kronecker tensorial product. Their matrix representation in 
the basis ordered as i3 = {|0i ■ ■ ■ OAr_iOAr), |0i ■ ■ ■ Oat.iIat), ■ ■ ■ , |1i ■ ■ ■ lAr_ilAr)} 
is just the matrix tensor product of the corresponding 2x2 matrices defined by 
Eq. 12.11 For two different qubits, these operators commute: 

fe^'] =OVjVA;. (2.5) 

In order to describe a genèric operation on the quantum register, it is also nec- 
essary to consider products of the Pauli operators cr^. Remarkably, every uni- 
tary (reversible) operation acting on the quantum register can be decomposed 
in terms of single qubit rotations K'^{'d) = e~*2<^í' and two-qubit gates, such 

as the Ising gate R^j^k^uj) = e~*^'^^°"í = cos(cj/2)]1 — isin{uj/2)aia^, with 
G M (rBBC931 DÍV95I ). The operations R^j^k{Lü) are also unitary, satisfy- 

ing R^j^k{u)[Rziz*'i^W — with [_R^j^fe(ci;)]'l" = R-^j .^k{—üj). Together with the 
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Figure 2.3: Circuit representation of the elementary gates in the conventional 
model. The top picture indicates a single-qubit rotation while the bottom one 
indicates the two-qubit Ising gate. Any quantum algorithm can be represented by 
a circuit composed of these elementary gates. 

single-qubit rotations they define a universal set of quantum gates. Their quantum 
circuit representation is shown in Fig |2.3[ 

Every (logical) state of the quantum register has associated a mathematical 
object denoted as bra in the foUowing way 

\n) ^ {n\, (2.6) 

and can be linearly extended for a general state \ip) by conjugation as 

^a„|n) ^ ^a;(n|, (2.7) 

n n 

where a* denotes the complex conjugate of a„. The product between a logical bra 
and a logical ket defines the inner-product in the associated Hilbert space given 
by 

(m\\n) = {m\n) = ó^n, (2.8) 

with 5mn being the Kronecker delta. In this vectorial space, two vectors (states) 
\tp) and 10) are orthogonal if their overlap, that is, their inner product (i^lcp) given 
by Eq. 12.81 vanishes. Moreover, the bra-ket notation allows one to represent every 
(Pauli) operator. For example, the single qubit flip operator is represented as 
Cx = |0)(1| + |1)(0|. This notation is very useful when computing, for example, 
expectation vàlues. 
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A measurement is defined as the action that gives some classical information 
about the state of the quantum register. In quantum mechanics, a measurement 
is considered to be a probabilistic process that collapses the actual quantum state 
of the system IIPer98ll . For example, a measurement of the polarization in the 
logical basis of every qubit (i.e., the measurement of the observables aí) when 
the state of the register is [tp) = J2n '^nln), projects it onto a certain logical state 
\m) with probability la^P (Bohm's rule). This is a von Neumann measurement. 
In particular, in a general von Neumann measurement of an observable M (= 
Aí^), the probability pm that the outcome m is obtained, where m is one possible 
eigenvalue of M, is given by INCOOI 

Pm = {^\MlM^\^), (2.9) 

where Mm is the projector onto the subspace of states with quantum number m. 
Moreover, if m is the actual outcome, the state after the measurement is given by 

1^') = , ^"^^ 1^)- (2.10) 

For example, when measuring the operator for a single qubit state |0), the 
two possible outcomes are m = ±1 (i.e., the eigenvalues of a^). Since |0) = 
75(1+) + the corresponding probabilities are 

Pi=p_i = l/2, (2.11) 

and if the outcome is +1 (—1), the state is projected onto |+) (|— )). Therefore, 
to obtain accurate information about the actual state of the quantum system, one 
needs to prepare many copies of the state and perform many different measure- 
ments. 

The expectation value of a measurement outcome is the expectation of the 
outcomes of many measurement repetitions. It can also be expressed in the bra- 
ket notation. If the state of the quantum system is |^), the expectation value of M 
is given by 

(M) = {^\M\^). (2.12) 

In the conventional model, any observable M can be written as a combination 
(sums and/or products) of the identity and Pauli operators. Therefore, if líp) is 
known, the expectation (M) can be algebraically computed by obtaining first the 
state M\íp), and by projecting it onto the bra {íjj \ using the inner-product relations 
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of Eq. lTHl For example, if a two-qubit State is givenby = :^[|0i02) + II1I2)], 
then 

(^k.V) = \ ((O1O2I + (I1I2I) (IO1O2) - II1I2)) = 0, (2.13) 

and 

{M^WM = \ ((O1O2I + (I1I2I) (II1I2) + IO1O2)) = 1. (2.14) 

Equations|TT3and|TT3|havebeen obtainedby noticing that a^|lil2) = — II1I2), 
a^a2|0i02) = |lil2),and(T^(T2|lil2) = IO1O2), together with Eq.El 

Nevertheless, certain quantum computations and QSs are done by evolving 
mixed states instead of pure states. A quantum register in a probabilistic mixture 
of pure states can be described in the bra-ket notation by a density matrix p = 
^sPsPs^ with Ps = \ips){'ips\ representing the quantum register being in the pure 
State \ijjs), with probability Ps (}2gPs = í;Ps > 0). Equivalently, every density 
operator p can also be written as a combination (sums and/or products) of the 
Pauli operators cx^ (a = x, y, z) and the identity operator t. These mixed states 
are useful when performing quantum computation with devices such as the NMR 
QC, where the state of the quantum register is approximated by the average state 
of an ensemble of molècules at room temperature; that is, an extremely mixed 
state. The expectation value of a measurement outcome over a mixed state is 
given by 

(M) = Tr(pM), (2.15) 

where M is the measured observable, p the density operator of the mixed state, 
and Tr denotes the trace. 

In brief, the conventional model allows one to describe every step in a QS by 
means of Pauli operators. The idea is to represent any quantum algorithm (QA) 
as a circuit composed of elementary single and two-qubit gates, together with the 
measurement process. The complexity of a QA is then determined by the amount 
of resources required, given by the number of qubits needed, the number of uni- 
versal single and two-qubit operations (Fig. 12.31) . and the number of measurements 
needed to obtain an accurate result (e.g., the number of times that the algorithm 
needs to be performed). For this purpose, a procedure to decompose an arbitrary 
operation in terms of elementary gates has to be explained. In the foUowing sub- 
section, I present some useful techniques and examples. 
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2.1.2 Hamíltonían Evolutíons 

When simulating a physical system on a QC it is necessary, in general, to perform 
a Hamiltonian (unitary) evolution to the quantum register HOGKÜ 1 [ IS 0GÜ2I . of 
the form 

U{t) = e-'^' = l-iHt+ ^{-iHtf + ■■■, (2.16) 

where H = H^^ is a physical Hamiltonian and t is a real parameter (e.g., time). A 
common H is given by 

H = H, + Hy = aal a'}j ai + P al al^ a^ , (2.17) 

where a and (3 are real numbers. From Eqs. 12.41 and l2.5l one obtains [Hx, Hy] = 0, 
and therefore, U{t) = ^-■'■Hxtç^-iHyt _ 

To decompose U (t) into single and two-qubit operations, the foUowing steps 
can be taken. First, the unitary operator 

U, = e'^''y = ^[l + ial] (2.18) 

takes al al, i.e., UlalUí = al, so f/|e*"'^^f/i = e*"""^. Second, the operator 

U, = e^^^'^-'^ = ^[l + ^alal] 

takes al a\al, so Vlé^^'^-V^ = e'^'^l"^-. Then, 

Us = e^ï'^^-^ 

takes (jyCr^ —ala'^a'^. By successively similar steps the required string of 
operators can be easily built: ala'^ ■ ■ ■ ai~^ai and also exp[ia(T^(j^ ■ ■ ■ ai'^a^] 
(up to a global irrelevant phase): 

f/^ ■ ■ f/2^t/íe'""'f/if/2 ■ ■ ■ f/fc = exp [iaalal ■ ■ ■ at^al] (2.19) 

where the integer k scales linearly with j. The evolution e~^^y^ can be decom- 
posed similarly so U (í) is decomposed as the product of both decompositions. 
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2.1.3 Controlled Operations 

Alternatively, one can use the well known ControUed-Not, or CNOT, gate instead 
of the two-qubit Ising gate. Its action on a pair of qubits (1 and 2) is 



CNOT 



Here, qubit 1 is the control qubit (the controlled operation on its state | li) is repre- 
sented by a sòlid circle in Fig. 12.41 ). If the state of qubit 1 is |0i) nothing happens 
(identity operation) but if its state is |li), the state of qubit 2 is flipped. The de- 
composition of the CNOT unitary operation into single and two-qubit operations 
is 

CNOT: é-ie-'-i^^^e-'-i^-é-i"^^"^-, (2.20) 



which was obtained by noticing that ie *2<^== = a'l, i.e., the spin-flip operator 
acting on qubit 2 (Eq. l2.1l) : 

a^lKti^hnl' (2.21) 

i mU) IS1O2) . 

By using the techniques described in Sec. 12.1.21 the CNOT operation in terms of 
single and two-qubit Ising gates is 

CNOT: é^e-'^''^^e-'^"''-é^"ye-'^"^-''^^e~'^''i . (2.22) 

The circuit representation of this decomposition is shown in Fig. 12.41 Five ele- 
mentary single and two-qubit Ising gates are required to perform the CNOT gate. 

These results can be extended to other controlled unitary operations like the 
CU operation defined as 

CU: 10)3(01 + ® f/s. (2.23) 

The unitary operation above performs the transformation Ug (also unitary) over 
a set of qubits s if the state of the control qubit a is and does not act other- 
wise. For the transformation Us = U{t) = e"*'^*, with Q = Q\ the operational 
representation of the CU gate is 



f/(t/2)f/(í/2)-"^ 



(2.24) 
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Figure 2.4: CNOT gate decomposition and its matrix representation. The control 
qubit is 1. Note that the last circuit realizes the CNOT matrix operation up to a 
global phase e~*ï. 



(a) 



(b) 




s — U{t) =4 U{t/2)U{t/2)-< 



s - U{t) =4 U{t/2)U{t/2)< 



Figure 2.5: (a) CU operation with the state of the control qubit a being in | l)a and 
(b) CU' operation controUed with the state |0)a. 
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where f/(t/2) ""^ = e"^®'^^ [Fig. I2.5r a)1. Equivalently, one can define another 
controUed operation, CU', on the state |0)a [Fig-Ob)]: U{t/2)U{t/2yi . 

ControUed operations are widely used in quantum algorithms. In general, their 
decomposition into single and two-qubit gates require a large number of these 
elementary operations, so they should be avoided when possible. 



2.2 Deterministic Quantum Algorithms 

In a QS, a QC performs certain tasks which are expected to give some Information 
about the physical system being simulated. These tasks are communicated by 
means of a program or quantum algorithm (QA), which can be schematically 
represented as a quantum circuit. In this section, I present a particular type of 
QA that can be used to obtain relevant properties of a quantum physical system s, 
using the conventional model (Sec. 12.11 ). Nevertheless, the same techniques can 
be used to simulate physical systems with other particle statistics (e.g., fermionic 
or bosonic systems), if they can be described by Pauli operators after an algebraic 
mapping. 

A deterministic QA is based on three different steps: i) the preparation of a 
pure initial state, ii) its evolution, and iii) the measurement of certain property of 
the evolved state, in which the result of the algorithm is encoded. To preserve the 
features of the quantum world, the evolution step is performed through a unitary 
operation and the measurement step is described by a certain observable (i.e., 
Hermitian operator). Here, I present only the class of QAs that allows one to 
determine, in a register of N qubits, physical correlation functions of the form 

G=(0|W^,|0), (2.25) 

where Ws is a unitary (reversible) operator associated to the system to be simu- 
lated; that is, H^sW^s^ = 11 (s refers to the system). 

Indirect measurement techniques can be used to obtain such correlation func- 
tions on a QC. In addition to the qubits used to represent the physical system s 
to be simulated (i.e., the qubits-system) extra qubits, called ancillas, are required. 
These ancillas constitute the probes that contain Information about the qubits- 
system. In the foUowing section I describe different measurement techniques. 
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2.2.1 One-Ancilla Qubit Measurement Processes 

In this case a single ancilla qubit allows one to obtain the correlation functions of 
Eq. 12.251 with Wg = WV, and U, V unitary operators acting on s HOGKÜll . For 
this purpose, the ancilla qubit a is first initialized in the State |+)a = -^(|0)a+|l)a) 
by applying, for example, the unitary Hadamard gate to the state |0)a^ Second, 
one makes it interact with the qubits-system, initially in certain pure state 
through two controUed unitary operations V and U, associated to the V and U 
operations, respectively. The first operation V evolves the system by V if the 
ancilla is in the state |1): V = |0)a(0| (8> 11^ + |l)a(l| ® The second operation U 
evolves the system by í/ if the ancilla state is |0): Ü = |0)a(0| |l)a(l| ® Is- 

Notice that V and U are reversible and commute with each other. 

After such evolution, the final state of the quantum register, |^/), is 

IV'/) = VÜ\+),\(P) = [|0)a ® í/|0) + |l)a ® Vm . (2.26) 

Interestingly, the expectation value of the Pauli operator 2a^ = cr^+iay associated 
with the ancilla qubit, in this state, gives the desired correlation function: 

G={^j\al + tal\^Pj), (2.27) 

where I have used the orthogonalitv propertv (Sec. 12. 1 .11) . that is. (0|l)a = (l|0)a = 
0, (0|0)a = (1 1 l)a = 1, and the action of the operators over the state of the an- 
cilla qubit given by Eq. l2.1[ The corresponding circuit for this quantum algorithm 
is shown in Fig. 12.61 Due to the probabilistic nature of quantum measurements, 
the desired expectation value is obtained with variance C(l) for each instance. 
That is, in order to get an accurate value for {2a^), repetition must be used to 
reduce the variance below what is required (Sec. l2.1.TT i. 

Nevertheless, sometimes it is necessary to compute the expectation value of 
an operator W of the form 

M 

W = Y,(^i UjVi , (2.28) 

i=l 

where Ui and Vi are unitary operators, > G M (with no loss of generality). 
In principle, this expectation value can be computed by preparing M different 
circuits as the one represented in Fig. 12.61 such that each algorithm computes 

'The Hadamard gate in terms of single qubit rotations is íe^*T·^!/e~*ï'^- 
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Figure 2.6: Quantum algorithm for the evaluation of the correlation G = 

m^v\<i>). 

(f//VÍ). However, in most practical cases, the preparation of the initial state 
is a very difficult task. This difficulty can then be reduced by using a particular 
QA that requires only one circuit, but with L ancilla qubits, where L = J + 1 and 
J > log2 M. Such QA has been described in Ref. fiSOG02il . The idea is to extend 
the results described above using controlled operations with respect to different 
ancilla qubits. 

2.2.2 Quantum Algorithms and Quantum Simulations 

Based on the indirect-measurement methods described in Sec. 12.2. li I now present 
certain QAs for QSs. These are useful for obtaining relevant properties of quan- 
tum systems, like the evaluation of the correlation function 

G{t) = (0|Tt4Ti?,.|0). (2.29) 

Here, Ai and Bj are unitary operators (any operator can be decomposed in a uni- 
tary operator basis as A = ohAí, B = J2 PjBj), T = e~*^* is the time evo- 

i j 

lution operator of a time-independent Hamiltonian H associated to the physical 
system to be simulated, and |0) is a particular state of the physical system. Notice 
that Eq. l2.29l is a particular case of Eq. l2!25l In particular, the evaluation of spatial 
correlation functions can be obtained by replacing the evolution operator T by the 
space translation operator. 

The quantum circuit for the evaluation of G(t) is shown in Fig. 12.71 It is 
equivalent to the one shown in Fig. I2.6l bv choosing W = T'^Aí and V = TBj. 
This particular selection allows one to reduce the complexity of the problem in 
the sense that the operation T need not to be controlled by the state of the ancilla 
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Figure 2.7: Quantum algorithm for the computation of spatial and time-correlation 
functions. In this case, {2(j%) = {(j)\T^AlTBj\(f)). Notice the simplification 
achieved by reducing the controUed-T operations. The same convention of 
Fig. l2.6l has been used. 



qubit ( IIS0GÜ2I ). As mentioned in Sec. l2.1.3l controUed operations require a large 
number of elementary gates, so they must be avoided when possible. 

In brief, the computation of G{t) is performed as foUows: First, the ancilla 
qubit a is prepared in the state |+)a, and the system is prepared in the state 
Second, a controlled evolution on the state |l)a, given by C-B = |0)a(0| ® + 
|l)a(l| (8) Bj, is performed. Third, the time evolution T is performed. Fourth, a 
controlled evolution on the state |0)a, given by C-A = |0)a(0| (g) + |l)a(l| ® Is, 
is performed. Finally the observable (2cr^) = {a^ + iay) = G{t) is measured. 

Sometimes one is interested in obtaining the spectrum (eigenvalues) of a given 
observable Q ( Hermitian operator), associated with some physical property of the 
system to be simulated. Techniques for getting spectral Information can be used 
on the quantum Fourier transform IIKit95llCEM98l and can be applied to physical 
problems IIAL99II . Nevertheless, the idea here is to use the methods developed in 
Sec. EO 

For some Hermitian operator Q, such as the Hamiltonian H of the system to 
be simulated, a common type of problem is the computation of its eigenvalues 
or, at least, the lowest eigenvalue (related, for example, to the ground state of the 
system). For this purpose, the qubits-system needs to be initialized in a certain 
state 10), that has a non-zero overlap with the eigenstates corresponding to the 
eigenvalues that need to be computed. Such a state can always be decomposed as 
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Figure 2.8: Quantum algorithm for the computation of the spectrum of an observ- 
able Q. Inthiscase, (2(7;) = (0|e-*'3*|</)). 

a linear combination of eigenstates of Q, 

\<P) = Y.^nm, (2.30) 

n 

where 7„ are complex coefficients, \tpn) are the eigenstates of Q, and A„ the cor- 
responding eigenvalues. If interested in computing Am, it is then required that 
\^m\ ^OinEq.ESOl 

For this purpose, the correlation function 

^(í) = (0|f/(í)|0), (2.31) 

with U (t) = e~*'^*, needs to be computed for different vàlues of the real parame- 
ter t (usually related with time). For a particular t, the measurement of S{t) can 
be performed using the one-ancilla method (Sec. l2.2.TT) as foUows. First, the ini- 
tial State |+)a ® is prepared. Second, the unitary evolution exTp[iQalt/2] is 
performed. Finally, the expectation {2a^) = S{t) is measured. The circuit rep- 
resentation of this QA is shown in Fig. 12.81 It is equivalent to the one shown in 
Fig.ESlby replacing W = V = e'^'^í/a 

For a particular value of t, the function S{t) is 

S{t) = J2hn\'e-'^"'. (2.32) 

n 

Then, the eigenvalues A„ can be obtained by performing a classical Fourier trans- 
form to Eq.|I221(i·e., ^(A) = / S{t)e'^'dt) 

S{X) = J22n\^n\'Si\-\n). (2.33) 
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However, S{t) is only obtained for a discrete set of vàlues of t. It is needed, in- 
stead, to calculate the corresponding discrete Fourier transform (see AppendixIXt 
to obtain information about the A„'s. 



2.3 Quantum Simulatíons of Quantum Physics 

In the most general case, a QS requires the simulation of systems with diverse 
degrees of freedom, like fermions, anyons, bosons, etc. The associated Hilbert 
spaces (space of states) differ from the one defined for the conventional model. 
For example, in the case of fermionic systems, fermionic states are governed by 
Pauli's exclusion principle. Then, at most a single spinless (or two spin-1/2) 
fermion can occupy a certain (atòmic) quantum state at the same time. There- 
fore, all the features associated with the physical system to be simulated must be 
preserved when transforming its operators to the operators describing the compu- 
tational model of the QC. 

In this section, I present isomorphic mappings that allow one to simulate ar- 
bitrary quantum systems, regardless of their particle statistics, by using the QAs 
defined for the conventional model (Sec. l2.2l l. Fortunately, such mappings can be 
easily performed without breaking the efficiency of a QA. 

2.3.1 Simulatíons of Fermionic Systems 

The systems considered here consist mainly of a lattice with modes (sites), 
where spinless fermions can hop between sites. These results can be easily ex- 
tended for the case of spin-1/2 fermions or higher spin fermions. 

In the second quantization representation, the (spinless) fermionic operators 
and Cj are defined as the creation and annihilation operators of a fermion in the 
j-th mode (j = 1, ■ ■ ■ , A^), respectively. Due to the Pauli's exclusion principle and 
the antisymmetric nature of the fermionic wave function under the permutation of 
two fermions, the fermionic àlgebra is given by the foUowing anticommutation 
relations 

{d, Cj} = 0, {4, Cj} = 5ij (2.34) 

where {, } denotes the anticommutator (i.e., {A, B} = AB + BA). 

The Jordan-Wigner transformation [JW2§| is the isomorphic mapping that al- 
lows the description of a fermionic system by the conventional model. It is per- 
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formed in the following way: 




(2.35) 



where the Pauli operators cr^ were previously introduced in Sec 12.11 If these op- 
erators satisfy the su(2) commutation relations (Eqs. 12.41 and 12.51) . the operators 
cj and c- obey the anticommutation relations of Eqs. 12.341 This is an isomorphic 
mapping between operator algebras and is independent of the Hamiltonian of the 
fermionic system to be simulated. 

Different Hamiltonians establish different connections (connectivity) between 
fermionic modes. Historically, Eqs. 12.351 and 12.361 correspond to lattices in one 
space dimension. Nevertheless, it is also vàlid for lattice systems in any dimen- 
sion, when the set of modes j is countable. In particular, the set of all ordered 
p-tuples of integers can be placed in one-to-one correspondence with the set of 
integers. For example, the simulation of a two dimensional fermionic lattice sys- 
tem can be done by re-mapping each mode (/,m) into a new set of modes as 
j = m + {l — l)Nx, where [/ = 1 • ■ ■ Ny] and [m = 1- ■ ■ N^] are integer numbers 
that refer to the position of a site in the lattice, and A^^, and Ny are the number of 
sites (modes) in the x and y direction, respectively. 

In order to compute physical properties of a fermionic system on a QC de- 
scribed by the conventional model, every step of the quantum simulation has to 
be expressed in terms of Pauli operators. For a (spinless) fermionic system with 

modes, a QC must contain, besides the ancilla qubit a, N qubits to represent 
the system. In the following, I describe how certain fermionic initial states can be 
prepared and how they can be evolved under a particular fermionic Hamiltonian 
evolution. 

Preparation of Initial Fermionic States 

Associated to each fermionic mode, there are two levels which correspond to the 
mode being empty or being occupied by a spinless fermion. The state-state map- 
ping is then trivial. Basically, the logical state \ lj) is associated to the jth mode if 
it is empty, and the logical state |0j) (up to a phase) is associated if the jth mode 
is occupied. In this way, the vacuum or no-fermion state |vac), which satisfies 
Cj|vac) = Vj, is mapped to the logical A^-qubit state II1I2 ■ ■ ■ Iat). 
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However, when simulating a fermionic system, more complex states need to 
be prepared. A general state of N^, fermions is a linear combination of Slater 
determinants (i.e., fermionic product states), 

L 

lV') = $^í?«l0a), (2.37) 
where the Slater determinants |0q,) are 

\<Pa) = l[c] \yac). (2.38) 

Due to the anticommutation relations of Eqs. 12.341 the fermionic operators satisfy 

c^ = -c]4íí^^J^ (2-39) 

implying that the Slater determinants \(pa) are antisymmetric wave functions under 
the permutation of an even number of fermions. 

Every state can be prepared on a QC made of qubits, by noticing that the 
quantum gate (i.e., unitary operator) 

= e^5('='"+^™) (2.40) 

creates a particle in the m-th mode when acting on the vacuum state. In other 
words, f/m|vac) = e'^^clJvac). Then, making use of the Jordan-Wigner transfor- 
mation (Eqs. I2.35l and l2.36l) . the operators Um in the spin language are 

771—1 

í|<T- n -<^i 

Um = e -1 . (2.41) 

The operators Um can easily be decomposed into elementary single and two-qubit 
gates as described in Sec. 12.1.21 The successive application of Ne similar uni- 
tary operators to the state |lil2 ■ ■ ■ Iat) generates the mapped state |0q,), up to an 
irrelevant global phase. 

The general fermionic state of Eq. 12.371 can be prepared by using L ancilla 
qubits, performing unitary controlled-í7m evolutions on the state of the ancillas, 
and finally, performing a measurement (projecting) on the ancillas. For example, 
if one is interested in preparing the state \^) = -^[\<Pi) + \<p2)], one needs to add 
an extra ancilla to the system. This ancilla is prepared in the state |+)a and a 
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controUed evolution to obtain the state ;^[|0)a ® |</>i) + |l)a ® |02)], is performed 
later. If the Hadamard gate is applied to the ancilla, this state evolves into 

10)3 ® + 102)] + |l)a ® ^Ul) + 102)]. (2.42) 

Therefore, the ancilla qubit is measured and projected, with probability 1/2, into 
|0)a or |l)a. If the former is obtained, the desired state is prepared. However, if 
the ancilla is projected into | l)a, the whole method needs to be applied again from 
the begining. 

In general, the probability of successful preparation of \ip) (Eq. 12.371) using 
this method is 1/L. Then, the order of L trials need to be performed before a 
successful preparation. A detailed description of this method can be found in 
Ref. IIOGKOll . 

Nevertheless, an important case consists of the preparation of Slater determi- 
nants (product state) in a different basis mode than the one given before: 

10/3) = 1^3^)- (2.43) 

i=i 

The fermionic operators d^'s are sometimes related to the operators through the 
foUowing canonical transformation 

_ ^a?^t^ (2.44) 

with = (c?í, ■ ■ ■ , c^lv)' "c^^ = (4) 4' ■ ■ ■ 5 cjv), and M being a x 
Hermitian matrix. (Sometimes the operators S are combinations of both, the 
creation and annihilation operators c\ and Cj.) 

Thouless's theorem states that one Slater determinant evolves into the other as 



|0/3) = í/|0a), (2.45) 

where the unitary fermionic operator 

U = e"^^^*'^ (2.46) 

can be written in terms of Pauli operators using the Jordan-Wigner transformation 
("Sec 12. 3. ÏT) . and can also be decomposed into elementary gates as described in 
Sec. l2T2l 
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In brief, the described fermionic product states can be prepared on a QC de- 
scribed by the conventional model, if the Jordan-Wigner transformation is per- 
fomied. Interestingly, the preparation can be done efficiently: the number of el- 
ementary single-qubit and two-qubit gates required scales polynomially with the 
system size N. In Chap.Hl I present another class of fermionic states that can also 
be prepared efficiently. 



Fermionic Evolutions 

The evolution of a quantum state is the second step in the realization of a QA. The 
goal is to decompose a genèric evolution into the elementary gates -R^(í9) and -R^j^fe (uj) 
(Sec. 12. ll i. Sometimes, the evolution step is associated to a Hermitian operator H 
which is, for example, the Hamiltonian H of the fermionic system to be simu- 
lated in terms of Pauli operators after Eqs. l2.35l and l2.36l have been performed. In 
this case, the corresponding evolution unitary operator is U{t) = e~*^* (i.e, the 
solution to the Schròdinger's evolution equation). 

In general, a fermionic Hamiltonian can be decomposed as H = K+V ,where 
K represents the kinetic energy of the fermions and V their potential energy. Usu- 
ally, [A', \/] 7^ and the decomposition of U (í) in terms of elementary gates is 
a complicated task. To avoid this difficulty this operator is approximated by, for 
example, using a first order Trotter decomposition IISuz93ll . That is, 

J\í 

Ü{t) = HÜiAt), (2.47) 

9=1 

Ü{At) = e™ = e*(^+^)^* = e'^^V^^* + (9(Aí2), (2.48) 

where K and V are the terms K and V in Pauli operators, respectively. Therefore, 
for At 0, Ü{At) ~ e'^^'e'^^K 

The potential energy V is usually a sum of commuting diagonal terms, and 
the decomposition of e*^^* into elementary gates is simple. However, the kinetic 
energy K is usually a sum of noncommuting hopping terms of the form CjCf, + c\cj 
(bilinear fermionic operators), and its decomposition is again approximated. A 
typical kinetic term e^^^í^'^'^^f'^i^^^ (j < k), when mapped onto the spin language 
gives 

e '=J+i = e '=J+i e . (2.49) 
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The decomposition of each term on the right hand side of Eq. l2.49l into elementary 
single and two-qubit gates was previously discussed in Sec. 12.1 .21 The amount of 
elementary gates required depends on the distance | j — A; | , and scales polynomially 
with that distance. Moreover, since H represents a physical system, it is a linear 
combination of a polynomially large (with N) amount of fermionic operators. 
Then, U (t) can be performed efficiently by applying a polynomially large amount 

of elementary gates. In the same way, the unitary operation U = e~* of 
Eq. l2.46l can also be efficiently implemented. 

Obviously, the accuracy of approximating U (í) using the Trotter decomposi- 
tion increases as Aí decreases. Then, a large amount of gates might be required to 
perform the desired evolution with small errors. To overcome this problem, one 
could use a Trotter approximation of higher order in At IISuz93l . All these approx- 
imation methods do not destroy the efficiency of the QA. Moreover, the evolution 
step induced by fermionic physical Hamiltonians with higher order products of 
creation and annihilation operators can also be efficiently implemented using the 
same techniques. 

2.3.2 Simulations of Anyonic Systems 

The concepts described in Sec. 12.3. li can be extended to other and more general 
particle statistics, namely hard-core anyons [BOOl |. These are partides that also 
obey the Pauli's exclusion principle: At most one (spinless) anyon can occupy a 
single mode. Assigned to each mode of the lattice are the creation and annihilation 
anyonic operators aj and a^, respectively. Their commutation relations are given 



where Uji = cLjidji is the number operator, [A,!?]^ = AB — e^^BA, and 6 is 
the statistical angle. In particular, 9 = n mod(27r) represents canonical spinless 
fermions, while 9 = mod(27r) represents hard-core bosons. 

In order to simulate this problem on a QC described by the conventional 
model, the foUowing isomorphic mapping between algebras can be performed: 



by (j < f) 




[a],a],]e = , 
5,y(l-(e-^ + l)n,0, 



(2.50) 
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where the Pauli operators were introduced in Sec. 12.1.11 Since they satisfy 
the commutation relations of Eqs. 12.41 and l2. 51 the commutation relations for the 
anyonic operators (Eqs. 12.501) are satisfied. 

As in the fermionic case (Sec. 12.3.11) . an anyonic evolution operator can be 
written in terms of Pauli operators using Eq. 12.511 and can be decomposed into 
single and two-qubit elementary gates. Therefore, the same procedure described 
in the previous section can be foUowed. 

Anyon statistics have fermion and hard-core boson statistics as limiting cases, 
satisfying always the Pauli's exclusion principle. In the next section this hard-core 
condition is relaxed and the important case of canonical bosons is considered. 

2.3.3 Simulations of Bosonic Systems 

Quantum computation is based on the manipulation of quantum systems that pos- 
sess a finite number of degrees of freedom (e.g., qubits). From this point of 
view, the simulation of bosonic systems appears to be impossible, since the non- 
existence of an exclusion principle implies that the Hilbert space used to represent 
bosonic quantum states on a lattice is infinite-dimensional; that is, there is no limit 
to the number of bosons that can occupy a given mode j. However, sometimes 
it is necessary to simulate and study properties such that the use of the complete 
Hilbert space is unnecessary, and only a finite sub-basis of states is sufficient. This 
is the case for A^-mode (e.g., A^ sites) lattice systems with interactions given by 
the boson-preserving Hamiltonian 

N 

H=Y^ «i/ b]hj, + (3jj. üjfij,, (2.52) 
i,i'=i 

where the operators bj (bj) create (destroy) a boson at site j, and fij = b'jbj is the 
number operator; that is 

b]\ni,n2,- ■ ■ ,nj,- ■ ■ ,nN) = v^n~TI ng, ■ • • , + 1, ■ • ■ , tíat), 
bj\ni,n2,- ■ ■ ,nj,- ■ ■ ^un) = ^nj |ni, n2, ■ ■ ■ , - 1, ■ ■ ■ , Wat), 
nj\ni,n2, - ■ ■ ,nj, - ■ ■ ,nN) = rij «2, ■ ■ ■ , n^, ■ ■ ■ , «Ar), (2.53) 



n 



'e'^ + 1 



(2.51) 
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n 



, njv) represents a quantum state with 



where the bosonic state |ni, n2, ■ ■ 
rij bosons in the j-th mode (site). 

The space dimension of the lattice is encoded in the parameters ajji and Pjji 
of the Hamiltonian. Since H contains pairs of creation and annihilation operators, 
the total number of bosons Np in the system is preserved and the idea is to work in 
this finite sub-basis of states (where the dimension of the associated Hilbert space 
depends on the magnitude of Np). 

The corresponding bosonic commutation relations (in an infinite-dimensional 
Hilbert space) are IICDG98J 



(2.54) 



However, if the operators bj are restricted to the finite basis of states represented 
by {|?T,i, 71.2, ■ ■ ■ , npf) with max(?7,j) = Np}, that is, Np is the maximum number 
of bosons per site, they acquire the foUowing matrix representation (see Eqs. 12.531) 



> ^ ( 
jti^ factor 



(2.55) 



where the symbol ® indicates the usual tensorial product between matrices, and 
the {Np + 1) X (Np + 1) dimensional matrices 11 and are given by 
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(2.56) 



It is important to note that in this finite basis, the commutation relations of the h\ 
differ from the Standard bosonic ones (Eq. 12.541) UB 0021 



(2.57) 



and clearly (0])^^+^ = 0. 

Since the goal is to simulate the bosonic system on a QC described by the con- 
ventional model, a corresponding mapping between both operator algebras must 
be given. Nevertheless, Eqs. I2.57l implv that the linear span of the operators fe] and 
hj is not closed under the commutator, and a mapping between the bosonic oper- 
ators and the Pauli operators like the Jordan-Wigner transformation (Sec. 12.3.11) 
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is not possible. Therefore, such isomorphic mapping needs to be found by first 
mapping quantum bosonic states onto quantum logical states in the conventional 
model (i.e., a Hilbert space mapping). 

The idea is to start by considering only the jth mode in the chain. Since 
this mode can be occupied with at most Np bosons, it is possible to associate an 
{Np + l)-qubit quantum state to each particle number state, in the foUowing way: 



|0), - 


^ |Oolll2·· 








^ lloOils- 


■ ■ ^Np)j 




|2), - 


^ lloliOs- 


■ ■ lAfp)i 


(2.58) 



\Np), ^ |lolil2···0,Vp)i 

where \n)j denotes a quantum state with n bosons in jth mode. Therefore, N(Np+ 
1) qubits for the simulation (where is the total number of modes) are needed. 
An example of this mapping for a quantum state with 7 bosons in a chain of 5 
sites, where the maximum number of bosons per site is Np = 3, is shown in 

Fig.m 

By definition (see Eqs. 12.531 12.551 and 12.561) \n)j = y/n + 1 \n + and 
this operator in the conventional model maps to 

h] ^h]=Y, V^ayal^^'\ (2.59) 

n=0 

where a pair (n, j) refers to the nih. qubit in the chain of qubits representing the 
jth bosonic mode. The Pauli creation and annihilation operators a\ were previ- 
ously defined in Sec. l2.1[ The operator òj acts then on the {Np + l)-qubit chain 
representing the jth bosonic mode as 

h]\lo ■ ■ ■ l„_iO„l„+i ■ • - Inp)] = Vn + 1 |lo ■ ■ ■ lnOn+lln+2 " " ■'^Np)j, (2.60) 

SO its matrix representation in this basis is analogous to the matrix representation 
of bj in the basis of bosonic states. 

Similarly, the number operator is mapped as 

% ^nj = J2^ ' (2.61) 

n=0 
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\(t>a) = iim>i <s> \im)2 ® i;m>3 <s> \unu ® ituDs 



Figure 2.9: Mapping of the bosonic state {(pa), of a chain with 5 sites and 7 bosons 
{Np = 3), into a four-spin-1/2 (or four-qubit) state. The convention is |tj) = |0j) 
and lij) = \ 

so its action over the corresponding logical states is 

^j|lo ■ ■ ■ In-lOnln+l ' ' ' '^Np)j = n \ lo ■ ■ ■ l„0„+il„+2 " ' ' '^Np)j- (2.62) 

Since the commutator [òj, J2n=o '^"'"'l ^ ^ operators òj (6j) always keep states 
within the same subspace. 

The Hamiltonian of Eq. I2.52l in terms of Pauli operators is then 

N 

^= ^rï + f^n' ^j^j'^ (2-63) 

where the operators (bj) are given by Eqs. 12.591 and üj by Eq. 12.611 In this 
way, physical properties of the bosonic system such as the spectrum of H can 
be obtained using a QC made of qubits. The same methods can be used when 
simulating any other type of boson-preserving quantum system. 
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Preparation of Initíal Bosonic States 

As in the fermionic case, the most general bosonic state of an A^-mode lattice sys- 
tem with a maximum of Np bosons per site can be written as a linear combination 
of bosonic product states like 

= Kiblriblr ■ ■ ■ (èjv)"" |vac), (2.64) 

where K is a normalization factor, uj is the number of bosons at site j (max{nj) = 
Np), and |vac) is vacuum or no-boson state, that is, òj|vac) = Vj. 

Using the mapping described in Eq. 12.581 the vacuum state in the conventional 
model maps as 

|vac) ^ lOoli ■ ■ ■ l7Vp)i ® ■ ■ ■ ® lOoli ■ ■ ■ 1np)n, (2.65) 

and the product state of Eq. 12. 641 maps as 

\<Pa) = llo-'-On, ■■■Up)i®···® |lo···0„^···U^)^ (2.66) 

(see Fig. l2.9l for an example). 

The preparation of the mapped bosonic state |0q) on a QC made of qubits 
is then performed by flipping the states of the corresponding qubits from a fuUy 
polarized state (i.e., the logical state with all qubits in |1)), using for example the 
flip operations cr" -' . Nevertheless, more general bosonic states like 

L 

IV') = (2.67) 

can also be realized as in the fermionic case. Again, the idea is to add L ancillas 
(extra qubits), perform controUed evolutions on their states, and finally perform 
measurements on the state of the ancillas. The state is successfuUy prepared with 
probability 1 /L llÜÜKOlJ . 

Bosonic Evolutions 

Again, the idea is to represent certain bosonic unitary evolution operator U (í) = 
e~*^*, where H is some boson-preserving Hermitian operator such as the Hamilto- 
nian of the system to be simulated (Eq. l2.52l) . in terms of Pauli operators ("Eq. 12.631) . 
Usually, a first order Trotter approximation IISuz93l also needs to be performed to 
separate those terms in H that do not commute. 
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In general, H = K + V (Eq. l2.52l) . where K is a kinetic term and V a potential 
term. The kinetic term is a linear combination of terms like 6|.ò; + ò^ò^.. Therefore, 
a single-step evolution operator 6*'^^*'''+''?^'=^^* is mapped onto Pauli operators as 



exp 



Np~l 



n,n'=0 



.1 n'+l,l\ I f „n,k n+l,k „n,k n+l,k\ /n' .1 n' +1,1 ^n' ,1 n' 
^y '^V ) + K^y -^y^x A^x '^y ' ^y ^x )\ 



(2.68) 



where = Aí/8 and Np is the maximal number of bosons per site. The terms 
in the exponent of Eq. l2.68l commute with each other, so the decomposition into 
elementary gates can be done using the methods described in Sec. 12.1.21 As an 
example, consider a system of two sites with maximal one boson per site {Np = 
1). Thus, 2(1 + 1) = 4 qubits are needed for the simulation, and Eq. I2.59l implies 
that b\ = cr°'V^'^ and b\ = cr^'^cr]l^. The mapped bosonic operator e*'^^i^2+''2^i)^* 
in terms of Pauli operators is 



exp(z^a°'V^'V°'V^'^) X exp(z^a°'V^'V;'VJ'^) x (2.69) 
expíz^aJ'VJ'Vfal'^) x explz^aJVJ'VfaJ'^) x expíz^aJV^'VJV^'^) x 
exp(-z^aJV^'V°'VJ'^) X exp(-^^a°'V;'VJV^'^) x exp(^^a°'V;'V°'V;'^), 



and the decomposition of each of the terms in Eq. 12.691 in terms of single and 
two-qubit elementary gates can be done, again, using the methods described in 
Sec. 12.1.21 An example of the decomposition of the term exp (f 0"°' V.^' V°'^cr^.'^í) , 
where the qubits were relabeled as (n, j) = n + 2j — 1 (e.g., (0, 1) — > 1) is shown 
in Fig.lTÏÜl 

Contrary to the fermionic case, the number of elementary operations involved 
in the decomposition is not related to the distance between sites, \k~l\. Neverthe- 
less, a physical bosonic operator H, such as the Hamiltonian of Eq. 12.521 involves 
a polynomially large number (with respect to N) of bosonic terms. Therefore, the 
corresponding evolution U (t) = e~*^* can be efficiently performed on a QC by 
applying a polynomially large number of elementary single and two-qubit gates. 
Again, when using approximate methods like the Trotter decomposition, the num- 
ber of operations needed increases with the desired accuracy. However, such an 
approximation does not destroy the efficiency of the simulation. 
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Figure 2.10: Decomposition of the unitary operator U{t) = e's'^^'^y^'y'^^ into ele- 
mentary gates as described in Sec. 12.1.21 The labeling convention is (0, 1) = 1, 
(1,1)^2, (0,2) = 3, and (1,2) =4. 

2.4 Applications: The 2D fermionic Hubbard model 

To clarify the methods described previously, in this section I present, as an ex- 
ample, the QS of the finite two-dimensional fermionic Hubbard model by using 
a CC that imitates a QC; that is, a quantum simulator. Since this is a classical 
simulation, the CC must keep track of an exponentially large number of quantum 
States, associated with the Hilbert space of the quantum system. Nevertheless, this 
simulation provides a good example for understanding the advantages and results 
that can be obtained when using a real QC. 

The physical system to be simulated consists of a rectangular lattice (Fig. 12.1 iT) . 
with Nx X Ny sites, where spin-1/2 fermions hop from site to site under the inter- 
action Hamiltonian 

(2.70) 

where the operators c|. -y^ (c(ijy,a) create (annihilate) a fermion with spin compo- 
nent denoted by a (=t or J,), at the site located at {x = i,y = j). and ty are the 
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Figure 2.1 1: Two-dimensional lattice in the Hubbard model. 

hopping terms in the x and y directions, respectively, and n{i^j).„ = c'(^i jy^c^- 
is the corresponding number operator. {H.C. denotes the Hermitian conjugate). 
Periòdic boundary conditions (PBC) are assumed: [{i,j);a] == [(i + Nx,j);a] = 
[{t,j + Nyy,a]. 

To use the QA for the obtention of the spectrum of H, described in Sec. 12.2.21 
it is necessary first to map the fermionic operators onto Pauli operators using, for 
example, the Jordan-Wigner transformation (Sec. 12.3.11 ). Considering that these 
are spin-1/2 fermions, the QA needs 2(N,j. x Ny) qubits to represent the system 
(qubits system), plus the ancilla qubit. 

Assuming that one is mainly interested in obtaining the lowest energy of 
Eq. 12.701 the prepared initial state should be the ground state of Eq. 12.701 How- 
ever, since no algebraic methods exist to exactly diagonalize Eq. l2.70l for large A^^, 
and Ny, such a state is not known, and therefore impossible to prepare. Neverthe- 
less, the ground state of the associated mean-field Hamiltonian 

Hmf = -Y1 [t-i^,,)■,aHi+lJy,- + + ^•^·] + 

^ Xl[('^(íj);T)^(íJ)a + n(i,jy,^{n^i,j)-i) - {n(ijy^) {n^ijyi)] , (2.71) 



is known to be a fermionic product state (Slater determinant) and its corre- 
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sponding state in the conventional model can be efficiently prepared by using the 
methods described in Sec. 12.3.11 that is, it can be prepared by applying a polyno- 
mially large (with respect to N^^Ny) set of elementary gates to the fuUy polarized 
State. For finite small lattices, l^/j) is a good approximation to the ground state of 
Eq. 12. 701 and can be used to obtain the ground state energy. 

The second step of the QA is to apply the unitary operator U (t) = e'^^^l"^ 
using single and two-qubit gates (see Fig. l2.8l for Q = H)), where, in this case, 
H is the Hamiltonian of Eq. 12.701 in terms of Pauli operators, t is a real (fixed) 
parameter, and al is the Pauli operator associated with the ancilla qubit. Since H 
is a linear combination of non-commuting terms (Eq. l2.70t . the operator U (t) = 
Ylj U ( Aí) can be approximated by using, for example, the first order Trotter de- 
composition IISuz93l . That is, H = K^+Ki + V, where denotes the kinetic 
energy associated with the fermions of spin a and V denotes the potential energy. 
Then, 

where K„ and V are the corresponding terms in Pauli operators. Also, the same 
approximation can be used to decompose each term e^'^'^'^'^*/^. Such approxi- 
mation leads to operators that can easily be decomposed in terms of elementary 
gates by using the methods described in Sec. 12.1.21 

The energy spectrum of the Hubbard model for a 4 x 2 lattice is shown in 
Fig. 12.121 It has been obtained after running the classical simulation for many 
different vàlues of t, and performing the Fourier transform on the data. The peaks 
show the eigenvalues and the results are compared to those of an exact diagonal- 
ization method. 

The algorithm described allows one to easily obtain the energy spectra of a 
finite small lattice. Nevertheless, the spectra of a large lattice cannot be efficiently 
obtained using the same methods. The problem is that the ground state of the 
mean-field approximation differs more from the ground state of Eq. 12.701 as the 
system increases, and the QA needs to be performed an exponentially large num- 
ber of times (with respect to A^^, x Ny) to obtain the desired property IISOG02I . 

2.5 Quantum Algorithms: Efficiency and Errors 

A QA for a physical simulation is considered efficient if the total number of opera- 
tions involved for the initial state preparation, the evolution, and the measurement 
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Figure 2.12: Energy spectrum of the Hubbard model obtained simulating the QA 
of Fig. lTHl on a CC. The lattice has 4x2 sites (which requires 16-1-1 qubits). Here, 
tx = iy = and W = 4. The time steps used in the Trotter approximation to pre- 
pare the initial state and apply the evolution are Aíi = AÍ2 = 0.05, respectively. 
The numbers in brackets are the results obtained from the exact diagonalization 
using the Lanczos method. 
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process, scales at most polynomially with the system size and with 1/e, where e 
is the maximal tolerable error in the measurement of a relevant property. 

While the decomposition of the operator U (í) = e^*^* can be done efficiently 
(e.g., when using the Trotter approximation) if is a physical operator, the prepa- 
ration of a general initial state could be inefficient. Such inefficiency would arise, 
for example, if the state {ip) defined in Eq. 12.371 or Eq. 12.671 is a linear combina- 
tion of an exponentially large number of elementary product states. In this case, 
L ~ , with the number of modes in the system and a; > 1, so an exponen- 
tially large number of trials need to be performed before successful preparation. 
However, if L < poly(A^), it can be prepared efficiently. This construction gen- 
eralizes to more general coherent states (see Chap.El). 

The three main reasons for the existence of errors e in the outcome of the 
quantum computation are gate imperfections, the use of the Trotter approxima- 
tion in the evolution operator, and the statistics in measuring the polarization 
of the ancilla qubit (Sec. 12.2.21) . Gate imperfections are very common in quan- 
tum Information because, contrary to classical Information, quantum gates are 
usually dominated by a continuous parameter. This problem can be solved by 
using quantum error correction methods and fault tolerant quantum computa- 
tion rtSte96. ■Kit97[ INCOOL IGot97l . According to the accuracy threshold theorem, 
provided that the physical gates have sufficiently low error, it is possible to quan- 
tum compute accurately in an efficient way. 

The type of error introduced by the discretization of the evolution operator 
U (í) (e.g., by using the Trotter decomposition or other approximations) is very 
similar to the error obtained when performing a classical simulation, such as when 
using Monte Cario methods. This error can be estimated by a detailed analysis of 
the discretization and can also be arbitrarily reduced in an efficient way. 

Finally, when using the QAs described in Sec. 12.2.21 the step corresponding 
to the measurement process can also be performed efficiently because it only in- 
volves the measurement of a single (ancilla) qubit, regardless of the number of 
qubits needed for the simulation. Nevertheless, repeatedly many same-simulations 
need to be performed to get an accurate value of such measurement. This is an in- 
herent property of the quantum mechanichs where a single measurement projects 
the quantum state (Sec. l2.1.TT l and does not give sufficient Information. If the rel- 
evant signal at the end of the quantum computation is small, like when obtaining 
the spectra of the two-dimensional Hubbard model on a large lattice (Sec. 12.41) . it 
is necessary to run the algorithm a larger number of times and the efficiency can 
be destroyed. 

In brief, a QC can only be more efficient than a CC when simulating quan- 
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tum physical systems if the three main steps of the corresponding QA can be 
performed efficiently. For example, the evaluation of certain correlation functions 
over a quantum state that can be easily prepared, can be efficiently done with a 
QC. In general (i.e., for non-integrable Hamiltonians), there is no known way to 
evaluate such correlation efficiently with a CC JSÜÜÜllSÜKÜa . 

2.6 Experimental Implementations of Quantum Al- 
gorithms 

In this chapter, I have shown that if a large QC existed today, some simulations 
of quantum systems could be performed more efficiently on it than on a CC. Nev- 
ertheless, I did not discuss how the corresponding QAs could be experimentally 
implemented. Although numerous proposals for implementing quantum Informa- 
tion processors (QIPs) are found in the literature IICZ95L lÜLKÜÜl IKLMCT . only 
few of them have been successfuUy implemented to process more than one qubit. 
In particular, liquid-state NMR devices allow one to simulate several systems by 
manipulating, nowadays, up to ten qubits IIRBC04I . 

The physical implementation of a large scale QC still remains one of the most 
important challenges for today's physicists. The problem is a QC should be de- 
signed such that the interaction between its constituents and the environment is 
small enough to keep coherence of the quantum state. But if such interaction is 
too small, the manipulation and control processes using external sources becomes 
impracticable. For this reason, quantum decoherence has been one of the most 
important subjects of study during the last decade. In general, decoherence phe- 
nomena is hard to predict due to the infinite degrees of freedom associated with 
the environment. Nevertheless, a QC is reliable whenever the time required to 
perform a certain task is much smaller than the corresponding decoherence time. 

In this section, I describe the experimental setting of a liquid-state NMR QIP 
and show how such devices can be used to execute the QAs described previously. 
Later on, I will show the experimental NMR simulation of a particular fermionic 
system, where some correlation functions and energy spectra have been obtained. 

2.6.1 Liquid-State NMR Quantum Information Processor 

Liquid-state NMR methods allow one to physically implement a slightly different 
version of the conventional model of quantum computation, with respect to the 
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initial State preparation and the measurement process. In this set-up the quantum 
register is represented by the average state of the nuclear spin-1/2 of an ensemble 
of identical molècules. Each nuclear spin is a two-level physical system and can 
then be considered a possible qubit. Thus, the idea is to perform single and two- 
qubit elementary gates by extemal radio-frequency (rf) pulses that interact with 
the nuclear spin state. In the foUowing, I present a bàsic analysis about how these 
processors can be used as possible QCs. 

In a liquid NMR setting, the molècules are placed in a strong magnètic field 
B(z) ^ 10 T, so that the spin of the j-th nucleus of a single molecule precesses 
at its Larmor frequency uj (Fig. 12.131) . In the frame rotating with the jth spin, its 
qubit state can then be rotated by sending rf pulses in the XY plane at the resonant 
frequency Uf. ^ Uj. If the duration of this pulse is 5t, the corresponding evolution 
operator in the rotating frame is IILKC02I 



where A is the amplitude of the RF-pulse and (p is its phase (i.e., orientation) in 
the XY plane {h = 1). Then, one can induce single spin rotations^ around any 
axis in that plane by adjusting 5t and ip. 

Single-qubit rotations around the z axis can be implemented with no exper- 
imental imperfection or physical duration simply by changing the phase of the 
abstract rotating frame with which one is working. One has then to keep track of 
all these phase changes with respect to a reference phase associated with the spec- 
trometer. Nevertheless, these phase tracking calculations are linear with respect 
to the number of pulses and spins, and can be efficiently done on a classical com- 
puter. Together with the rotations around any axis in the XY plane, the z rotations 
can generate any single qubit rotation on the Bloch sphere. 

Two-qubit gates, like the Ising gate i?^j^fe(a;) (Sec. I2.1.n ). can be performed 
by taking advantage of the spin-spin interactions (i.e., nuclei interaction) present 
in the molecule, and then achieve universal control. To first order in perturbation, 
this interaction, named the J-coupling, has the form 



where j, k denote the corresponding pair of qubits and Jjk is their coupling strength 
Under typical NMR operating conditions, these interaction terms are small enough 



—iA{cos{ifi)crx+sin{ip)a-y)St 



(2.73) 



(2.74) 



^One actually is restricted to 90 and 180 degrees rotations for experimental calibration issues. 
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Figure 2.13: Bloch's sphere representation of a single nuclear spin-1/2 state pre- 
cessing around the quantization axis determined by the extemal magnètic field B. 
The precession frequency is given by Uj = HjB, with i^ij the magnètic moment 
of the jth nucleus. Due to the chemical environment, each nucleus precesses at a 
different Lannor frequency i^j. 

to be neglected when performing single-qubit rotations with rf pulses of short du- 
ration. Nevertheless, between two pulses they are driving the evolution of the 
system. By cleverly designing a pulse sequence, i.e., a succession of pulses and 
free evolution periods, one can easily apply two-qubit gates on the state of the 
system. Indeed, the so-called refocusing techniques' principle consists of per- 
forming an arbitrary Ising gate by flipping one of the coupled spins (vr -pulse), 
as shown in Fig 12.141 The interaction evolutions before and after the refocusing 
pulse compensate, leading to the effective evolution 



where the effective coupling strength a = Jjk{Sti — 5t2) is being determined by 
the difference between the durations 5ti and 5t2. 

Although the physics on a single molecule has been analyzed, liquid-state 
NMR uses an ensemble of about 10^'^ molècules in a solution maintained at room 
temperature 300-fí). For typical vàlues of the magnètic field, this thermal state 
is extremely mixed. Clearly, this is not the usual state in which one initializes 
a quantum computation since qubits are nearly randomly mixed. Nevertheless, 





(2.75) 
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Figure 2.14: Circuit representation for the refocusing scheme to control J cou- 
plings. The Ising-like coupling Jjk between spins can be controUed by perform- 
ing flips on one of the spins at times ti = Ati and Í2 = + AÍ2, respectively. 
The effective coupling is a = ai — «2 = Jjki^h — ^h), and vanishes when 
Aíi = Aí2. 

known NMR methods IILKC02I can be used to prepare the so-called pseudo-pure 
State (ppp)·^ 

Ppp = -^-^jv"^^ + ^Ppu'·e' (2-76) 

where ï is the identity operator, pp^re is a density operator that describes a pure 
State, and e is a small real constant (i.e., e decays exponentially with the number 
of àtoms in the solution due to the Boltzmann's distribution). 

Under the action of any unitary evolution U, this state evolves as 

plf = Up,,U^ = ^^^11 + í/epp,,eí/t. (2.77) 

The first term in Eq. I2.77l did not change because the identity operator is invariant 
under any unitary transformation. Therefore, performing quantum computation 

^Even though efficient techniques to prepare a pseudo-pure state exist in theory ISV98I . they 
are very hard to implement in practice, and one instead uses non-efficient methods that suffer an 
exponential decay of the observed signal with respect to the number of qubits in the pseudo-pure 
state. 
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on the ensemble is equivalent to performing quantum computation over the initial 
State represented only by Ppure- 

At the end of the computation, the orthogonal components of the sample po- 
larization in the XY plane, M,. = Jr{pl"^-' Y.L· 4). and My = Jr{plf Y.L· 4) 
are measured (Eq. 12.151) . Note that the invariant component of Ppp^' does not con- 
tribute to the signal since Tr(]lcr^ j^) = 0. Because the polarization of each single 
spin, Mi = Tr(ppp^'(T^) and = Tr(ppp^'(j^), precesses at its own Larmor fre- 
quency Vj, a Fourier transformation of the temporal recording (called FID, for 
Free Induction Decay) of the total magnetization needs to be performed. By do- 
ing so, one obtains the expectation value of the polarization of each spin (averaged 
over all molècules in the sample). 

Summarizing, a liquid-state NMR setting allows one to initialize a register of 
qubits in a pseudo-pure state, apply any unitary transformation to this state by 
sending controUed rf pulses or by leaving free interaction periods, and measure 
the expectation value of some quantum observables (i.e., the spin polarization). 
Hence, these systems can be used as QIPs. 

2.7 Applications: The Fano-Anderson Model 

I now present the experimental QS of the fermionic one-dimensional (ID) Fano- 
Anderson model using a liquid-state NMR IINSO05L by manipulating the state of 
the spin nuclei as described in Sec. 12. 6. li Such simulation shows then how reliable 
these experimental methods are and how well the elementary gates (Sec. l2.1.Tl) can 
be implemented using NMR techniques. 

The ID fermionic Fano-Anderson model consists of an n-sites ring with an im- 
purity in the center (Fig. 12.151) . where spinless fermions can hop between nearest- 
neighbors sites with hopping matrix element (overlap integral) r, or between a site 
and the impurity with matrix element V/ ^Jn. Taking the single-particle energy of 
a fermion in the impurity to be e, and considering the translational invariance of 
the system, the Fano-Anderson Hamiltonian can be written in the wave vector 
representation as lOGKOll 

^ = E ^^A^'h + ^^^^ + ^(4o^ + ^^CfcJ, (2.78) 

where the fermionic operators c|,^ and 1)^ (c^^ and ò) create (destroy) a spinless 
fermion in the conduction mode ki and in the impurity, respectively. Here, the 
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Figure 2.15: Fermionic Fano-Anderson model. Fermions can hop between nearest 
neighbor sites (exterior circles) and between a site and the impurity (centered 
circle), with hopping matrix elements r and V/ ^/n, respectively. The energy of 
localization in the impurity is e. 

wave vectors (modes) are /c; = ^ (/ = [0, .., — 1]) and the energies per mode 
are = — 2r cos ki. 

In this form, the Hamiltonian in Eq. l2.78l is almost diagonal and can be exactly 
solved: There are no interactions between fermions in different modes ku except 
for the mode k^, which interacts with the impurity. Therefore, the relevant physics 
comes from this latter interaction, and its spectrum can be exactly obtained by 
diagonalizing a 2 x 2 Hermitian matrix, regardless of n and the number of fermions 
in the ring, N^. Nevertheless, its simulation in a liquid-state NMR QIP is the 
first step in QSs of quantum many-body problems and consti tutes a proof of the 
principies described throughout this thesis. 

In order to successfuUy simulate this system in a liquid-state NMR QIP, the 
fermionic operators need to be mapped onto the Pauli operators (Sec. l2.3.TTl . This 
is done by using the foUowing Jordan-Wigner transformation 



b = a\ 



fet = a 



.1 
+ 




(2.79) 
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Figure 2.16: Mapping of the fermionic product state b'^cl^cl^ |vac) into a five-qubit 
State, using the Jordan-Wigner transformation. The convention is = |0j) 
(fiUed) and = \ (empty). 

In this language, a logical state \0j) (with |0) = |t) in the usual spin-1/2 notation) 
corresponds to having a spinless femiion in either the impurity, if j = 1, or in the 
mode kj^2, otherwise (Fig. l2.16T ). (Again, the fermionic vacuum state |vac) maps 
onto |vac) = jlilg ■ ■ ■ l„+i).) 

The algorithms described in Sec. 12.2.21 can be used, for example, to evaluate 
the probability amplitude of having a fermion in mode ko at time t, if initially 
(í = 0) the quantum state is the Fermi sea state with Ne fermions; that is, |FS) = 

iVe-l 

Yl cl |vac). This probability is given by the modulus square of the foUowing 

1=0 ' 

dynamical correlation function: 



ko 



impurity 



G(t) = (FS|ò(í)6t(0)|FS), 



(2.80) 
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where b{t) = T^b{0)T, T = e~*^* is the time evolution operator, and 6^(0) = . 
Basically, G{t) is the overlap between the quantum state 0^^(0)1 FS), which does 
not evolve, and the state ò^(í)|FS), which does not vanish unless the evolved state 
T|FS) already contains a fermion in the impurity site {(b'^(t))'^ = (61^(0))^ = 
0), i.e., contains the fermion which initially was in the ko mode. In temis of 
Pauli operators (see Eq. 12.791) . this correlation function reduces to a two-qubit 
problem IIÜÜKUll : 

= (0|f Vf , (2.81) 

where T = e~'^* is an evolution operator arising from the interaction terms in Eq. 
ITtHI with 

H = l<yl + + ^{alal + alal) , (2.82) 

and 10) = II1O2) in the logical basis (i.e., the initial state with one fermion in the 
ko mode). 

In order to use the quantum circuit depicted in Fig. 12.71 all operators in Eq. 
I2.81l must be unitary. Because of the symmetries of H, such as the global 'k/2-z 
rotation that maps (cr^, cr^) {a^y, leaving the state |0) invariant (up to a 

phase factor), then (^TV^T a,^|0) = (0|f tf^if (j^l^) = and (0|f tcrif a^l^) = 
(0|TV^Tcr^|0). Therefore, Eq. I2.81l can be written in terms of unitary operators 
as 

G{t) = (0|e^^V^e-^^V^.|0). (2.83) 

Figure ITT71 shows the quantum circuit used to obtain G{t). It is derived from 
Fig. l2.7l bv making the foUowing identifications: T — > e~*^*. Ai al, and Bj — > 
al- The corresponding controUed operations C-A and C-B transform into the 
well-known controUed-not (CNOT) gates (Sec. 12. 1.31 . All the unitary operations 
appearing in Fig. l2.17l were decomposed into elementary NMR gates (single qubit 
rotations and Ising interactions). In particular, the decomposition of e~*^* can be 
found in Ref. LQGKOli . obtaining 

e"*^* = [/e-^^i'^'*e~*^2'"'*í/^ , (2.84) 

where Ai(2) = T v^A^TV^), with E = and A = The unitary 

operator U is decomposed as (Fig. 12.171 ) 

U = e^^''''^e~^^''ye~^i''^-''''-é^"yé^''^''e~^^''^'^e~^^ (2.85) 



with the parameter 6 satisfying cos 6 = 1/ Vl + 5^, and 5 = ( A + VA^ + V^) /V . 
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Figure 2.17: Quantum circuit for the evaluation of G{t) (Eq. 12.801) in terms of 
elementary gates directly able to be implemented with liquid-state NMR methods. 
The controlled operations 6'^^= and Al°^= correspond to the operations C-B and C- 
A of Sec. 12.2.21 respectively. 
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The CNOT gates C-A and C-B can also be decomposed into elementary gates, 
as explained in Sec. l2.1.2l Therefore, G{t) can be obtained using an NMR QIP by 
applying the appropriate rf pulses (Sec. l2.6.TT) . Remarkably, only three qubits are 
required for the simulation (Fig. l2.17T) : The ancilla qubit a, one qubit representing 
the impurity site (qubit- 1), and one qubit representing the ko mode (qubit-2). 

Similarly, the algorithm depicted in Fig. 12.81 can be used if interested in ob- 
taining the spectrum of the Hamiltonian H of Eq. 12.781 replacing Q ^ H. In 
particular, when n = 1 (one site plus the impurity), Eq. 12.781 in terms of Pauli 
operators reduces to H = lí|^ + with H defined in Eq. 12.821 In this case, the 
two eigenvalues Xi (i = 1,2) of the one-particle subspace can be extracted from 
the correlation function 

S{t) = {(j)\e-'^'\(f)) = e-^(^+='=o)*(0|e-^^*|0), (2.86) 

which can be obtained by measuring the polarization of the ancilla qubit after the 
quantum circuit shown in Fig. l2.8l has been applied. Since |0) = II1O2) is not an 
eigenstate of H, it has a non-zero overlap with the two one-particle eigenstates, 
called |lPí) (see AppendixlBl). 

Again, the operator e*^*^^*/^ (Fig. 12. 81) needs to be decomposed into elementary 
gates for its implementation in an NMR QIP. Noticing that [al^H] = [cr|, í/] = 0, 
then 

where the unitary operator U is decomposed as in Eq. 12.851 Figure l2T8l shows the 
corresponding circuit in terms of elementary gates. Again, qubits 1 and 2 represent 
the impurity site and the mode, respectively. a denotes the ancilla qubit. Since 
the idea is to perform a DFT on the results obtained from the measurement (see 
Appendixini), this circuit needs to be applied for several vàlues of the parameter t 
(Sec. irT3l . 

2.7.1 Experimental Protocol and Results 

The NMR setting used for the evaluation of G{t) and S{t) is based on an ensem- 
ble solution of trans-crotonic acid and methanol dissolved in acetone (Fig. l2.19t . 
This molecule can be used as a seven-qubit register, where methanol is used to 
perform rf-power selection and accurately calibrate the rf pulses. In this way, the 
decoherence time (Tg*) is in the range from several hundreds of milliseconds to 
more than a second, allowing one to perform around 1000 single-qubit gates and 
around 100 two-qubit (Ising) gates. 
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Figure2.18: Quantum circuit for the evaluation of ^(í) rEq. l2.86ll . The parameters 
Al and A2 are defined in Sec. 12.71 and a = (e + £ko)/'^· The decomposition of the 
operator U in NMR gates is shown in Fig. 12.171 
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Figure 2.19: The transcrotonic acid molecule is a seven-qubit register. The methyl 
group is used as a single qubit HKLMOOj . The table shows the vàlues of the Chem- 
ical shifts (main diagonal) and the J couplings (off-diagonal) between every pair 
of nuclei or qubits, in hertz. 
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The relevant nuclei of the molecule are denoted as Ci, C2, C3, and C4, corre- 
sponding to the Carbon àtoms, Hi and H2 corresponding to the Hydrogen àtoms, 
and M corresponding to the methyl group. Akhough this is a seven-qubit register, 
the simulation of the Fano-Anderson model requires only three qubits. Consid- 
ering different practica! issues, such as the nuclei-nuclei interaction, the selection 
has been made as follows: The spin nucleus Ci represents qubit-1, the spin nu- 
cleus M represents qubit-2, and the spin nucleus C2 represents the ancilla qubit 
a. 

The idea is then to apply the desired elementary gates by sending appropriate rf 
pulses. Nevertheless, designing a pulse sequence to implement exactly the desired 
unitary transformation would require very long refocusing schemes to cancel out 
all the unwanted naturally occurring J couplings (free evolutions). Then, the 
overall duration of the pulse sequence increases and decoherence effects could 
destroy the signal. Therefore, a pulse sequence compiler is used to approximate 
the evolution and numerically optimize the delays between pulses to minimize the 
error introduced by such approximation. 

The approximate evolution is then applied to the corresponding initial state 

C P C Ml 
11^ + (T^2)1^0 , which is the density operator corresponding to 

the pure state |+)a 11102)- Since the identity part ]1 ^ is not an observable, the 

pseudo-pure State p^t = cr^-" l^'O^, with 1 = |1)(1| andO = |0)(0|, canbe used 
equivalently. Then p^^^, which is a deviation of the completely mixed state due 
to the high temperature of the ensemble (Sec. 12.6.11) . can easily be prepared using 
already developed NMR techniques IILKC02I . 

As mentioned in Sec. 12.2.21 the desired result of the QA is encoded in the po- 
larization (2(t^^) of the nucleus C2 (i.e., ancilla qubit). This component precesses 
at the C2-Larmor frequency uq^. To measure it, a Fourier transformation on the 
measured/ree induction decay FID (i.e., the decay in the polarization due to the 
contribution of different Larmor freqüències) must be performed, integrating only 
the peak located at uq_^ . Nevertheless, the absolute value of this signal is irrelevant 
since it depends on many experimental parameters such as the solution concen- 
tration, the probe sensitivity, and the gain of the amplifier. The relevant quantity 
is its intensity relative to a reference signal given by the observation of the ini- 
tial state Pinit- To get a good signal-to-noise ratio, each experiment (or scan) was 
done several times and the corresponding experimental data were added. More- 
over, to average over small magnètic fluctuations occurring within the duration 
of the whole experiment the scans of the reference experiment (i.e., the measure- 
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ment of the reference signal) are interlaced with scans of the actual complete pulse 
sequence. To increase the spatial homogeneity of the field over the sample, sev- 
eral automated shimming periods, consisting of fine tuning small additional coils 
located around the sample, have been inserted. 

In Fig. 12.201 I show the experimental results obtained for the evaluation of 
G{t) for Eko = -2, e = -8, V = and Eko = -2,e = 0,V = A, and for 
different vàlues of t. The duration of the optimized pulse sequences from the 
beginning of the initialization step to the beginning of the data acquisition, was 
97 ms. For comparison, the analytical form of G{t), as well as the simulated 
data points (i.e., data points obtained by simulating the quantum algorithm on a 
conventional computer) are also shown in the figure. 

In Fig. 12.211 I present the experimental results obtained for the evaluation 
of S{t), to obtain the corresponding eigenvalues for the Hamiltonian, with = 
—2, e = —8, and V = 0.5. The pulse sequence applied is the one corresponding to 
the quantum circuit shown in Fig. 12.181 with the corresponding refocusing pulses. 
In Fig. 12.211 I also show the analytical and simulated data points. The DFT of 
the experimental data is shown in Fig. 12.221 revealing the expected peaks at the 
frequency corresponding to the two one-particle eigenvalues of Eq. 12.781 for the 
above parameters. 

The close agreement between the experimental results and the corresponding 
simulations, using the refocusing pulses, suggests that the main contribution to 
errors comes from the incomplete refocusing scheme used in the optimization 
procedure. Therefore, increasing the number of refocusing pulses might have led 
to more accurate results even if they would have increased the overall duration of 
the pulse sequence. 

2.8 Summary 

Throughout this chapter, I have addressed several broad issues associated with the 
simulation of physical phenomena with QCs. In particular, I presented efficient 
ways to map the algebras of operators associated to the physical system to be 
simulated onto the àlgebra of Pauli operators. These mappings are sufficient to 
establish the equivalence of the different physical models to a universal model of 
quantum computation: The conventional model. 

I also explored various issues associated with simulations (Sec. 12.31) . remark- 
ing that every step in a QA must be performed efficiently to have an efficient 
simulation. Although a QC can simulate some quantum physical systems more 



2.8 Summary 



55 




1 2 3 1 2 




Time (s) 



Figure 2.20: Real and imaginary parts of the correlation function G{t) of Eq. lZSÜl 
The top panels show the results when the parameters in Eq. l2.78l are e^o = ~2, e = 
—8, V = A. The corresponding parameters Ai, A2, 6* can be determined using Eqs. 
12.841 and 12.851 The bottom panels show the results for = —2, e = 0, V = 4. 
The (black) sòlid line is the analytic solution, the red circles are obtained by the 
numerical simulation (including the refocusing pulses), and the blue circles with 
the error bars are experimental data. 
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Figure 2.21: Real and imaginary parts of S{t), for Eko = —2, e = —8, and 
V = 0.5 in Eq. 12.781 The (black) sòlid line corresponds to the analytic solution. 
The red circles correspond to the numerical simulation (using refocusing pulses) 
and the blue circles with the error bars are experimental data. S{t) has been 
measured using the network of Fig. 12. 181 with a = (e + eko)/2. 
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Figure 2.22: Discrete Fourier transform of the real part of the experimental data 
of Fig. 12.211 The position of the two peaks corresponds to the two eigenvalues of 
the Hamiltonian of Eq. 12.781 for Sko = —2, e = —8, and V = 0.5. Numbers in 
parentheses denote the exact solution. The size of the dots representing experi- 
mental points is the error bar (see Appendix|Bl). An upper bound to the error in 
the frequency domain is ^ 0.5, which was determined by the resolution of the 
spectrum due to the time sampling of the simulation. 
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efficiently than its classical analogue, I showed that many challenges still remain 
to prové this statement in a general way (Sec. l2.5l) . As an example, I presented the 
QS of the two-dimensional fermionic Hubbard model, where there is no known 
way to obtain its ground state energy efficiently using a QC. 

Finally, I described the experimental implementation of a QA using a liquid- 
state NMR QIP (Sec. l2.7l ). This experiment allows one to understand the advan- 
tages and disadvantages of simulating physical systems with today's QCs, and, in 
particular, to understand the power of quantum computation. 



Chapter 3 



Quantum Entanglement as an 
Observer-Dependent Concept 

"...Maximal knowledge ofa total (quantum) system does not neces- 
sarily include total knowledge ofall its parts, not even when these are 
fully separated from each other and at the moment are not influencing 
each other at all." 

E. Schródinger ( 1935). 

Quantum Entanglement (QE) is referred to the existence of certain correla- 
tions in a quantum system that have no classical interpretation. This concept was 
first introduced by E. Schródinger IISch35l as the essence of quantum mechanics, 
and is responsible for many counterintuitive physical processes like the violation 
of the local realism. Naturally, it has been the main focus of philosophical dis- 
cussions since the early days of quantum mechanics, and it is now known that 
entanglement is the defining resource that allows one to perform certain proto- 
cols like quantum cryptography, quantum teleportation, and even more efficient 
computation. For this reason, QE has been one of the most important subjects of 
study in QIT during recent years. Nevertheless, its analysis requires a good under- 
standing of the conceptual foundations of the quantum theory. For this purpose, a 
historical introduction to the subject is given in the foUowing sections. 

The EPR Paradox 

In 1935, Einstein, Podolsky and Rosen designed an experiment, the so called 
EPR paradox IIEPR35I . to prové that quantum mechanics was an incomplete de- 
scription of physical reality. Neither they nor others agreed that the (probabilistic) 
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outcome of a measurement performed on a quantum system was not uniquely de- 
termined by its quantum state |^). They believed, instead, that the result of a 
measurement was a property associated with the quantum system right before the 
measurement was performed. 

To prové this statement, they considered a system composed of two quantum 
partides where, ahhough the position and momentum of each particle were not 
well defined (uncertainty), the sum of their positions and the difference of their 
momenta were. For pedagogical purposes, I consider here a simpler version of 
this experiment by means of the conventional model IIBohSlI where the qualitative 
results are equivalent to those of the EPR paradox. This simplified model consists 
of two qubits initially prepared in some state (Fig. l3.1l) 

IBell) = -^[|0^1e)-|U0B)]. (3.1) 

[Such a state represents, for example, the singlet state (total spin 0) of two-spin 
1/2.] Assuming that a measurement (in the logical basis) is performed on qubit A, 
quantum mechanics says that the outcome of such measurement projects the state 
|Bell) onto the state |0^1b) with probability 1/2, and onto the state ll^Og) with 
the same probability (Sec. 12.1.11) . Thus, a measurement on qubit A projects also 
the state of qubit B, such that the outcome of a later measurement performed on it 
can be predicted with certainty. 




Figure 3.1: Bohm's representation of J. Bell's experiment. Two qubits in a pure 
quantum state are separated and their polarizations are locally measured along 
three different directions. The measurement outcomes violate the Bell inequality 
of Eq. lXÏÏl Here M denotes the polarizers. 

What worried Einstein and others is the fact that, in principle, both partides 
(qubits) could be very far apart in physical space. Therefore, a person located 
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close to qubit A could gain information about the state of qubit B almost im- 
mediately, violating the principle of locality; that is, some information could be 
propagated at infinitely large speed. If such were the case, the existence of these 
quantum correlations would be against the relativity theory, and the understand- 
ing of the physical world would be completely different. Nevertheless, a deeper 
analysis shows that such violation does not exist and that quantum mechanics is 
right. Below, I present the experiments suggested by Bell to show that quantum 
correlations exist in nature and cannot be explained by means of classical theory. 

Bell's Inequalities 

The idea of existence of hidden variables which determine the outcome of a 
measurement in a quantum system was proved to be incompatible with quantum 
mechanics by J. Bell IIBel64l in 1964. Basically, he proposed a similar experiment 
to the one described above but such that the measurements could be performed in 
any direction, that is, in a possible basis different to the logical one. For example, 
when working with photons, such measurements could be performed by using 
polarizers rotated by different angles. 

Assuming that qubit A is measured in some basis denoted by the vector l^i 
and qubit B is measured in the basis denoted by 1^2, such measurements project 
the State of either qubit in the corresponding direction with eigenvalues ±1 (Sec. l2.1.n) . 
In other words, the state after the measurement of the jth qubit can be represented 
as a vector pointing in the direction l^j (+1) or in the opposite direction —l^j 
(-1). 

The idea of Bell was to obtain the function A{l^i, 1^2), which denotes the 
average product of the outcomes of the corresponding measurements for different 
directions T^j. It can be proved that if T^i = 1^2, then = — Ifora 

Bell state of the form of Eq. 13.11 This is because such state can be written in the 
same form for other locally rotated basis. For example 

^ [\IaOb) - \OaIb)] = ^ [\+a-b) - hA+B)] , (3.2) 

where \±j) = -^[|0j) ± are the eigenstates of the Pauli flip operator cr^ with 
eigenvalue ±1. Then, a local measurement of the Bell state in the x direction 
shows the same results as a measurement in the z direction (logical basis). 

Equivalently, for the Bell state |Bell) one obtains A{l^i, —l^i) = +1, and 
for arbitrary directions, A{l^i, 1^2) = —^1-^2; that is, the usual scalar product 
between vectors in real space. Interestingly, such a result is associated with the ex- 
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istence of correlations of quantum nature and cannot be obtained by any classical 
theory, including hidden variables. 

To show this, J. Bell first assumed that the complete state of both qubits was 
characterized by the existence of uncontroUable hidden variables, denoted by A. 
He also assumed that the measurement outcome of qubit B was independent of 
the oríentation l^i where qubit A is measured (i.e., locality). Therefore, there 
exist two functions My^{l^i, A) and Mq{1^2, which correspond to both mea- 
surements outcomes, with vàlues 

MaÍ^u A) = ±1 ; Msir^, A) = ±1. (3.3) 

In particular, both functions satisfy the desired result if "r*i = 1^2 (see Eq. l3.ll) : 

MA^u A) = -MbÍ^i, A) va. (3.4) 

Nevertheless, if p{\) is the probability distribution for the hidden variable, with 
/ p{\)dX = 1, the average of the product of both outcomes should be 

A{r,, ^2) = j p(A)M^(ri, A)MB(r2, X)dX, (3.5) 

and fromEq.EH 

^(^1, ^2) = - y" p{X)MA^^, A)M^(r2, X)dX. (3.6) 

Considering another unit vector and considering that the three integers M_4 = 
±1 satisfy 

M^(ri)[M^(r2) - M^(r'2)] = ±[1 - M^(r2)M^(r;)], (3.7) 

and after some simple algebraic manipulations of Eq. l3?6l together with [M_a{1^i, A)]^ = 
1, the foUowing Bell inequality is obtained: 

|A(^i,^2) - A(^i,^'2)| < 1 + ^(^2,^2)- (3.8) 

Remarkably, this Bell inequality sometimes is violated in a quantum mechan- 
ical system due to the existence of non-classical correlations. For example, if the 
vectors 'r'i and 1^2 are orthogonal, and "r^2 between them making a 45° angle 
("Fig. l3.ll) . then for the quantum state of Eq. 13.11 

^(r'i,r2) = 0, ^(^1,^3) = ^(^2,^3) = -cos(7r/4), (3.9) 
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and 



cos(7r/4) ^ 1 



cos(7r/4). 



(3.10) 



Interpretation 

The violation of Bell's inequalities was experimentally proved in the 1980s 
using photon atòmic transitions IIAGR82I . The results obtained were in excellent 
agreement with the predictions of quantum mechanics, assuring that no hidden 
variable theory could overcome the problem of nonlocality. However, an analysis 
of the measurement process telis one that the existence of non-classical correla- 
tions does not imply a violation of causality, as Einstein, et. al., stated. In other 
words, the almost instant influence of the measurement outcome of one qubit in 
the measurement outcome of the other, cannot be regarded as the existence of 
some sort of communication (i.e., travelling Information) between both qubits: 
The outcome of the measurement (projection) is completely random. 

Many different analysis about the nonclassical properties of quantum systems 
are discussed in almost any book about quantum mechanics. In most cases, the 
problem of measurement, which is not analyzed in detall in the present work, is 
also considered. In this chapter, I mainly focus on the study of these nonclassical 
correlations which define QE. In the first place, I review the usual (or traditional) 
notion to prepare for a more general concept denoted as generalized entangle- 
ment, which is also presented here. Finally, I present some examples to show how 
generalized entanglement can be used in a more general framework. 

3.1 Quantum Entanglement 

The Standard setting for studying entanglement usually involves a quantum system 
S composed of two distinguishable subsystems A (for Alice) and B (for Bob), 
with Hilbert spaces denoted by Ha and Hq, respectively. Then the total Hilbert 
space for the joint system is Hs = ® Hb, with (g) the usual tensor product. A 
quantum state p (i.e., its density operator) is said to be separable if and only if 



where J2sPs = 1' '^i^h Ps (> 0) being the corresponding probabilities. The 
density operators pf = \^%){^%\ and pf = \í>^)Í^b\' which describe pure states 




(3.11) 



s 
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of Alice and Bob, respectively, need not describe orthogonal states for different 
integers s. A quantum state of the joint system is then separable whenever it 
can be written as a probabilistic (convex) combination or a mixture of product 
(separable) states. For example, if both A and B are single qubits, the pure state 

1^) = l [Mb) + \OaIb) + \IaOb) + lUlfí)] (3.12) 
is separable because it can be rewritten as 

W = M^«M£^. (3,3) 
V 2 V 2 

Interestingly, separable states can be prepared by means of local operations 
(e.g., rotations, local measurements, etc.) and classical communication between 
the different parties of the system. Then, they do not possess correlations of quan- 
tum nature and cannot be distinguished from classical states. This property leads 
to the foUowing definition: A quantum state p of a composite system is said to 
be entangled if it is not separable and unentangled otherwise. For example, the 
Bell state of Eq. 13. ll is entangled because there is no basis where it can be writ- 
ten as a product state. Naturally, this definition of entanglement applies not only 
to bipartite systems but to all systems composed of many (i.e., more than two) 
distinguishable subsystems. 

Finding a decomposition like Eq. 13.1 ll for a given quantum state p of a joint 
system, or even showing that p describes an entangled state (i.e., no such decom- 
position exists) is, in general, a very difficult task. In fact, algebraic methods to 
check if p is entangled or not only exist for the case of bipartite systems, and for 
Alice and Bob being single qubits or two level subsystems. However, if the state 
is pure (i.e., p = |V^)(?/'|), simple methods to check separability exist. 

3.1.1 Separability and von Neuman Entropy 

The Schmidt decomposition of a pure quantum state \íJj) provides a useful tool to 
check its separability. When using this decomposition, any state {i/j) of a bipartite 
quantum system can be written as IIPer98l 



i=i 



(3.14) 
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where the pure states and are orthonomial states of subsystems A and 
B, respectively; that is 

{<P'Atò = = ^n' VJ,/ e\l---R]. (3.15) 

Without loss of generality, the coefficients Cj 7^ of Eq. 13. 141 are real (i.e., the 
phases have been absorbed in the corresponding states) and, from the normaliza- 
tion of 1^), they satisfy Ylif=ii.^jY — 1- '^a denote the dimensions of the 

Hilbert spaces Ha and Hb, respectively, the integer R satisfies 

R<mm{da,db}. (3.16) 

Obviously, the pure state \íp) is entangled whenever i? > 1 in the decompo- 
sition and unentangled otherwise {R = 1). Moreover, [tp) is entangled whenever 
it looks mixed from the point of view of either observer (Alice or Bob). In other 
words, when obtaining the reduced density operator of one subsystem by tracing 
out the other, such density operator describes a mixed state when \íp) is non- 
separable or entangled. This is only vàlid for pure states of the joint system. For 
example, the reduced density operator p-^ associated to Alice is 

R FI 

= Tr^(p) = Y.i'í^'Mim) = T.^'^^rww- (s-iv) 

Then, if (cj)^ < 1 Vj (or equivalently R > 1), the reduced density operator p-^ 
describes a mixed state and {ip) is entangled. Similar results are obtained from the 
point of view of Bob. 

A useful way to quantify the entanglement of a pure bipartite state \ip) is given 
by the von Neuman entropy of either p-^ or p'^. Its definition is 

R 

EsW = -Tr {p^. \og, p-^) = -Tr (p^. log^ p^) = - ^(c,)^ log^lc,)^. 

i=i 

(3.18) 

Es is zero (minimal) for a product state and takes its maximum value (Es = 1) 
for maximally entangled states, such as the Bell state of Eq. I3TT1 Moreover, Es 
remains invariant or decreases under any operation on performed locally by 
Alice or Bob. As I will show, this is an important property for a measure of 
entanglement. 

All the concepts described in this section can be naturally extended for mul- 
tipartite systems. Then a pure quantum state IíJj) is entangled whenever it looks 
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mixed to, at least, one of the observers associated to one of the parties. A von 
Neuman measure can be defined accordingly, by measuring the entanglement of 
every single party with the rest of the system. 



3.1.2 Mixed-State Entanglement and the Concurrence 

As mentioned, checking the separability of a mixed quantum state is a difficult 
task because of the many equivalent ways that the density operator can be decom- 
posed. Nevertheless, certain (classical) algorithms to calculate the entanglement 
of a mixed state exist, but their complexity scales with the dimension ds of the 
Hilbert space H associated to the joint quantum system. Such dimension is known 
to scale exponentially with the system size. Thus, these algorithms can only be 
applied to study the entanglement of small systems. 

For a mixed state of abipartite system p = J2sPs\^s) {i's], that is, the system 
being in the pure state {ips) with probability Ps < l, the entanglement of formation 
E{p) is defined as 



E{p) = min 



(3.19) 



where Es{^s) is the von Neuman entropy defined in Eq. 13.181 Each possible 
decomposition of p corresponds to a certain amount of entanglement, so the min- 
imum needs to be obtained. (Equation 13. 191 is trivially extended for multipartite 
systems.) For example, the mixed state of a two-qubit system 

p = ^|Bell)i(Bell|i + ^|Bell)2(Bell|2, (3.20) 

where |Bell)i = |0^0b) + and |Bell)2 = |0^0b) — are maximally 

entangled states (i.e.. Bell states), seems to be maximally entangled. However, p 
is actually separable and thus unentangled : 

P = 1\0a0b){0a0b\ + ^\1aIb){IaIbI (3.21) 

that is, a mixture of product states. 

As mentioned, the entanglement for a mixed state of two qubits {A and B) can 
be exactly computed without calculating the minimal value of Eq. l3.19IIIWoo98l . 
For this purpose, the time reversal operation needs first to be applied. That is, if p 
is the state of the two-qubit system, then 

P=KX)P*KX), (3.22) 
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where p* is the complex conjugate of p, is the corresponding spin-flipped state. 
Remarkably, the entanglement of formation (Eq. 13.191 1 for this system can be ex- 
prés sed as 

E{p)=S{C{p)), (3.23) 

where the real function S{C) is monotonically increasing with C, and C{p) is 
denoted as the concurrence of p IIWoo98l . (Due to this property, one can use 
the concurrence as a measure of entanglement in the two-qubit system, instead.) 
Defining the operator 

R = \IV~PPVp^ (3-24) 

the concurrence is defined by 

C{p) = max{0, Al - A2 - A3 - A4}, (3.25) 

where the A^s are the eigenvalues of R in decreasing order (Aj > 0). Then, C{p) 
takes its maximum value [C(p) = 1] for Bell states (Eq. l3.1l) and vanishes for any 
separable (mixed) state. 

3.1.3 Measures of Quantum Entanglement 

A good measure of entanglement E(p) for a multipartite quantum system, such 
as the entanglement of formation of Eq. 13.191 needs to satisfy certain require- 
ments flVidOOil . First, such a measure must take its minimum value, E(p) = 0, 
whenever p describes a separable state. Second, unitary local operations, local 
measurements, and classical communication between different parties in the sys- 
tem (usually referred to LOCC operations) cannot increase E(p). It is reason- 
able that LOCC operations do not transform, for example, a separable state into 
a nonseparable one. Moreover, local measurements project the state lowering its 
entanglement. 

Other requirements also have to be considered when defining a good measure 
of entanglement. These include continuity, convexity, additivity, and subadditiv- 
ity. For example, if two density operators describe almost the same state, they 
must have similar amounts of entanglement. That is, E(p) must be a continu- 
ous function of p. Also, the entanglement of a linear combination of two density 
operators pi and p2, that is, p = ppi + (1 — p)p2 (with < p < 1) must satisfy 



E(p)<pE(pi) + (l-p)E(p2). 



(3.26) 
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A function satisfying Eq. 13 .261 18 said to be convex. Basically, the reason for such 
convexity is because an operator like p tends to have less entanglement than the 
corresponding operators pi and p2· 

It is not very clear how important the requirements of additivity and subad- 
ditivity are. The additivity property states that if the system is composed of m 
identical subsystems, all being in the state given by the density operator r, then 



The subadditivity property states, however, that if the density operator p of the 
total (multipartite) system can be expressed as p = ri ® T2, then 



In fact, the entanglement of formation for multipartite systems, as defined by 
Eq l3.191 does not satisfy the additivity property. 

3.2 Generalized Entanglement 

The main purpose of generalized entanglement (GE) is to extend the concepts of 
traditional entanglement to a more general setting by defining new measures that 
can be applied to any quantum system, even when a nontrivial subsystem decom- 
position exists. In this way, GE is considered as an observer-dependent property 
of a quantum state and, as in the usual case, it is determined by a preferred set of 
observables of the quantum system under study. It is expected that this new con- 
cept allows for a better understanding of non-classical correlations in quantum 
mechanics. 

In this section, I rigorously define GE HBKOOBl IBKO04II in terms of reduced 
states and I analyze the important case where the preferred set of observables 
belong to a certain Lie àlgebra (i.e., an àlgebra closed under commutation). Some 
examples to show how GE works in different quantum systems will be presented 
and, in particular, I will show how the traditional notion can be recovered. 

3.2.1 Generalized Entanglement: Definition 

Assume that f) is a finite set of observables; that is, \] = {Òi, O2, ■ ■ ■ , Òm}, with 
Òj = Ój the Hermitian operators that linearly map quantum states of a Hilbert 
space H into quantum states of the same space. Such set of observables is usually 



E(p) = H®T=ir) = mE(r). 



(3.27) 



E(p) < E(ri) + E(r2). 



(3.28) 



3.2 Generalized Entanglement 



69 



intrinsically associated to the quantum system S under study and depends on its 
nature, control access, superselection rules (e.g., particle number conservation), 
etc. For a quantum state with corresponding density matrix p = J^sPsli's) {i's] 
(Ps > 0; J2s Ps = 1)' its f)-state (reduced state) is defined to be a linear functional 
A on the operators of f) according to 

A(Ò,) = MpÒj) = (Ò,) Vj G [1 ■ ■ ■ M], (3.29) 

with Oj E [), and (Oj) its expectation value over the state p. In particular, if f) 
denotes the set of all linearly independent observables associated to the quantum 
system, the set of f)-states completely determine the state of the system, since p 
can be exactly recovered from those expectations. 

Considering the set of all f)-states, it can be shown that such set is closed 
under convex (or probabilistic) combination. In other words, if are f)-states 
associated to different density matrices pk, then Pk^k is also a possible [)-state 
for any probabilistic distribution {p^ > 0, J2kPk = ^ compact 

then all f)-states can be obtained as combinations of extremal states. Extremal 
States are those defined as the reduced states that cannot be written as a convex 
combination or mixture of two or more í)-states. In general, extremal states are 
least uncertainty states so they are also referred to í)-pure states. 

With these definitions and properties, a pure state (i.e., p = \íjj){tp\)ís said 
to be generalized unentangled relative to the distinguished set of observables f) if 
its reduced state, that is, the state defined by the expectations of the operators in 
f), is pure or extremal. Otherwise, [tp) is said to be generalized entangled relative 
to f). 

Although such definition can be applied to any set of observables, relevant 
physical applications are usually found when the set f) is a Lie àlgebra of operators. 
In the foUowing, I describe this important case while giving a short description of 
Lie àlgebra theory. 



3.2.2 Generalized Entanglement and Lie Algebras 



I now focus on the case when i) = {Òi, Ó2, ■ ■ ■ , Òm} is a finite real M-dimensional 
Lie àlgebra of observables acting irreducibly on H (the Hilbert space), with bracket 
given by 

M 



(3.30) 



k=l 
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(Here, no distinction between the abstract Lie àlgebra isomorphic to f), and the 
concrete matrix Lie àlgebra f) of observables acting on the Hilbert space Tí is 
made.) The corresponding induced Lie group of f) is given by the map X — > e*"^, 
with X G f). Also, [) is usually assumed to be semi-simple, that is, with no 
commutative ideals ^ 

The projection map of a quantum state, with density matrix p, onto f) can be 
uniquely defined by the trace inner product as 

M 

P^V^{p) = J2iÒj)Ò„ (3.31) 

i=i 

where {Oj) = Tr{pOj) is the corresponding expectation value. Notice that only 
the [) -state associated to p must be known to build Vi,. Similarly, the relative purity 
or [)-purity is defined as the squared length of the projection; that is 

A(P) = Tr (Pj(p)) , (3.32) 
and if the operators are Schmidt orthogonal satisfying Tr (^ÒjÒk^ = 5jk, then 

N 

P^{p) = Kj2{Ò,)'. (3.33) 

Here K is a real constant for normalization requirements. The P)-purity as defined 
by Eq. l3.33l is a group invariant function: 

P^(p) =P^(e*V"^''); Xei). (3.34) 

Interesting consequences from the definition of the relative purity are obtained 
when the matrix Lie àlgebra f) acts irreducibly on the Hilbert space H associated 
to the system, that is, the representation of i) is irreducible. In such a case, a 
pure quantum state \íJj) is extremal or [)-pure if and only if the [;)-purity takes 
its maximum value IIBK0Ü3II . Such states are the so-called generalized coherent 
States (GCSs) of f) and are constructed as 



|GCS) = exp 



■ M 



|ref), (3.35) 



' An ideal of f) is a subalgebra 3 C f) invariant under the commutation with any element of f). 
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where Q E M, and |ref) is a reference state corresponding to the highest (or 
lowest) weight state of (). GCSs are a generalization of the traditional coherent 
States of the harmònic oscillator or the radiation field IISZ02I . where in those cases 
I ref ) is the vacuum or no-excitation state. (The first term on the right hand side of 
Eq. l3.35l is a general displacement operator.) 

In general, the reference state | ref) is a well defined object in finite semi-simple 
Lie algebras, after a particular Cartan- Weyl (CW) decomposition is performed. In 
a CW basis, the àlgebra is written as 



(3.36) 



Here, Í)d = {hi, /i2, ■ ■ ■ , K}, with hj = h)-, is a Cartan subalgebra (CSA) of f) 
defined as the biggest set of commuting operators in f), and the integer r defines the 



rankoff). The sets f)+ = {Ea^.Ea^ 



with En 



t 



Ecu) and f)_ = {E_ai,E_a2, ■ ■ ■ , E_a^}, 



E-ai , are usually referred to the sets of raising and lowering op- 



erators, respectively. By notation, r 
sponding commutation relations are 



21 



M, the dimension of f). The corre- 



hk, hk' 



E E 



En,, En 



0, 

r 
k=l 



(3.37) 
(3.38) 

(3.39) 
(3.40) 



where the vectors a,- 



■ , G W are defined as the roots of the 



àlgebra, and Njj' depends on the corresponding roots. 

Equation 13.371 states that the operators in f)^) can be simultaneously diago- 
nalized. Their eigenstates in a given representation are the corresponding weight 
States. A weight state \(f)P) of i) satisfies the eigenvalue equation 



hkm = eim VA;G[l···r], 



(3.41) 



where the eigenvalues e]^ G M are defined as the weights. Moreover, the vectors 



are defined as the weight vectors. 



Due to the commutation relations of Eq. 13.381 a raising operator acting over a 
weight state either annihilates it (i.e., Ea 10'') = 0) or maps it to a different weight 
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State: 

hkÈ^^m = {È^^hu + = [a] + el)È^^m. (3.42) 

Similar results are obtained for the lowering operators È^aj from Eq. 13.391 The 
simplest case is given by the àlgebra su(2) = {a^, 0"+, o'_}, where the spin-1/2 
representation is given by the Pauli matrices defined in Eq. l2.1l 

A weight vector "e^*' = (ef , 63, ■ ■ ■ , e^) is said to be positive if the first non- 
zero component is positive. Then, a weight vector "e*^ is said to be larger than a 
weight vector "e^^ if they differ on a positive weight vector. In particular, the root 
vectors aj = (a], a^, ■ ■ ■ , ap in Eq. 13.361 are weight vectors of a representation 
of i) called adjoint representation (AppendixO. With no loss of generality, these 
vectors have been chosen to be positive (or equivalently, —Uj to be negative) in 
the CW decomposition. 

The definition of positivity naturally extends to any weight vector in any rep- 
resentation of í). Therefore, a highest weight state \ HW) of f) is defined as the state 
with highest weight vector in the representation. Similarly, a lowest weight 
state I LW) is defined as the state with lowest weight vector in the representation. 
These states are annihilated by every raising and lowering operator, respectively: 

E„JHW) =Í_,JLW) = 0. (3.43) 

Both, |HW) and |LW), are possible reference states |ref) to generate the GCSs or 
generalized unentangled states of f) (Eq. 13.35b . With no loss of generality, I will 
mainly use |HW) as the reference state. Since its weight vector is positive, this 
state is then the unique ground state of a Hamiltonian 

r 

H = Y,lkhk; ikeR, (3.44) 

k=l 

where the sign of the coefficients 7^ depends on the weight. Therefore, every GCS 
is the unique ground state of some Hamiltonian in f) IIBKO03I . 
A common uncertainty measure for the àlgebra f) is given by 

M 

(AF)2 = ^(Ò|) - (Ò,)^ (3.45) 

Remarkably, GCSs of f) are also minimal uncertainty states. When f) acts irre- 
ducibly on H, the first term on the rhs of Eq. 13.451 is a Casimir, i.e. invariant. 
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while the second term is the í)-purity. Because GCSs have maximum purity, the 
corresponding uncertainty is minimal. 

It is important to realize that the relationships just mentioned between max- 
imal purity, generalized coherence, and generalized unentanglement established 
for a pure state relative to an irreducibly represented Lie àlgebra i) do not automat- 
ically extend to the case where f) acts reducibly on H. In fact, for a semi-simple 
Lie àlgebra f), a genèric finite dimensional representation of i) may be decomposed 
as a direct sum of irreducibly invariant subspaces, H = (Bel·le, with each of the 
Hi being in turn the direct sum of its weight spaces. Then, the àlgebra acts irre- 
ducibly on each He, with the corresponding irreducible representation (irrep). In 
particular, every irrep appearing in the decomposition has a highest (and lowest) 
weight state, and for each of these irreps there is a manifold of GCSs constructed 
as the orbit of a highest weight state for that irrep. Therefore, these GCSs will not 
satisfy in general the extremality property that defines generalized unentangled 
States. Indeed, the extremal weight vectors, which correspond to generalized un- 
entangled States, need not all have the same length (í)-purity). Also, the minimal 
uncertainty property could be lost. Maximal purity remains then as a sufficient, 
though no longer a necessary, condition for generalized unentanglement. Nev- 
ertheless, all the statements obtained for the irreducible case still apply for the 
examples and problems presented in this thesis. 



3.2.3 Generalized Entanglement and Mixed States 

For mixed states on H, the direct generalization of the squared length of the pro- 
jection onto f) as in Eq. (13.321) does not give a GE measure with well-defined 
monotonicity properties under appropriate generalizations of the LOCC set trans- 
formations as explained in Sec. l3.1.3l[|BKO03l . A proper extension of the quadratic 
purity measure defined by Eq. 13.321 for pure states to mixed states may be natu- 
rally obtained via a Standard convex roof construction. If p = J2Ps\^s) {''Ps], 

s 

with J2Ps = ^ ^iid J2p's < 1' the latter is obtained by calculating the maxi- 

s s 

mum fi-purity (minimum entanglement) over all possible convex decompositions 
{ps, \ips)} of the density operator p as a pure-state ensemble. In general, similarly 
to what happens for most mixed-state entanglement measures, the required ex- 
tremization makes the resulting quantity very hard to compute. Nevertheless, the 
cases studied in this thesis are always related with pure states and no extremization 
process is required. 
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3.3 Relative Purity as a Measure of Entanglement 
in Quantum Systems 

In order to understand its meaning as a measure of entanglement for pure quantum 
States, I now apply the definition of relative purity to the study of non-classical 
correlations in different physical systems. First, I will focus on spin systems, 
showing that for particular subsets of observables, the f)-purity can be reduced to 
the usual notion of entanglement: The pure quantum states that can be written 
as a product of states of each party are generalized unentangled in this case. For 
this purpose, and because they will also be needed in future cases, I introduce the 
foUowing representative quantum states for spins of magnitude S: 



where is'l, 5*2, ■■■ , Sn) = \Si)i ® 182)2 ® ■ ■ ■ ® \Sn)n is a product state, and 
\Si)i denotes the state of the ith party with z-component of the spin equal to Si 
(defining the relevant computational basis for the ith subsystem). 

However, for other physically natural choices of observable sets a different no- 
tion of QE is obtained. In this way, one can go beyond the usual (distinguishable) 
subsystem partition. For this reason, I will focus later on the study of the í)-purity 
as a measure of entanglement for fermionic systems. This is a good starting point 
to understand how such a measure can be applied to quantum systems obeying 
different particle statistics and/or described by different operator languages. 

3.3.1 Two-Spin Systems 

The simplest system to study a measure of entanglement is composed of two 
spins. For simplicity, I begin by studying the GE of a two-qubit system (two- 
spin-1/2, A and B), where the most general pure quantum state can be written as 
= alO^Og) + 6|0_41b) + c|1_40b) with the complex numbers a, h, c, 

and d satisfying |ap + |6p + |cp+ = 1. For a spin-1/2 system the associations 
|0) = IT) = l|) |1) = li) = are commonly considered. The traditional 




S, S, ■ ■ ■ , S) , 



(3.46) 
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measures of pure state entanglement in this case are well understood (Sec. 13.11 ). 
indicating that the Bell States |GHZ|) = ^ [|0^0b) + |U1b)] IIEPR35I (andtheir 
local spin rotations such as the state of Eq. 13.11) are maximally entangled with re- 
spect to the local Hilbert space decomposition Hb- On the other hand, 
calculating the purity relative to the (irreducible) Lie àlgebra of all local observ- 
ables i) = su(2)^ © su(2)g = {cr^; i : A,B; a = x,y,z} classifies the pure 
two-qubit States in the same way as the traditional measures (Fig. 13.21) . Here, the 
operators a-^ = ^ and crf = 11 ® cTq are the Pauli operators acting on qubits 
A and B, respectively. In this case, Eq. 13.321 simplv gives 



where Bell's states are maximally entangled (P^ = 0) and product states of the 
form lip) = \(f)X)yx (8) \4>b)b (GCSs of the àlgebra [) = su(2)_4 © su(2)g) are gen- 
eralized unentangled, with maximum purity. Therefore, the normalization factor 
K = 1/2 may be obtained by setting P(, = 1 in such a product state. As ex- 
plained in Sec. 13.2.21 Pf, is invariant under group operations, i.e., local rotations 
in this case. Since all GCSs of f) belong to the same orbit generated by the ap- 
plication of group operations to a particular product state (a reference state like 
lO^Ofí) = |i, |)), they all consistently have maximum P)-purity (Pf, = 1). 

Another important insight may be gained by calculating the purity relative to 
the àlgebra of all observables, f) = su(4) = {(Tq,(J^ ® cr|; i = A,B; a, (3 = 
X, y, z} in the case of two-qubit system. One finds that any two-qubit pure state 
\^) (including the Bell's states of Eq. 13.11 ) is generalized unentangled (P[j = 1; 
See Fig. 13.21) . This property is a manifestation of the relative nature of GE, since 
considering the set of all observables as being physically accessible is equivalent 
to not making any preferred subsystem decomposition. Accordingly, in this case 
any pure quantum state becomes a GCS of su(4). 

In Fig. 13.21 1 also show the GE for systems of two parties of spin-S" relative to 
different algebras. It is observed that the purity reduces again to the traditional 
concept of entanglement for higher spin if it is calculated relative to the (irre- 
ducible) Lie àlgebra of all local observables t) = su(25' + 1)^ © su(2S' + l)^. 
For example, when interested in distinguishing product states from entangled 
states in a two-spin-1 system, the purity relative to the (irreducible) àlgebra f) = 
su(3)^ © su(3)b = {A^ © 11^, 11-^ © (1 < a < 8)}, needs to be calculated. 
Here, the 3x3 Hermitian and traceless matrices A^ are the well known Gell-Mann 




(3.47) 
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Figure 3.2: Purity relative to different possible algebras for a two-spin-S' system. 
The quantum states |GHZ|) and |F|) are defined in Eqs. l3.46l 

matrices IIGeo99ll : 




X, = j^{ ol;Afi = 4^ 



A7 = 75 I -z I ; Ag 
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which satisfy Tr[AQ,A/3] = 5a^p. In this basis, the spin-1 states are represented by 
the 3-dimensional vectors 






|1) ' = j 1-1) = I • (3.48) 

Then, the relative purity for a genèric pure state becomes 

a=l i 

where (A^) denotes the expectation value of A^ in the state In this way, 
product States like = are generalized unentangled (Pf, = 1) and 

States like | GHZ^) (and states connected through local spin unitary operations), are 
maximally entangled in this àlgebra (P(, = 0). 

Different results are obtained if the purity is calculated relative to a subalgebra 
of local observables. For example, the two-spin-1 product state |0, 0) = |0) ® |0), 
where both spins have zero projection along z, becomes generalized entangled 
relative to the (irreducible) àlgebra su(2)^ © su(2)e of local spin rotations, which 
is generated by {5^; i : A,B] a = x,y, z}, the spin-1 angular momentum oper- 
ators Sa for each spin being given by 

5. = ^ I 1 1 I , ^„ = ^ ( z -z I , = ( I (3.50) 






Notice that access to local angular momentum observables suffices to operationally 
characterize the system as describable in terms of two spin-1 partides (by imagin- 
ing, for instance, the performance of a Stern-Gerlach-type of experiment on each 
particle). Thus, even when a subsystem decomposition can be naturally identified 
from the beginning in this case, states which are manifestly separable (unentan- 
gled) in the Standard sense may exhibit GE (see also Appendix|Dl). On the other 
hand, this is physically quite natural in the example, since there are no SU(2) x 
SU(2) group operations (local rotations from exponentiating f) = su(2)^©su(2)0) 
that are able to transform the state |0, 0) into the unentangled product state |1, 1). 
In fact, |0, 0) is maximally entangled with respect to su(2)^ © su(2)b (i.e., Pt, = 
0). 

The examples described in this section together with other examples of states 
of bipartite quantum systems are shown in Fig. 13.21 It is clear that calculating 
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the purity relative to different algebras gives information about different types of 
quantum correlations present in the system. 



The f)-purity distinguishes pure product states from entangled ones if it is calcu- 

N 

lated relative to the (irreducible) àlgebra of local observables f) = su{2S + l)j 

(see AppendixiDll because of the group invariance of the relative purity (Eq. l3.34l) . 
which, in this case, constitutes all local rotations. Here j denotes every subsys- 
tem (spin) of the total system. For example, in the previous section I denoted the 
subsystem 1 to be ^ (or Alice) and the subsystem 2 to be i3 (or Bob). Therefore, 
product States are GCSs and generalized unentangled relative to the set t) of all 
local observables. 

For example, the usual concept of QE in an A^-qubit quantum state can be 

N 

recovered if the purity is calculated relative to the local àlgebra f) = su(2)j = 

{al^al^al, - ■ ■ ,a^}, where the Pauli operators (a = x,y,z) were 

introduced in Sec. 12. 1.11 The local purity is then 



where the normalization factor was obtained by setting P^, = 1 in any prod- 
uct state of the form = 102)2®" ■ ■^\4'n)n (i-e-, a GCS of this àlgebra). 

With this definition, states like |GHZ^), [(|01) - |10))/V^]®" (with obvious no- 

2 

tations), and the well known cluster states introduced in Ref. HBROll . are 
maximally entangled (P(, = 0). Also, Eq. l3.51l can be shown to be equivalent to 
the measure of QE introduced by Meyer and Wallach in Ref. IIMW02I . 

In Fig. 13.31 1 present some examples of the purity relative to the local àlgebra 

N 

[) = 3u(2)j for a A^-spin-S system. I also show the purity relative to the àlgebra 

i=i 

of all observables su([2S' + 1]^), where any pure quantum state is a GCS, thus 
generalized unentangled (Pf, = 1). 



3.3.2 A^-Spin Systems 




N 



(3.51) 
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Figure 3.3: Purity relative to different algebras for a A^-spin-S* system. The quan- 
tum States |GHZ^), |W^),and |F^) are defined in Eqs . 13461 

3.3.3 Fermionic Systems 

The case of fermionic systems is important because it shows how the concept 
of GE can be widely used. The system considered here consists of N (spinless) 
fermionic modes j, where each mode is described in terms of canonical creation 
and annihilation operators cj, Cj, respectively, satisfying the anti-commutation 
rules of Eqs. l234l For instance, different modes could be associated with different 
sites in a lattice, or to delocalized momentum modes related to the spatial modes 
through a Fourier transform (i.e., wave vectors). In general, for any N x N unitary 
matrix V = ||%||, any transformation t-^ maps the original modes 

into another possible set of fermionic modes (Bogolubov transformation IIBR86II '). 
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The commutation relations of quadratic fermionic operators can be obtained 
using Eqs. l2.34l finding that 

[c\cj. c\ci] = ójkcjq - óaclcj . (3.52) 

Thus, the set of bilinear fermionic operators {CjCj,; 1 <]■,]' < N} provides a re- 
alization of the Lie àlgebra u(A^)^ in the 2^-dimensional Fock space Típock of the 
system. The latter is constructed as the direct sum of subspaces 7í„ correspond- 
ing to a fixed fermion number ri = 0, . . . , A^, with dim(7Y„) = A^!/[?7,!(A^ ^ '^)']· 
Here, it is more convenient to express u(A^) as the linear span of a Hermitian, 
orthonormal operator basis, which can be chosen as 

{{c]c^, + c],Cj) with 1 < j < j' < 
i{c]c^, - c],Cj) with 1 < j < j' < A^ , (3.53) 
V2{c]cj-l/2) withl<j<A^ 

(the large left curly bracket means "is the span of"). The action of u(A^) on Típock 
is reducible, because any operator in u(A^) conserves the total number of fermions 
n = {J2f=i '^iS)' turns out that the irrep decomposition of u(A^) is identical to 
the direct sum into fixed-particle-number subspaces Hn, each irrep thus appearing 
with multiplicity one. 

Using Eq. 13.321 the [)-purity of a genèric pure many-fermion state relative to 
u(A^) is 



j<j'=l 3=1 

(3.54) 

For reasons that will become clear shortly, the normalization factor was chosen to 
be K = 2/A^. In this case, the fermionic product states (Slater determinants) of 
the form 

I0) = n^-l^^^)' (3.55) 

m 

with |vac) denoting the reference state with no fermions and m labeling a partic- 
ular set of modes, are the GCSs of the u(A^) àlgebra IIGil74[|Per85ll . 

Because a Slater determinant carries a well defined number of partides, each 
GCS belongs to an irrep space 7í„ for some n\ states with different n belonging to 

basis for the matrix Lie àlgebra u{N) is given by all real N x N matrices. 



N . N 
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different orbits under u(A^). A fixed GCS has maximum f)-purity when compared 
to any other state within the same irrep space. Remarkably, it also turns out that 
any GCS of f) = u(A^) gives rise to a reduced state which is extremal (thus gener- 
alized unentangled) regardless of n; the f)-purity assuming the same (maximum) 
value in each irrep. Using this property, the normalization factor K = 2/N was 
calculated by setting P[j = 1 in an arbitrary Slater determinant. Thus, the purity 
relative to the u{N) àlgebra is a good measure of entanglement in fermionic sys- 
tems, in the sense that P(, = 1 in any fermionic product state, and P(j < 1 for 
any other state, irrespective of whether the latter has a well defined number of 
fermions or not. 

Due to the invariance of P(, under group transformations (Eq. l3.34l) . the prop- 
erty of a state being generalized unentangled is independent of the specific set of 
modes that is chosen. That is, if |0) = Ylj 4 |vac) is generalized unentangled with 

respect to u(A^), so is the state = Ylm '^Jnl^^'^)' "^^^^ ^In = J2j 'Vmjc\· There- 
fore, the u(A^)-purity is a measure of entanglement that goes beyond a particular 
subsystem decomposition in this case and only distinguishes fermionic product 
States from those which are not. 

For example, if the system has only = 4 sites (modes), then a fermionic 
state like |0) = -^(c|c2 + C3C4) |vac) is maximally entangled relative to the àlge- 
bra u(4) (that is, Pu(4) = 0) because there is no basis where it can be written as 
a fermionic product state. However, the state |0) = -^(cjcg + 4c4)|vac) is un- 
entangled with respect to u(4) (i.e., Pu(4) = 1) because it can also be written in a 
certain basis as the fermionic product state |0) = c|c||vac), with c\ = ■:^{cl + c\). 

3.4 Summary 

In this Chapter I have introduced a generalization of entanglement which provides 
a unifying framework for defining entanglement in an arbitrary physical setting. 
For this purpose, I first presented the usual notion of entanglement which referres 
to a particular subsystem decomposition of the total system. The usual concept 
can be naturally extended to a more general case by means of generalized entan- 
glement. The latter is regarded as an observer-dependent property of quantum 
States, being definable relative to any physically relevant set of observables for the 
system under study. 

In particular, I implemented some steps for the purpose of associating the 
theory of entanglement with the theory of generalized coherent states in the Lie 
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algebraic setting. I have shown that whenever the preferred set of observables 
constitutes a Lie àlgebra acting irreducibly on the Hilbert space associated with 
the system under study, GCSs are unentangled relative to such a set, and possess 
minimal uncertainty. Such properties are also obtained in the usual framework 
with respect to local observables. That is, the usual notion of QE is recovered 
from GE when choosing the àlgebra of all local observables, corresponding to a 
particular subsystem decomposition. 

Finally, some useful examples were presented to realize the dynamics of this 
innovative approach. 



Chapter 4 



Generalized Entanglement as a 
Resource in Quantum Information 

Because of the interesting non-classical features of entanglement, physicists have 
been studying this quantum mechanical property almost since the early years of 
quantum mechanics IISch351 IEPR35II . Entanglement is considered nowadays to 
be a fundamental resource in quantum Information processing, where it can be 
exploited to perform certain tasks like quantum cryptography, quantum telepor- 
tation, quantum simulations, etc, which are difficult or impossible in a classical 
setting. Nonetheless, it is not yet fuUy understood how or when QE can be used to 
perform more efficient computation. In Chap. |21 for example, I have shown that 
certain properties related with the simulation of physical systems, can be obtained 
more efficiently using a QC than a CC. As I will explain later, this is vàlid only 
if entangled states are involved in the simulation. However, if no entanglement is 
created at any step of a deterministic QA (i.e., preparation of a pure initial state, 
evolution, and measurement), the quantum state carried over the whole process 
remains a product state (unentangled) and such simulation can be performed on a 
CC using the same amount of resources. 

On the other hand, it is also known that the creation of entanglement is not a 
sufficient condition to claim that there is no classical algorithm able to perform the 
simulation efficiently IGKl IVidÜ^ . In fact, certain QSs which involve entangled 
States can be imitated on a CC with the same efficiency. So, which problems can 
be solved more efficiently using a QC? The answer to this question is not yet 
known: The lack of a general theory of QE that can be applied to any system 
and any pure or mixed quantum state is the principal reason behind the difficulty 
of understanding and comparing quantum with classical complexity in a general 
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case. 

In Chap. |3l I introduced the notion of GE which consti tutes a main step to- 
wards the developing of a general theory of QE. In the foUowing sections I will 
show how this novel concept leads to a better understanding about when QE can 
be used as a resource for more efficient computation. In particular, I will show 
that the creation of generalized entangled states relative to every small dimen- 
sional Lie àlgebra is a necessary condition for a QA to be more efficient than the 
corresponding (known) classical one. This result goes beyond the idea of creating 
entangled states in the traditional way to gain efficiency. 

This chapter is organized as foUows: First, I present how QE can be exploited 
to perform interesting protocols in quantum Information processing. Second, I 
focus on the study of the relation between the concept of GE with classical and 
quantum complexity, showing that a wide class of problems in quantum mechan- 
ics can be easily solved with a CC. The problem of efficient state preparation is 
also considered. 

4.1 Quantum Entanglement and Quantum Informa- 
tion 

Here I present some examples of how the usual notion of QE can be exploited to 
perform interesting tasks in quantum Information. These tasks involve the well 
known processes of quantum cryptography and quantum teleportation. For sim- 
plicity, the setting used here usually involves two qubits (or many copies of it), A 
and B for Alice and Bob, which are initially prepared in some maximally entan- 
gled Bell State. 

4.1.1 Quantum Cryptography 

Quantum cryptography IIBB84I is the process that allows Alice and Bob to share 
a secret conversation through the encryption of their messages. Such an idea 
has been widely used in classical protocols. In a classical framework, the secret 
conversation between Alice and Bob is done by sending and receiving arrays of 
Os and Is. The encryption and decryption of these arrays are done by using a key 
which is assumed to be known only by them. This key is a random string of Os 
and Is and its length is usually the same that the length of the message to be sent. 
For example, assume that Alice wants to send Bob the binary number '10'. Then 
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Alice creates a random key, such as the array '11', and sends to Bob the encrypted 
message '10' + '11' = '01' (i.e., sum of binary numbers). Bob then receives the 
encrypted message and by using the key it proceeds to its decryption by adding 
them: 'Ol'+'ll' = '10', corresponding to the original message. This cryptography 
method is trivially extended to longer messages. 

The type of encryption described is more secure if the key used is changed in 
every sent message. Otherwise, a potential eavesdropper E (for Eve) could easily 
gain information about the key and therefore 'hear' the conversation between A 
and B. Nevertheless, more secure efficient methods exist in a quantum setting 
by using entangled states. The reason is simple: If Eve tries to figure out the 
key, which is now an entangled state, she needs to perform a measurement, and 
therefore destroys the entangled state by projecting it into some other unentangled 
one. 

In more detail, the scheme for quantum cryptography consists in a source that 
creates maximally entangled two-qubit Bell states |Bell) = ;^[|0^0b) + jl^lg)] 
(e.g., entangling photon polarizations, etc), where one qubit is given to Alice and 
the other is given to Bob (Fig. l4.1l) . In addition, each party can measure the state of 
the corresponding qubit in a different basis, if they want, from the logical one. (As 
mentioned in Chap.|3l such freedom is necessary to take advantage of the quantum 
correlations.) The source then emits many copies of the Bell state and the parties 
measure their corresponding qubits several times and different orientations. After 
the measurements are performed, Alice and Bob have a classical conversation 
where they teli each other about the chosen orientations but they never talk about 
the results of the measurements. Because of the nature of the Bell state, whenever 
they choose the same orientation, the results of the corresponding measurements 
are the same. For example, if Alice measures 1 1^) in the logical basis and Bob 
measures in the same basis, Alice knows that Bob measured the qubit state |le). 
In this way, they only keep the results obtained whenever the measurement is 
performed in the same orientation. Such results constitute the key, which was 
never discussed between them, for the encryption and decryption of the message. 

Although the scheme presented seems simple and successful, an eavesdropper 
could still hear the conversation by, for example, stealing entangled qubits used to 
build up the key. Many different cryptography protocols have been designed for 
different situations and can be found in the literature. Here I only discussed the 
bàsics of quantum cryptography to show that entangled states play an important 
role in information processes. 
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Figure 4. 1 : Quantum cryptographic protocol. A source L emits pairs of maximally 
entangled photons. Alice and Bob measure their polarization in random bases and 
they only keep those measured in the same direction. In this way they build a key 
used to encrypt their message. 

4.1.2 Quantum Teleportation 

Quantum teleportation IIBBC93II consists of the process of teleporting the state of 
a distant quantum system by means of local quantum operations, including mea- 
surements, and classical communication (LOCC). The simplest case considers the 
teleportation of the state of a single qubit. For this purpose assume that Alice pos- 
sesses two qubits denoted by Ai and A2, and Bob possesses only one, denoted 
by B. Moreover, assume that qubits A2 and B are in a maximally entangled state 
which was given by an external source. Such an entangled state can then be used 
to teleport the state of qubit Ai to qubit B. 

The idea is the foUowing: The global three-qubit state of Alice and Bob to- 
gether is given by = (a|0^,) + h\U,)) ® ^(|0^,0b) + lU^lfí))- This state 
can also be written as 




(4.1) 




1 



(|0^,0^,)±|UU))and 



(4.2) 



V2 



4.2 Quantum Entanglement and Quantum Computation 



87 



1 



(|0^,U,)±|UO^,)), 



V2 



respectively. Remarkably, in this basis (Eq. l4.2l) the state of Bob's qubit looks very 
similar to the state to be teleported (i.e., the state of ^i). Then, Alice performs a 
measurement in her two qubits projecting them into one of the four possible Bell 
States. (Such a measurement can be done using local operations only but I do not 
describe it here.). 

The state obtained by Alice corresponds to one of the four possibilities shown 
in Eq. 14.21 Alice then contacts Bob and telis him about the result of the measure- 
ment, after which Bob acts on its qubit to recover the state of Ai. For example, if 
Alice projects its two-qubit state into |,^+), then she telis this result (by means of a 
classical communication) to Bob. He thereby transforms the state of his qubit by 
using a flip operation transforming 



which is the desired teleported state. Similar operations can be done for any other 
state measured by Alice. 

4.2 Quantum Entanglement and Quantum Compu- 
tation 

Perhaps the most important practical case where QE can be used as a resource is 
in computational tasks. As mentioned before, parallelism is one of the properties 
of the quantum world that needs to be exploited in a quantum computation, and 
such a property is naturally associated to QE. Shor's factoring algorithm IISho94ll 
constitutes a nice example where entangled states are used to find the prime fac- 
tors of a whole number. It has been shown that using a QC, this algorithm can 
be efficiently performed; that is, with a number of steps that scales at most poly- 
nomially with the integer P to be factorized. However, it is not known how to 
perform such an algorithm efficiently on a CC: In order to find the factors of P it 
is necessary to divide it by all the whole numbers between 1, and P/2, and so on, 
constituting a time consuming task. 

To demonstrate Shor's algorithm, I present the factorization of the number 
P = 15. (Although the solution to this problem is immediately found, P = 15 = 
3x5, this method can be extended to the factorization of larger numbers like 
P' = 12319, where finding the solution is more complicated: P' = 97 x 127.) 



a\lB) + b\OB) 



a\OB) + b\lB), 



(4.3) 
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As I will show, the problem of finding the factors of P is equivalent to the order- 
finding problem. Thus, a random whole number m between 1 and P — 1 needs 
to be chosen first. Assume that this number is m = 8. Second, a QA needs to 
be performed to find the order n of m, modulo P. Such order is defined as the 
least positive integer such that mP" = 1( mod P), where a = b{ mod P) if a — 6 is 
divisible by P. 

Then, the first step of the QA consists of performing r Hadamard gates to the 
initial (r + 4)-qubit polarized state |0i ■ ■ ■ 0^) ® \OaObOcOd)- The number r of extra 
qubits has to be big enough to find the order n with high accuracy IISho94l . Here, 
r = 11. The other 4 qubits are necessary to encode Information about the factors 
of 15(i.e., 2^ = 16 > 15). 

After the Hadamard gates have been applied (i.e., the gates that transform 
|0) \+)), the evolved state is 




(4.4) 



where |j) is the state corresponding to the binary decomposition of the integer j 
in the logical basis; for example, \j = 3) = |0i ■ ■ ■ l^-iOr). The next step of the 
algorithm is to apply a unitary operation that transforms 

\j) ® \OaObOcOd) ^ |j) ® mod P)). (4.5) 

This operation can be efficiently performed using controUed operations on the 
state of the r qubits, but I will not explain it here. After this evolution, the evolved 
state is 

= —= \j) ® \m^{ mod P)) (4.6) 

= 11) + 11) ® |8) + |2) ® |4) + |3) ® |2) + ■ ■ ■] , 

where again I have chosen the integer representation of the states in the logical 
basis. The state of Eq. 14. 61 is a highly entangled state between the r = 11 qubits 
and qubits a, b, c, and d. 

By defining the states for qubits a, ò, c, d, as 

m^XmodP)), (4.7) 



^ Vexp 



—i2'7iqk 



n 



4.2 Quantum Entanglement and Quantum Computation 



89 



where again n is the not yet known order, Eq. l4.6l can be written as 

n-l 2^-1 



^ n-l T--\ 

exp 



q=0 j=0 



.2nqj 



n 



(4.8) 



A (inverse) quantum Fourier transform INCOOII performed to the r qubits in the 
above state maps it to 



E 

j=Q 



exp 



.2vrgj 



n 



\j) = q/n), 



(4.9) 



where the state \K = q/n) of the r qubits is a product state which depends on the 
coefficient q/n as in the usual Fourier transform. This transformation can also 
be efficiently implemented in a quantum circuit. Therefore, the transformed state 
reads 

^ n— 1 

\^) = —Y^ \K = q/n) ® I/i,). (4.10) 

* q=ü 

In this way, the register of r qubits is measured in the logical basis to project the 
state \il)) and obtain one possible state \K = q/n)® After applying the same 
algorithm several times and by applying a continued-fractions algorithm, it is pos- 
sible to obtain the order n from the statistics of the measured states \K = q/n). 
In this case one would obtain n = A, such that 8"^ = 4096 = 1( mod 15). Such 
a solution could be obtained if the result of the measurement in the register of 
r qubits gives, for example, \K = q/n) = |1536) (in the binary representation), 
where I have used again the integer representation of states in the logical basis. If 
this is the case, q/n = = |, resulting in n = 4. From this Information, one 
can obtain that 15 = 3 x 5 as it is explained in Ref. IISho94l . 

A factoring QA such as the one presented is more efficient than a classical 
one only because entangled states have been created at some step. Otherwise, no 
advantages can be gained, in general, from using a deterministic QC. As a simple 
example, suppose (with no loss of generality) that the initial pure state of a register 
of qubits is the usual polarized state |0i ■ ■ ■ Oat), and that the QA never does 
create QE. That is, the set of unitary gates Uj applied to the initial state are only 
single qubit rotations transforming product states into product states. Therefore, 

Uj = e'^^'^-ü) =cose + i sin Oaf,^.}^ , (4.11) 
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where 7r(j) = [1 ■ ■ ■ A^] denotes the qubit rotated at the jth step, and a{j) = 
[x,y,z] is the axis of rotation. The evolved state is then 

llUj\Oi---ON). (4.12) 

After the evolution is performed, certain qubit(s) is(are) measured in, for example, 
the logical basis. As usual, this is the last step of a deterministic QA (Sec. 12.21) . If 
the A;th qubit is measured, the result obtained is 

(^') = (Oi ■ ■ ■ ^i) ■ ■ ■ 0^)- (4.13) 

The fact is that such QA can be simulated classically and the expectation value 
of Eq. HTT^ can be computed in a CC with the same amount of resources (i.e., same 
efficiency). The idea is to only keep track of the Information about [9j, a{j)] if 
7r(j) = k, because the other qubits are not measured. If the QA is performed 
efficiently, only a polynomially large number of pairs [9j, a{j)] are kept. A clas- 
sical algorithm to obtain the result of Eq. l4.13l would consist of updating at each 
step the State of the A;th qubit only. For example, if at the jth step the state is 
Oj|0fc) + bj\lf:) and the gate f/j+i = 6*^°^== is performed after, the evolved state 
becomes 

[cos(7r/4) + ism{7r/4)a^]{aj\0k) + bj\lk)) = a^+ilOk) + (4.14) 

and since is the corresponding flip operation, cij^i = :^(% + ibj) and bj^i = 
-^{idj + bj). Updating the vàlues of the complex coefficients aj and bj can be 
easily done with a conventional computer. At the end of the total evolution the 
state of the kth qubit is aM\Ok) + ^A/llfc)- Therefore, the result of Eq. l4T3l is 

{a':) = \aM\'-\bM\', (4.15) 

which can be efficiently computed with a CC. 

However, if some of the gates Uj create entanglement, there is no known way, 
in general, to classically evaluate Eq. I4.13l efficientlv. Naturally, the elementary 
two-qubit gates, like the R^j^k{uj) or the CNOT gate introduced in Sec. 12.1.31 
are responsible for creating entanglement. For example, if the initial state is the 
product (unentangled) state ^|0 + 1) ® |0), applying the CNOT gate evolves it 

into ;^(|00) + | 11)), which is known to be one of the Bell states, that is, maximally 
entangled. Nonetheless, I will show in the next section that creation of QE is not a 
sufficient condition to claim that Eq. l4.13l cannot be obtained efficiently by using 
a classical algorithm. 
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4.3 Efficient Classical Simulations of Quantum Physics 

The purpose of this section is to generalize the concepts described previously, 
presenting a wide set of problems which can be solved efficiently with both, a 
QC or a classical one. To show this, I will exploit the notion of GE, which was 
presented in Sec. 13.21 Again, an algorithm is said to be efficient with respect to 
a certain variable N, whenever the amount of resources required to perform it is 
bounded by a polynomial in N. Otherwise, the algorithm is said to be inefficient. 
It is important to remark that here I only make a comparison between efficient 
and inefficient algorithms, with the previous definition, and do not compare the 
total amount of resources needed. For example, if a QA can be performed in 
operations and the corresponding classical one in A^^ operations then, despite 
being more convenient to use the former, both are considered efficient. 

In most cases, the type of problems one encounters in quantum mechanics are 
related with the evaluation of a certain expectation value or correlation function 

(W) = (4.16) 

where W is an operator acting on a finite quantum system S in some state \ (p) E H. 
(Ti. is the Hilbert space associated to S.) In Sec. 12.21 I presented some QAs that 
allow one to compute Eq. 14.161 in a QC described by the conventional model, 
whenever the operator àlgebra associated to the system S could be mapped onto 
the Pauli operators (Sec. 12.31) . The main steps of those QAs consist of the prepa- 
ration of the state \ip) = |+)a ® from some initial (boot-up) reference state 
like |0)a ® |0 ■ ■ ■ 0), then an evolution, and finally the measurement of the ancilla 
qubit a. Similarly, in the rest of this section I assume that the àlgebra of operators 
used to describe the quantum computational model and/or the àlgebra associated 
to the system S are not necessarily described by the conventional model, but a 
one-to-one efficient mapping between both exists. 

Since I am mainly interested in comparing the efficiencies of obtaining Eq. 14.161 
using a QC or a classical (conventional) one, the efficiency is defined relative to 
a certain variable A^ which here is, in general, the number of quantum mechan- 
ical elements that compose the system S. I assume that there exists a QA that 
evaluates Eq. 14. 161 efficientlv. that is, using at most poly(A^) resources. These are 
the so-called bounded polynomial quantum algorithms (BPQA). This is usually 
vàlid whenever the QC can be efficiently initialized and every step, including the 
preparation of can be efficiently performed. These results were previously 
discussed in Sec. 12.51 
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The most common classical algorithm to evaluate Eq. l4.16l consists of writing 
the matrix representation of W, and obtaining the corresponding matrix element. 
This is usually a hard task because of the dimension of the Hilbert space H grow- 
ing exponentially with A^. However, in some cases this complexity can be highly 
reduced. The simplest case is when the state |(^) of Eq. 14. 161 is known to be an 
eigenstate of W, so that the corresponding expectation is simply the known eigen- 
value. 

Nonetheless, a more general result can be obtained in a Lie algebraic frame- 
work. For this purpose, I start by pointing out that, with no loss of generality, the 
state 10) G H of Eq. 14. 161 can be obtained by transforming any reference state 
I ref ) G H; that is 

|0) = [/|ref), (4.17) 

where W = (unitary). Moreover, the reference state can be chosen to be 
the highest (or lowest) weight state | HW) of certain finite semi-simple Lie àlgebra 
t) = {Òi, Ó2, ■ ■ ■ , Òm}, with Òj = Ó] (Sec. 13.2.21) . In this case, f) also admits 
a Cartan-Weyl decomposition: i) = {Í)d, ^-}, with Í)d the Cartan subalgebra 
(CSA) of f), and f)+ and the set of raising and lowering operators, respectively 
(Sec. IXO) . 

Now suppose that the transformation U of Eq. 14.171 is a group operation in- 
duced by [); that is 



U = exp 



M 



(4.18) 



Then l'í/;) is a GCS and generalized unentangled state (extremal) of [). Moreover, 
suppose that W of Eq. RTT^ can be decomposed as 



M 



W = J2KjÒj, (4.19) 

i=i 

where kj G C Vj, and Oj G (). In general, W does not have the form of Eq. l4~T9l 
and it belongs to a Lie àlgebra other than [). Then, f) must be considered as the 
Lie àlgebra that contains both. (The important case when W E H = e*'', i.e. the 
group induced by í), is discussed in the next section.) 

In this way, computing classically Eq. l4.16l is equivalent to evaluating 

(ly) = ^K,(HW|e Oje^=' |HW). (4.20) 
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The operators 

i'=i 

are also elements in f) because they are transformed by a group operation. The co- 
efficients v^y of Eq. l4.21l can be computed by means of the representation theory. 
A matrix representation of \) is the mapping $ : f) — C^^ such that 

$(Ò,) = 0„ (4.22) 

mò,.òr]) = [<í>(ò,),<í>(ó,o] = [o,,o,v], 

where [A, B] is the commutator and the complex matrices Oj are p x p dimen- 
sional. A representation $ is said to hefaithful when maps linearly independent 
elements in f) to linearly independent p x p matrices. If the Lie àlgebra f) is com- 
pact, one can always find a faithful representation satisfying the inner product 

Tr[0,0,,] = 6,,,, (4.23) 

with Jr[Q] being the trace of Q. For example, the well known adjoint represen- 
tation (AppendixO, which associates an M x M matrix to every operator in f) 
(M is the dimension of {)), is faithful and satisfies Eq. l4.23[ Nevertheless, smaller 
dimensional faithful representations (p < M) can usually be found. Then, if $ is 
faithful, 

0; = e 0,e-^ =Y^u^^,0,,, (4.24) 

i'=i 

which is the matrix representation of Eq. I4.21[ Since this is a matrix operation, 
the matrices Ò'^ of Eq. l4!24l can be classically computed with C(p^) computational 
operations^ (addition and multiplication of complex numbers). Each coefficient 
Ujji is then classically obtained by using the inner product equation, that is 



n 



MO'prl (4.25) 



which also requires 0{p'^) computational operations^. 
'Here O denotes the order of required operations. 

^To obtain Vjji at accuracy e the order of poly (1/e) operations are needed. Nevertheless, in this 
section I consider that the calculation is done at the computational accuracy given by the particular 
type of CC used. 
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Similarly, in the CW basis, 

r l 

O'j = -fjkhk + hj'^aj' + '*jj'^-o^j' (4-26) 
fc=i j'=i 

where hk e Í)d (CSA), E^^, G [)+, and E^a^, G f)_ (Sec. l3.2.2l) . The coefficients 
"fjk and tjj/ of Eq. l4.26l can also be obtained classically with 0{p'^) computational 
operations by the trace inner product: The operators hk, E^.,, and E_a^, are usu- 

ally a simple linear combinations of the operators Oj. 

By definition, | HW) is an eigenstate of the operators in the CSA, that is 

/ifc|HW) = efc|HW), (4.27) 

where the eigenvalues are assumed to be known. Then, 

(HW|Í«^,,|HW) = (HW|È_„^,,|HW) = V/. (4.28) 

In other words, Eq. l4.20l can be classically computed by only keeping the elements 
of Ò] (Eq.i2D belonging to fi^: 

r r 

(HW|Ò;.|HW) = 5^7,fc(HW|/ifc|HW) = Y^^^^Ck, (4.29) 

k=l k=l 

and 

M r 

(^) = ZE'^i^^·'^^'=· (4.30) 

i=i k=i 

Since the coefficients Ck are assumed to be known, and the coefficients kj and jjk 
can be classically obtained with C(]?^), Eq. l4.30l can be computed with C(Mrp^) 
computational operations. If M < poly(A^), then r < poly(A^). Also, the exis- 
tence of the adjoint representation (Appendix O guarantees that one can always 
find p such that p < M < poly(A^). Therefore, Eq. l4.30l can be efficiently com- 
puted with a conventional computer. A more detailed proof about the number 
of operations required, including the approximation of the exponential matrix, is 
given in the Appendix|E| 

In brief, when the state |0) of Eq. 14.161 is a GCS of a polynomially large 
dimensional Lie àlgebra t), and the operator W is an element of f), Eq. l4.16l can be 
efficiently computed with both, a classical and a quantum device. 
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However, in some cases W or U (Eq. 14.171 1 are associated with an exponen- 
tially large dimensional Lie àlgebra and the corresponding classical simulation is 
inefficient. Nevertheless, it is important to remark that for mixed states of the form 

L 

P = '^Psps, (4.31) 

s=l 

where ps = |0s)(0s|, and \(ps) being GCSs and generalized unentangled states of 
a polynomially large dimensional Lie àlgebra, the computation of 

(W) = Tr(piy), (4.32) 

can also be efficiently performed with a CC if L < poly(A^). 



4.3.1 Higher Order Correlations 

The previous results can be extended to the computation of higher order correla- 
tions. For simplicity, I start by studying the case of two-body correlations of the 
form 

m = (0|W^i#2|0), (4.33) 

where W is the product of the observables Wi = Wl and W2 = Wl- Again, Wi 
and 14^2 are considered to be elements of a Lie àlgebra [) = {Oi, ■ ■ ■ , Om}, and 
10) is assumed to be a GCS of f); that is 



with U = exp 
Eq. l4.33l reads 



|0) = t/|HW), (4.34) 
i Y^- (jOj , and |HW) the highest weight state of f). In this way, 

{(t)\WiW2\(t)) = {HV\l\U^WiW2U\H\N). (4.35) 
Since UW = 1, 

{(l)\WiW2\(l)) = {l·\\N\W{W^\H\N), (4.36) 

where Wl = WWiU. 

Equation 14.331 can be efficiently computed with a classical computer if the 
dimension M of i) scales at most polynomially with the variable N; i.e. M < 
poly(A^). To show this, I first decompose the corresponding operators as 

M 

Wi = J2'^'A^ (4.37) 
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M 



(4.38) 



which transfomi under the action of U as 



M 



(4.39) 



In a CW basis of f), the transformed observables Ò'j = WÒjU can be decom- 
posed as in Eq. 14.261 Such decomposition can be done efficiently with a CC (i.e., 
computing the coefficients 7jfc and Ljji) by working in low-dimensional faithful 
representation of f) (AppendixO. Then, 



M 



with 



(HW|Ó:Ó;| HW)= E l^kljk'{hkh') + E in'ijfiÈa.È ) 

k,k'=l i',j'=l '- 



(4.40) 



CC. + Lii'L*j/{Ea.,E_a^,) + L*^, Ljji {E_a^, Ea.,) 



r l 



(4.41) 



+ Y. likLjj'ihkÈaJ + 'JikL*i>{hkÈ_o. ) + CC. 

k=lj'=i L j JJ 1 , 

where CC denotes the complex conjugate, and the expectation vàlues are now 
taken over the highest weight state, i.e. (A) = (HW|y4|HW). Since raising oper- 
ators annihilate | HW) only a few expectations do not vanish in Eq. 14.411 These 



are 



{hkhk') — GkCk', 

r r 



(4.42) 
(4.43) 



k=l 



k=l 



where the coefficients aj are the components of the root vector aj (Sec. 13.2.21) 
and Cfc are the known weights of | H W) . 
In brief, Eq. l4.33l takes the form 



\W,W2\ 



M 



4 



k'^k^k' 



k.k'=l 



ji=i k=i 



Lij't'jj'Oíj,ek 



(4.44) 
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which can be efficiently computed with a CC if M < poly(A^). 
Higher order correlation functions of the form 

{W) = {(t)\W'P ■■■W^\(t)) (4.45) 

can also be computed efficiently for a fixed integer R, whenever there is a polyno- 
mially large dimensional Lie àlgebra [) associated to the problem. Obviously, the 
complexity of this problem increases as p does. Again, the idea is to keep only 
the nonzero expectations of the transformed operator W' = WW^ ■ ■ ■ W^U = 
W'p . . -W'^, where |0) = ?7|HW). For this purpose, all the raising operators 
Éaj appearing in W' must be taken to the right side, by using the commutation 
relations of the àlgebra, so that they destroy the state |HW). This problem is 
analogous to the one given by Wick's theorem HWicSOl for fermionic operators. 
In Appendix |F1 I describe an efficient method for the classical computation of 
Eq.lífel 

4.4 Generalized Entanglement and Quantum Com- 
putation 

In Sec. 14.31 1 presented a class of classical algorithms that allows one to efficiently 
obtain certain correlation functions of quantum systems in a conventional Com- 
puter. Here, I analyze some advantages of having a QC. For example, assume 
then that one is interested in the evaluation of the correlation function of Eq. H~TH 
when there is no polynomially large dimensional Lie àlgebra f) associated with 
the problem. That is, either \(p) is not a GCS of [) or VÍ^ is not an element or 
group element induced by f) . Then, if a bounded polynomial quantum algorithm 
(BPQA) like the ones introduced in Sec. 12.21 exists in this case, Eq. 14. 161 can be 
efficiently computed using a QC. Usually, no known classical algorithm can be 
used to compute such a correlation efficiently and this represents, in some cases, 
an exponential speed-up of the quantum simulation with respect to the classical 
one. 

In Sec. 12.51 I discussed the complexity of simulating physical systems us- 
ing QCs. Although certain QSs have been shown to be efficient, others, like the 
obtention of the ground state properties in the two-dimensional Hubbard model, 
remained inefficient. The problem relied in the complexity of preparing the de- 
sired initial state When simulating quantum systems, this state is usually 
not completely known (e.g., it is the ground state of some non-solvable Hamilto- 
nian) and an approximated state, having a non-zero overlap with needs to be 
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prepared. The purpose of the foUowing section is then to show a wide class of 
quantum states that can be prepared efficiently by quantum networks. Again, the 
resuhs obtained are independent of the physical representation of the QC when- 
ever an efficient mapping between the operators used to describe the system to be 
simulated and the operators describing the computational model, exists. 

4.4.1 Efficient Initial State Preparation 

I start by describing the simple case when the initial state to be prepared is given 
by 

= f/|HW), (4.46) 

where |HW) is the highest weight of some Lie àlgebra f), and is considered to 
be the boot-up state of the computer, which can be efficiently initialized. If the 
operator U of Eq. 14.461 is a group operation induced by f|, that is, U = e'^ and 
A G f) = {Oi, ■ ■ ■ , Om}, and M < poly(A^), then |0) can be prepared efficiently 
with arbitrary accuracy by using, for example, a first order Trotter decomposi- 
tion IISuz93ll . In this case, the state \(p) is generalized unentangled and GCS of t). 
Then, À = ^^li CjÓj and 

gíA^-Qg^Aí^ (4.47) 

M 

Assuming that every operation e*^^'^^^* is either an elementary gate or can be 
performed by applying a small number of them, the evolution e*"^ can be then 
approximated by applying a polynomially large number of these gates. Particular 
examples were discussed in Sec. l2.3[ 

Another interesting case is when the state is not given in the form of 
Eq. 14.461 but some Information, like the expectation vàlues of a set of observ- 
ables (i.e., its reduced state), {Oj), is known or can be easily obtained. If f) = 
{Oi, ■ ■ ■ , Om} is a semi-simple Lie àlgebra, with M < poly(A^), and |0) is shown 
to be a GCS of f) by calculating, for example, the í)-purity (Sec. 13.2.2b . then it can 
be prepared efficiently with a QC. To show this, one needs to obtain first the op- 
erator A = X],=i CjOj such that U = e*^ is the transformation of Eq. l446l 
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For this purpose, I define the fictitious Hamiltonian 



M 



HF = -J2{Ój)Ó„ (4.49) 
i=i 

whose expectation value over the state {ip) is 

{Í;\Hf\^) = -Y.(Òj)'- (4.50) 

i=i 

If = Vlíp) is another GCS of f), with V a unitary group operation that trans- 
forms as 

M 

V^Ò,V = Y,'^n'Òr, (4.51) 

j'=i 

then 

M 

{^P\HF\^P) = - ^ (Ò,>,y(Ò,.). (4.52) 

Since the M x M matrix defined by the real coefficients Ujji is orthogonal, one 
obtains (Schur's inequality) 

(^I^fI^) > (mpl^)- (4.53) 

Therefore, \ip) is the unique GCS and ground state that minimizes such expecta- 
tion value. 

Assume now that Hp can be efficiently diagonalized with a conventional com- 
puter, obtaining 

Hp = UHoU\ (4.54) 
where U isa unitary group operation and 

r 

Hq = Y^ Ikhk , (4.55) 

k=l 

is an unperturbed Hamiltonian (i.e., are elements in the CSA of [)). Such a 
diagonalization algorithm is described in Sec. 15.2.11 Defining \íp) = one 
obtains 

impi^) = imoi^), (4.56) 
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which takes its minimum value, by definition, for a certain weight state In 
particular, the coefficients 7^ of Eq. H35l can be chosen such that = | HW) (see 
Chap.|5ll. Therefore, = f/|HW), and U is the desired group operation. Again, 
U can be approximated by applying a polynomially large number of elementary 
gates so can be prepared efficiently with a QC. 

4.5 Summary 

In this chapter I have addressed several issues relating the concept of GE with 
quantum and classical complexity. I have shown that the creation of GE relative 
to every polynomially large dimensional Lie àlgebra is a necessary requirement to 
gain efficiency over the known classical algorithms. Otherwise, efficient classical 
simulations exist, which realize the computation by working in low-dimensional 
representations of the àlgebra. This result goes beyond the requirement of creating 
entangled states in the usual sense. For example, a quantum algorithm involving 
gates induced by the Lie àlgebra [) = so(2A^), which is presented in detall in the 
next chapter, creates entangled states in the Standard sense (i.e., non-separable). 
However, since the dimension of f) is polynomially large, such quantum computa- 
tion can be efficiently simulated with a CC. 

The efficient preparation of a wide class of states by quantum networks has 
also been studied. Again, any GCS or generalized unentangled state of a polyno- 
mially large dimensional Lie àlgebra can be efficiently prepared on a QC. 



Chapter 5 



Generalized Entanglement and 
Many-Body Physics 

Information is inevitably physical. 
Rolf Landauer 

Every QC has associated certain physical representation (e.g., spin- 1/2 system, 
etc). The fact that Information is not just an abstract entity and is always linked 
to a physical system, implies that it must be govemed by the laws of physics. 
In Quantum Infomiation Theory, for example, the manipulation and control of 
Information is based on the foundations and laws of the quantum mechanics. On 
the other hand, many concepts, such as QE, have been developed only from an 
information-theory point of view. It is expected then that these concepts are also 
useful for the analysis of physical phenomena. In this chapter, I apply the notion 
of GE (Chap.|31) to the study of different problems in many-body physics. 

First, I will focus on the characterization of quantum phase transitions (QPTs) 
in matter. QPTs are qualitative changes occurring in the properties of the ground 
State of a many-body system due to modifications either in the interactions among 
its constituents or in their interactions with an extemal probe IISac99l . while the 
system remains at zero temperature. Typically, such changes are induced as a pa- 
rameter g in the system Hamiltonian H{g) is varied across a point at which the 
transition is made from one quantum phase to a different one. Often some cor- 
relation length diverges at this point, in which case the latter is called a quantum 
critical point. Because thermal fluctuations are inhibited, QPTs are purely driven 
by quantum fluctuations. Thus, these are purely quantum phenomena: A classical 
system in a pure state cannot exhibit quantum correlations. Prominent examples 
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of QPTs are the quantum paramagnet to ferromagnet transition occurring in Ising 
spin systems under an extemal transversa magnètic field IILSM6 1 1 |Pfe701 |BM7 1 1 . 
the superconductor to insulator transition in high-Tc superconducting systems, and 
the superfluid to Mott insulator transition originally predicted for liquid helium 
and recently observed in ultracold atòmic gases |GMEQ2|. 

Since entanglement is a property inherent to quantum states and intimately 
related to quantum correlations (Chap.El), one would expect that, in some appro- 
priately defined sense, the entanglement present in the ground state undergoes a 
substantial change across a point where a QPT occurs. Thus, the concept of GE 
becomes specially well suited for this study because it is directly applicable to any 
algebraic language associated to the system under study. 

Another important problem in quantum mechanics is to exactly diagonalize 
and obtain the spectrum of a Hamiltonian of a many-body system. In this case, no 
approximate methods (e.g., mean field theory) are needed and one has complete 
knowledge of the physical properties of the system through algebraic methods. 
By using information-theory concepts such as the ones described for the efficient 
simulation of physical systems (Sec. 14.31) . I will show that whenever there is a Lie 
algebraic structure behind a problem in quantum physics, and the dimension of 
the associated Lie àlgebra is small enough, such problems can be solved easily 
by using a CC. This constitutes the final part of my thesis after which come the 
conclusions (Chap.0). 

5.1 Entanglement and Quantum Phase Transitions 

In this section, I characterize the QPTs present in the Lipkin-Meshkov-Glick 
(LMG) model and in the anisotropic XY model in an extemal magnètic field 
through the GE notion, relative to a particular subset of observables which will 
be appropriately chosen in each case. Interestingly, for both of these models the 
ground states can be computed exactly by mapping the set of observable opera- 
tors involved in the system Hamiltonian to a new set of operators which satisfy the 
same commutation relations; thus, preserving the underlying algebraic structure. 
In the new operator language, the models are seen to contain some symmetries 
that make them exactly solvable, allowing one to obtain the ground state proper- 
ties in a number of operations that scales polynomially with the system size. It is 
possible then to understand which quantum correlations give rise to the QPTs in 
these cases. 

Several issues should be considered when looking for an àlgebra [) of observ- 
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ables that may make the corresponding relative purity a good indicator of a QPT. 
First, one observes that in each of these cases a preferred Lie àlgebra exists, where 
the respective ground state would have maximum -purity independently of the 
interaction strengths in the Hamiltonian. The purity relative to such an àlgebra re- 
mains constant, therefore it does not identify the QPT. (In these cases, this àlgebra 
is in fact the Lie àlgebra generated by the parametrized family of model Hamil- 
tonians, as the parameters are varied.) Thus, one needs to extract a subalgebra 
relative to which the ground state may be generalized entangled, depending on 
the parameters in the Hamiltonian. A second, closely related observation is that 
the purity must contain Information about quantum correlations which undergo a 
qualitative change as the transition point is crossed: thus, the corresponding de- 
gree of entanglement, as measured by the purity, must depend on the interaction 
strengths goveming the phase transition. Finally, whenever a degeneracy of the 
ground state exists or emerges in the thermodynamic limit, a physical requirement 
is that the purity be the same for all ground states. 

Although these restrictions together tum out to be sufflcient for choosing the 
relevant àlgebra of observables in the foUowing two models, they do not pro- 
vide an unambiguous answer when solving a non-integrable model whose exact 
ground state solution cannot be computed efficiently. Typically, in the latter cases 
the ground states are GCSs of Lie algebras for each of which the dimension in- 
creases exponentially with the system size. Choosing the observable subalgebra 
that contains the proper Information on the QPTs (such as Information on critical 
exponents) then becomes, in general, a difficult task. 

A concept of general mean-field Hamiltonians (GMFH) emerges from these 
considerations. Given a Hilbert space of dimension (with p an integer > 1), 
I define a GMFH as the Hermitian operator 

iÍMF = ^ l·íjÒj , /tj G M , (5.1) 

a 

which is an element of an irreducibly represented Lie àlgebra of Hermitian opera- 
tors f) = {Oi, ■ ■ ■ , Om}, whose dimension scales at most polynomially in N that 
is, M < poly(iV). When the ground state of í/mf is non-degenerate, it turns out to 
be a GCS of f) IIBKO03L while the remaining eigenstates (some of which may also 
be GCSs) and energies can be efficiently computed. The connection between Lie- 
algebraic mean-field Hamiltonians and their efficient solvability deserves careful 
analysis in its own right, and will be presented in the foUowing section. 
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5.1.1 Lipkin-Meshkov-Glick Model 

Originally introduced in the context of nuclear physics IILMG65I . the Lipkin- 
Meshkov-Glick (LMG) model is widely used as a testbed for studying critical phe- 
nomena in (pseudo) spin systems. This model was shown to be exactly-solvable 
in Ref. IIOSD05II . In this section, I investigate the critical properties of this model 
by calculating the purity relative to a particular subset of observables, which will 
be chosen by analyzing the classical behavior of the ground state of the system. 
For this purpose, I first need to map the model to a single spin, where it becomes 
solvable and where the Standard notion of entanglement is not immediately appli- 
cable. 

The model is constructed by considering N fermions distributed in two A^-fold 
degenerate levels (termed upper and lower shells). The latter are separated by an 
energy gap e, which will be set here equal to 1 . The quantum number a = ±1{] or 
[) labels the level while the quantum number k denotes the particular degenerate 
state in the level (for both shells, k G {ki, . . . .k^}). In addition, I consider a 
"monopole-monopole" interaction that scatters pairs of partides between the two 
levels without changing k. The model Hamiltonian may be written as 

ií = /7o + V + # = i^acLc,, + ^S^cLcL^^c,^ (5.2) 

k,<T k,k',a 

\- t t 

9/V 2-^ ^ku^k'a^k'a^^ka ' 



where a = —a, and the fermionic operators c]^^ (Cka) create (annihilate) a fermion 
in the level identified by the quantum numbers {k, a) and satisfy the fermionic 
commutation relations given in Sec. 12.3. li Thus, the interaction V scatters a pair 
of partides belonging to one of the levels, and the interaction W scatters a pair of 
partides belonging to different levels. Note that the factor 1 /N must be present in 
the interaction terms for stability reasons, as the energy per particle must be finite 
in the thermodynamic limit. 

Upon introducing the pseudospin operators 

k 

J- = E^c.T, (5.4) 
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= IjY^ ^^a^ka = 2 - ^l) > (5-5) 

fc,(T 

which satisfy the su(2) commutation relations of the angular momentum àlgebra, 

[J„J±] = ±J±, (5.6) 

[J+.J^] = Jz. (5.7) 

the Hamiltonian of Eq. l5.2l mav be rewritten as 

H = J. + ^{Jl + J!) + ^{J+J- + J-J+) ■ (5.8) 



As defined by Eq. 15.81 H is invariant under the Z2 inversion symmetry operation 
K that transforms {Jx, Jy, Jz) ^ i—Jx,—Jy,Jz), and it also commutes with the 
(Casimir) total angular momentum operator = + Jy + J"^. Therefore, the 
non-degenerate eigenstates of H are simultaneous eigenstates of both K and J^, 
and they may be obtained by diagonalizing matrices of dimension 2 J+ 1 (whereby 
the solubility of the model). Notice that, by definí tion of Jz as in Eq. 15.51 the 
maximum eigenvalue of and J = | J| is N/2. In particular, for a system with 
fermions as assumed, both the ground state \g) and the first excited state |e) belong 
to the largest possible angular momentum eigenvalue J = N/2 IILMG65II (so- 
called half-fiUing configurations); thus, they can be computed by diagonalizing a 
matrix of dimension + 1 . 

The Hamiltonian of Eq. l5.8l does not exhibit a QPT for finite N. It is important 
to remark that some critical properties of the LMG model in the thermodynamic 
limit ^ 00 can be understood by using a semiclassical approach (note that 
the critical behavior is essentially mean-field): first, I replace the angular momen- 
tum operators in H/N (with H given in Eq. 15.81) by their classical components 
(Fig.EI]); that is 

J = {Jx,Jy,Jz) {J sinO cos(j), J sinOsmcj), J cos9) , (5.9) 
H/N ^ h,ij,e,(j)), (5.10) 

where hc is the resulting classical Hamiltonian and j = J/N, j = 0, . . . , 1/2. In 
this way, one can show that in the thermodynamic limit (AppendixiGb 

lim = lim ^ = min hJj,e,(j)) , (5.11) 
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so the ground state energy per particle Eg/N can be easily evaluated by minimiz- 
ing 



Figure 5.1: Angular momentum coordinates in the three-dimensional space. 

As mentioned, the ground and first excited states have maximum angular mo- 
mentum i = 1/2. In Fig. 15.21 1 show the orientation of the angular momentum 
in the ground states of the classical Hamiltonian h^, represented by the vectors J, 
Ji, and J2, for different vàlues of V and W. When A = \V\ — W < 1 we have 
9 = 71 and the classical angular momentum is oriented in the negative 2-direction. 
However, when A > 1 we have cos 6* = — A^^ and the classical ground state 
becomes two-fold degenerate (notice that hc is invariant under the transformation 
I— i> —0). In this region and for V < the angular momentum is oriented in the 
XZ plane (0 = 0) while for \/ > it is oriented in the YZ plane (0 = ±7r/2). The 
model has a gauge symmetry in the line V = 0, W < —1, where can take any 
possible value. 

First and Second Order QPTs, and Critical Behavior 

For the Hamiltonian of Eq. 15.21 the quantum system undergoes a second order 
QPT at the critical boundary Ac = \Vc\ — Wc = 1, where for A > Ac the ground 
and first excited states \g) and |e) become degenerate in the thermodynamic limit, 
and the inversion symmetry K breaks. The order parameter is given by the mean 





J 



y 
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Figure 5.2: Representation of the classical ground state of the LMG model. 

number of fermions in the upper shell (rij) = 1/2 + ( J^)/iV, which in the ther- 
modynamic limit converges to its classical value, 

lim („,) = i±|2ïí . (5.13) 

Obviously, for A < Ac one has = 0, and {n^) > 0, otherwise (Fig. 15.21) . The 
critical exponents of the order parameter are easily computed by making a Taylor 
expansion near the critical points (A 1+). Defining the quantities x = Vc — V 
and y = Wc — W, one obtains 

A™+^^^ \ (í/" + x^)/2 forV<0 ' 

where the critical exponents are a = 1 and P = 1. 

In Fig. 15.31 1 show the exact ground state energy per particle Eg/N (with Eg = 
{g\H\g)) as a function of V and W in the thermodynamic limit (Eqs. 15.1 ÏT ). One 
can see that also in the broken symmetry region (A > 1) the system undergoes 
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a first order QPT at 1^ = 0; that is, the first derivative of the ground state energy 
with respect to V is not continuous in this line. 




Y z -n 



Figure 5.3: Ground state energy per particle in the LMG model. 



Purity as an Indicator of the QPTs in the LMG Model 

The Standard notion of entanglement is not directly applicable to the LMG model 
as described by Eq. 15.81 as this is a single spin system and no physically natural 
partition into subsystems is possible. Therefore, using the [)-purity as a measure 
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of entanglement becomes an advantage from this point of view, since the latter 
only depends on a particular subset of observables and no partition of the system 
is necessary. The first required step is the Identification of a relevant Lie àlgebra 
of observables relative to which the purity has to be calculated. 

Since both the ground and first excited states of the quantum LMG model may 
be understood as states of a system carrying total angular momentum J = N/2, a 
first possible àlgebra to consider is the su{N + 1) àlgebra acting on the relevant 
(N + l)-dimensional eigenspace. Relative to this àlgebra, \g) is generalized un- 
entangled for arbitrary vàlues of V and W, thus the corresponding purity remains 
constant and does not signal the QPTs. However, the family of Hamiltonians of 
Eq. l5.8l do not generate this Lie àlgebra, but rather an su(2) àlgebra. 

Thus a natural choice, suggested by the commutation relationships of Eqs. l5^ 
and l5.7l is to study the purity relative to the spin- representation of the angular 
momentum Lie àlgebra f) = su(2) = {Jx, Jy, Jz}' 



where the normalization factor K = N'^/A is chosen to ensure that the maximum 
of P(, is equal to 1 (Eq. l3.32l) . With this normalization factor, P[j can be calculated 
exactly in the thermodynamic limit by relying on the semi-classical approach de- 
scribed earlier (Appendix [ü] and Eq. 15.91) . For V = and arbitrary W > 0, 
Ig) = \Jz = —N/2) which is a GCS of su(2) and has P^ = 1. For genèric interac- 
tion vàlues such that A < 1, the classical angular momentum depicted in Fig. 15.21 
is oriented along the z direction and is not degenerate: Because ( J^) = {Jy) = 0, 
only (Jz) contributes to P^. By recalling that \imiy^oo{Jz/N) = —1/2, this gives 
P(, = 1, so that as far as relative purity is concerned the ground state behaves 
asymptotically like a coherent state in the thermodynamic limit. Physically, this 
means that with respect to the relevant fluctuations, GCSs of 5u(2) are a good 
approximation of the quantum ground state for large particle numbers, as is well 
established for this model [GF78|. However, in the region A > 1 the ground state 
(both classical and quantum) is two-fold degenerate in the — * oo limit, and the 
value of P(, depends in general on the particular linear combination of degenerate 
states. This can be understood from Fig. 15.21 where different linear combinations 
of the two degenerate vectors Ji and J2 imply different vàlues of {J^) for < 
and different vàlues of {Jy) for V > 0, while {Jz) remains constant. With these 
features, the purity relative to the su(2) àlgebra will not be a good indicator of the 
QPT. 

An altemative option is then to look at a subalgebra of su(2). In particular, if 



4 - 

ÏV2 . 



{J.? + {Jy)' + {Jz)' 



(5.14) 
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I only consider the purity relative to the single observable P) = so(2) = { J^} (i.e., 
a particular CSA of su(2)), and retain the same normalization as above, then 




(5.15) 



This new purity will be a good indicator of the QPT, since i^f, = 1 only for A < 1 
in the thermodynamic limit, and in addition P^, does not depend on the particular 
linear combination of the two-fold degenerate states in the region A > 1, where 
P(, < 1. Obviously, in this case Pf, is straightforwardly related to the order param- 
eter (Eq. l5.13t : the critical exponents of Pf, — 1 are indeed the same (a = 1 and 
/5 = 1). 

In the region A < 1 where P(, = 1, the quantum ground state of the LMG 
model (Eq. 15.81) does not have a well defined ^-component of angular momentum 
except at V = ([H, Jz] ^ if V ^ 0), thus in general it does not lie on a 
coherent orbit of this àlgebra for finite N. However, as discussed above, it behaves 
asymptotically (in the infinite N limit) as a GCS (in the sense that P(j 1). 

In Fig. 15.41 1 show the behavior of Pf, as a function of the parameters V and 
W. The purity relative to is then a good indicator not only of the second order 
QPT but also of the first order QPT (the line ^ = 0, 1^ < -1). 



5.1.2 Anisotropic XY Model in a Transverse Magnètic Field 

In this section, I exploit the purity relative to the u(A^) àlgebra (introduced in 
Sec. 13.3.31) as a measure of GE which is able to identify the paramagnetic-to- 
ferromagnetic QPT in the anisotropic one-dimensional spin-1/2 XY model in a 
transverse magnètic field and classify its universality properties. 

The model Hamiltonian for a chain of N sites is given by (Fig. 15.51) 

N N 

^ = E [(1 + + (1 - ^K<'] + E (5.16) 

1=1 1=1 

where the operators (a = x,y,z) are the Pauli spin-1/2 operators on site i 
(Sec. 12. II) . g is the parameter one may tune to drive the QPT, and < 7 < 1 is the 
amount of anisotropy in the XY plane. In particular, for 7=1, Eq. l5.16l reduces 
to the Ising model in a transverse magnètic field IIPfe70l . while for 7 — > the 
model becomes isotropic. Periòdic boundary condi tions are considered here, that 
is cr!+^ = cr* , for all i and a. 
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Figure 5.4: Purity relative to the observable in the ground state of the LMG 
model. 



When g ^ l and 7 = 1 the model is Ising-like. In this limit, the spin-spin 
interactions are the dominant contribution to the Hamiltonian (Eq. I5.16I) . and the 
ground state becomes degenerate in the thermodynamic limit, exhibiting ferro- 
magneticlong-range order correlations in the a; direction: = lira n 00 {o'icrx^'^) > 
0, where is the magnetization in the x direction. In the opposite limit where 
g ^ 0, the external magnètic field becomes important, the spins tend to align 
in the z direction, and the magnetization in the x direction vanishes: = 
\im N ^ {ala x^"^) = 0. Thus, in the thermodynamic limit the model is subject 
to a paramagnetic-to-ferromagnetic second order QPT at a critical point g^ that 
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X 

5^ 



N 



Figure 5.5: Anisotropic one-dimensional XY model in an external transverse mag- 
nètic field B. 



will be determined later, with critical behavior belonging to the two-dimensional 
Ising model universality class. 

This model can be exactly solved using the Jordan-Wigner transformation IIJW28I . 
which maps the Pauli (spin 1/2) àlgebra into the canonical fermion àlgebra through 
(Sec. ICTT) 



4 



n( 



(5.17) 



In order to find the exact ground state, I first need to write the Hamiltonian 
given in Eq. 15. 161 in terms of these fermionic operators, 



Af-l 



H = -2g Y,{4c.+i+lclcl^i + h.c.)+2gK{c^j^c,+^c^j^c\ + h.c.)+2N , (5.18) 



1=1 



N 



where -^"=11 {~'^z) operator that commutes with the Hamiltonian, and 



N 



N = '^i'^i is the total number operator (here, I choose to be even). Then, 



í=i 



the eigenvalue of is a good quantum number, and noticing that K = e*'^^, one 
obtains K = 1) whenever the (non-degenerate) eigenstate of ií is a linear 

combination of states with an even (odd) number of fermions. In particular, the 
numerical solution of this model in finite systems (with N even) indicates that the 
ground state has eigenvalue K = +1, implying anti-periodic boundary condi tions 
in Eq.l5ÏÏ8l 
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The second step is to rewrite the Hamiltonian in terms of the fermionic opera- 
tors cl (c^), defined by the Fourier transform of the operators cj (Cj): 



N 



LE^-^'M' (5.19) 



where the set V of possible k is determined by the anti-periodic boundary con- 

-ÍL -I-Stt _|_ (jV-l' 
■Af'-'-Af'''''-'- TV 



ditions in the fermionic operators: V = V"+ + VI = [±^, ±— ■ ■ ■ -^Í^O)ZLl 



Therefore, 

H + N = -2^(-l + 25(cos/c)4cfc + Í5'7sinA;(c^fc4 + c_fcCfc) • (5.20) 

keV 

The third and final step is to diagonalize the Hamiltonian of Eq. l5.20l using the 
Bogolubov canonical transformation [BRSó | 

7fc = UkC^ - WkC^_f^ 
where the real coefficients Uk and Vk satisfy the relations 



with 



Uk = u-k, Vk = -v^k , and ul + vl = 1 , (5.21) 



Uk = cos ( Y j , Vk = sin ) ' (5.22) 



and (pk given by 

, , , 2gf7 sin k 
tan(0fc) = — . (5.23) 

— 1 + 2g cos k 

In this way, the quasiparticle creation and annihilation operators and 7;., 
also satisfy the canonical fermionic anti-commutation relations of Eq. 12.341 and 
the Hamiltonian may be finally rewritten as 

H = J2^^(^hk-m, (5.24) 

keV 



where ^k = 2^ (—1 + 2g cos k)"^ + 4(7^7^ sin^ k is the quasiparticle energy. Since 
in general ^k > 0, the ground state is the quantum state with no quasiparticles 
(BCS state IITak99l ). such that 7fc| BCS) = 0. Thus, one finds 

I BCS) = II {uk + ^VkèiclJ^) I vac) , (5.25) 

keV+ 
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where |vac) is the state with no fermions (c^lvac) = 0). 

Excited States with an even number of fermions (K = +1) can be obtained 
applying pairs of quasiparticle creation operators -yl to the |BCS) state. However, 
one should be more rigorous when obtaining excited states with an odd number of 
partides, since K = —1 implies periòdic boundary conditions in Eq. 15.181 and the 
new set of possible k's (wave vectors) is ^ = [— tt, ■ ■ ■ , — 0, ^, ■ ■ ■ , ^^^^^^] 
(different from V). 

QPT and Critical Point 

In Fig. 15.61 1 show the order parameter = Mm.N^ooWlcrx^'^) as a function 
of g in the thermodynamic limit and for different anisotropies 7 IIBM71II . Then, 
= for g < Qc and 7^ for (7 > g^, so the critical point is located at 
gc = 1/2, regardless of the value of 7. The value of gc can also be obtained by 
setting = in Eq. 15.241 where the gap vanishes. 

1 I ' 1 ' 1 ' 1 ' 1 ' 

Ml 
0.8 - 

0.6 - 

0.4 - 

0.2 ■ 

- 
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Figure 5.6: Order parameter in the thermodynamic limit as a function of g for 
different anisotropies 7. The critical point is at g^. = 1/2. 

The Hamiltonian of Eq. l5.16l is invariant under the transformation that maps 
{(^U (^i'^ (^z) ^ (-<•; (^2 symmetry), implying that (<.) = for all 
g. However, since in the thermodynamic limit the ground state becomes two-fold 
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degenerate, for g > , itis possible to build up a ground state where the discrete 
Z2 symmetry is broken, i.e. (a*.) ^ 0. This statement can be easily understood 
if we consider the case of 7 = 1, where for < g < gc the ground state has no 
magnetization in the x direction: For g = 0, the spins align with the magnètic 
field, while an infinitesimal spin interaction disorders the system and = 0. 
On the other hand, for (7 ^ cxd the states \gi) = ^[|— í>, ■ ■ ■ , ^) + ■ ■ ■ , ^)] 



and \g,) = ^[\^,...,^)-\^,..., ^)], with H) = M\^) + li)] and |^) = 



7^[|t) " li)]' become degenerate in the thermodynamic limit, and a ground state 
with (al) ^ can be constructed from a linear combination. 

Remarkably, this paramagnetic-to-ferromagnetic QPT does not exist in the 
isotropic limit (7 = 0). In this case, the Hamiltonian of Eq. 15.161 has a con- 
tinuous u(l) symmetry; that is, it is invariant under any z rotation of the form 
exp[i6' cr^]. Since the model is one-dimensional, this symmetry cannot be 
spontaneously broken, regardless of the magnitude of the coupling constants. 
Nevertheless, a simple calculation of the ground state energy indicates a diver- 
gence in its second derivative at the critical point gc = 1/2, thus, a second order 
non-broken symmetry QPT. For g < g^ all the spins (in the ground state) are 
aligned with the extemal magnètic field, with total magnetization in the z direc- 
tion = J2j{'^z) — and the quantum phase is gapped. For g > g^., the total 
magnetization in the z direction is > — A^, the gap vanishes, and the quantum 
phase becomes critical (i.e., the spin-spin correlation functions decay with a power 
law), with an emergent u(l) gauge symmetry |BO04 |. Then, in terms of fermionic 
operators (Eq. l5.18l l. an insulator-metal (or superfluid) like second order QPT ex- 
ists at gc for the isotropic case, with no symmetry order parameter. It is a Lifshitz 
transition. 

u(A^)-Purity in the BCS State and Critical Behavior 

The |BCS) state of Eq. 15^251 is a GCS of the àlgebra of observables í) = so{2N). 
This is spanned by an orthonormal Hermitian basis which is constructed by ad- 
joining to the basis of u(A^), given in Eq. 13.531 the foUowing set r of number-non- 
conserving fermionic operators: 



Then, the |BCS) state is generalized unentangled with respect to the so(2iV) àl- 
gebra and its purity P[j contains no Information about the phase transition: Pf, = 
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1 V(7, 7. Therefore, in order to characterize the QPT one needs to look at the pos- 
sible subalgebras of so(2A^). A natural choice is to restrict to operators which 
preserve the total femiion number, that is, to consider the u(A^) àlgebra defined in 
Sec. 13.3.31 relative to which the |BCS) state may become generalized entangled. 
(Note that as mentioned in Sec. 13.3.31 the u(A^) àlgebra can also be written in 
terms of the fermionic operators c|, and Cfc, with k belonging to the set V.) 

In the |BCS) state, {c\ci^,) 7^ only if k = k', thus using Eq. 13.541 the purity 
relative to f) = u(A^) is: 

A(|BCS)) = 1 J](cK - 1/2)^ = A J2{vl - 1/2)2 (5_2V) 

keV keV 

where the coefficients can be obtained from Eqs. I5.22l and l5.231 In particular, 
for (7 = the spins are aligned with the magnètic field and the fuUy polarized 
|BCS)ç,=o = li) i) • • • ) i) state is generalized unentangled in this limit (a GCS of 
u(A^) with P(, = 1). In the thermodynamic limit, the purity relative to the u(A^) 
àlgebra can be obtained by integrating Eq. 15.271 

A(|BCS)) = ^ j{vl-l/2fdk, (5.28) 


leading to the foUowing result: 

n(iBcs)) = 



if(?< 1/2 



^l_4g2(l_^2)_ 

^ \íg>l/2 



(5.29) 



1+7 

Although this function is continuous, its derivative is not and has a dràstic change 
dX g = 1/2, where the QPT occurs. Moreover, P^^ is minimum for g > 1/2 
implying maximum entanglement at the transition point and in the ordered (ferro- 
magnetic) phase. Remarkably, for g > 1/2 and N ^ 00, where the ground state 
of the anisotropic XY model in a transverse magnètic field is two-fold degener- 
ate, P(, remains invariant for arbitrary linear combinations of the two degenerate 
States. 

As defined, for large g the purity P(j approaches a constant value which de- 
pends on 7. It is convenient to remove such dependence in the ordered phase by 
introducing a new quantity = Pt, — (shifted purity). Thus, 
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The new function behaves like a disorder parameter for the system, being zero 
in the ferromagnetic (ordered) phase and different from zero in the paramagnetic 
(disordered) one. The behavior of P,^ as a function of g in the thermodynamic 
limit is depicted in Fig. 15.71 for different vàlues of 7. In the special case of the 
Ising model in a transverse magnètic field (7 = 1), one has the simple behavior 
= 1/2 - 2g^ for g < 1/2 and P^ = if ^ > 1/2. 




Figure 5.7: Shifted purity PÚ(Ar) of the |BCS) as a function of g for different 
anisotropies 7, Eq. 15.301 PÚ[n) behaves like a disorder parameter for this model, 
sharply identifying the QPT at = 1/2- 

The critical behavior of the system is characterized by a power-law divergence 
of the correlation length e, which is defined such that íox g < 1/2, 

lim |(aX)|^exp(-t^). (5.31) 

Thus, e ^ 00 signals the emergence of long-range correlations in the ordered 
region g > 1/2. Near the critical point (g 1/2^), the correlation length behaves 
as e ~ {gc—g)"", where is a critical exponent and the value u = \ corresponds to 
the Ising universality class. Let the parameter A2 = e^^/^ The fact that the purity 
contains Information about the critical properties of the model foUows from the 
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possibility of expressing P/ for (7 < 1/2 as a function of the correlation length, 



^;(iBcs)) 



7 



7^ 



1 + 



7 



(5.32) 



2^A2(l-7)-l. 

where a known relation between g, 7, and A2 has been exploited IIBM71I . Per- 
forming a Taylor expansion of Eq. 15.321 in the region g l/2~, one obtains 

(Fig.EU) 

P;~2((7,-(?)77 z/ = l,7>0. (5.33) 
Thus, the name disorder parameter for is consistent. 



y = a + iAii\g - g,\ ^J-/ 




V = 1 



-2 -1 

1" Ifl - flcl 



Figure5.8: Scaling properties of the disorder parameter for anisotropy 7 = 1. The 
exponent v = \ belongs to the Ising universality class. 

Some physical insight into the meaning of the ground-state purity may be 
gained by noting that Eq. 15. 271 can be written in terms of the fluctuations of the 
total fermion operator N 

P,(|BCS)) = l-^((iV2)-(iV)^) , (5.34) 

where the | BCS)-property (cJ,Cfc/) = Sk^k'^l has been used. In general, the pu- 
rity relative to a given àlgebra can be written in terms of fluctuations of observ- 
ables IIBKUU3II . 
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Since fluctuations of observables are at the root of QPTs, it is not surprising 
that this quantity succeeds at identifying the critical point. Interestingly, by re- 
calling that Pso(2Ar)(|BCS)) = 1, the u(A^)-purity can also be formally expressed 
as 

P„(^)(|BCS)) = 1- (5.35) 

where the sum only extends to the non-number-conserving so(2A^)-generators be- 
longing to the set r specified in Eq. l5.26l Thus, the purity is entirely contributed by 
expectations of operators connecting different u(A^)-irreps, the net effect of corre- 
lating representations with a different particle number resuking in the fluctuation 
of a single operator, given by iV = Xlfc ^l^^- ^^S- 15 91 I show the probability 
Í2(n) of having n fermions in a chain of = 400 sites for 7 = 1. For g > 1/2 
the fluctuations remain almost constant, and so does the purity. 



1 

Q.(n) 
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Figure 5.9: Distribution of the fermion number in the |BCS) state for a chain of 
= 400 sites and anisotropy 7 = 1. 

Again, the isotropic case (7 = 0) is particular in the sense that P(, = 1 
(or = 0, see Fig. 15.71) . without identifying the corresponding metal-insulator 
QPT. The reason is that in this limit, the Hamiltonian of Eq. ISTSl contains only 
fermionic operators that preserve the number of partides (i.e., H E u{N)), and 
the ground state of the system is always a GCS of the u(A^) àlgebra. Therefore, 
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in order to obtain information about this QPT, one should look at algebras other 
than u(A^), relative to which the ground state is generalized entangled. A more 
detailed analysis can be found in Ref. IISOB04I . 

5.2 General Mean-Field Hamiltonians 

A Lie-algebraic analysis of many-body problems leads to a powerful tool for find- 
ing the spectrum and eigenstates of the Hamiltonian that describes the interactions 
in the system: 

M 

H = Y^ KjÒj, Kj G M, (5.36) 

i=i 

where the Hermitian operators Oj are linearly independent Schmidt-orthogonal 
elements of a Lie àlgebra [) = {Ói, ■ ■ ■ 0^} (Sec. 13.2.2b . Each operator Òj corre- 
sponds toadxd matrix acting on the Hilbert space H (of dimension d) associated 
with the physical system. In general, d is infinite or scales exponentially with a 
certain variable N, like the volume of the system, etc. For example, for a A^-spin- 
1/2 system one obtains d = 2^ . However, if the system is bosonic or is composed 
of harmònic oscillators, then d ^ oo. Throughout this section, I will assume that 
N < poly[log(d)]. 

Without loss of generality, one can consider [) to be a semi-simple Lie àlgebra 
acting irreducibly on H: The Casimir elements in H behave as constants (symme- 
tries) and do not change the physics. If the dimension of f) satisfies M < poly(A^), 
then H is defined to be a general mean field. Hamiltonian and is denoted as Hmf- 
Since f) admits a CW decomposition (see Sec. 13.2.2b . one obtains 

r / 

Hmf = ^ Ikhk + ^ tjÈa^ + L*È_aj , (5.37) 
k=i j=i 

where 7^ G M and lj G C are known (or can easily be computed) coefficients, 
hk ^ (i-C-j the CSA of [)), E±aj are ladder operators, and (r, /) < poly(A^). 

The eigenstates \(pP) of f)/) (i.e., the weight states of the Standard representation 
given by H) and the corresponding eigenvalues , satisfying 

h,\r) = el\r), (5.38) 

are assumed to be known with high accuracy. In particular, every weight state can 
always be obtained by the (efficient) successive application of lowering operators 
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to a highest weight state | HW) G Tí: 

i 

\<jf) = N,l[Èl^^^\H\N), (5.39) 

i=i 

where A^p G M is a normalization factor (nj G Z). If (ei, • • • , Cr) is the highest 
weight vector (associated with |HW)), then 



hm 



i 



(5.40) 



Therefore, all the information about the weight states of f) can be obtained from 
|HW), its weights, and the root vectors aj. This applies to every irreducible rep- 
resentation of f). 

The diagonal form Hd G P)d of Eq. l5.37l is the Hermitian operator 

r 

HD = Y,^khk = UHU\ (5.41) 
fc=i 

where U = e*^ is a unitary operator (F = F^^). In principle, U and Hd can 
be obtained by a simple matrix diagonalization algorithm working in some low- 
dimensional faithful representation of i). However, this is not sufficient to show 
that the eigenvectors and eigenvalues obtained correspond to a group operation 
and an element in [) d, respectively. Nevertheless, an important result in Lie theory 
States that an operator such as U can always be chosen to be a group operator 
induced by f), that is 

M 

^ = EOÒ, (5.42) 

Because Hd G fiD, its eigenstates are the weight states 10^) of Eq.|5321 Then, 

i 

\r) = U^r) = N,U^ W I HW) (5.43) 

i=i 

are the eigenstates of H. In particular, F can be chosen such that |HW) is the 
ground state of Eq . l5.41f and 

\G) = U^\H\N) (5.44) 

'Equivalently, one can choose F such that the lowest weight state |LW) is a ground state of 
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is the ground state of Eq. 15.371 Therefore, \G) is generalized unentangled and 
GCS of t). Moreover, considering that t/ is a similarity transformation, the spec- 
trum of Eq. l5.37l is then given by 



Hm 



In, 



(5.45) 



where the weights can be obtained from the highest weight vector as in Eq. l5.40l 
In brief, Eq. 15.371 is diagonalized if the coefficients of Eq. 15.411 and Q of 
Eq. I5.42l are obtained. This can be done in any faithful representation of f). To 
show this, assume that Oj is any faithful matrix representation associated to the 
operator Oj. From the definition of the exponential mapping, one obtains 



é^He-'^ = H + t[F,H] - - [F, [F, H]] 



and, equivalently, 



H + z[F,H]-^[F, [F,H]] 



k=l 



k=l 



(5.46) 



(5.47) 



which is the matrix representation of Eq. 15.461 Since M < poly(A^), then Hmf 
can be exactly diagonalized^ using the classical algorithm given in the foUowing 
section. 



5.2.1 Diagonalization Procedure 

Any element H of a real semi- simple Lie àlgebra () can be diagonalized through 
a similarity group transformation induced by f). Based on Ref. IIWÍ193I . in this 
section I introduce a classical algorithm for the diagonalization when f) is also a 
compact Lie àlgebra. This algorithm is a generalization of the well known Jacobi's 
algorithm IIPTV92I used for general matrix diagonalization. 

Assume then that $ is a faithful matrix representation of [) = {Oi, ■ ■ ■ , Om}, 
satisfying 

Tr [Ó,Ór] = (5.48) 

with Oj = ^{Oj). Equation 15.481 is always satisfied in the case of compact Lie 
algebras if working with the adjoint representation. In general, an orthogonal basis 

^By exact diagonalization, I mean efficient diagonalization. 



5.2 General Mean-Field Hamiltonians 



123 



is given by the Hermitian matrices hk, [Èaj + É^^j], and i[Èaj — -E-qJ (up to 
normalization factors), which are representations of the operators in the CSA of () 
and in the set of ladder operators, respectively. Therefore, the trace inner product 
of Eq. l5.48l allows one to project any operator in i) onto the sets ()+, and 
That is, the coefficients 7a; and lj of Eq. 15.371 can be obtained by projecting the 
matrix Hmf, associated to Hmf, onto the representations of the corresponding 
operators. 

To find the group operation U of Eq. 15.411 one starts by searching classically 
the Índex t such that \Lt\ > Vj G [1 ■ ■ ■/]. Second, one needs to perform a 
particular group (unitary) operation Ut such that the transformed operator = 
UtH^Ul with H° = Hmf, is 

r l 

fc=i i=i 

The operator Ut can be obtained by noticing the existence in f) of the su*(2) Lie 
àlgebra 

su\2)l Sí = [E.(«?)T'^-- (5.50) 

where are the components of the associated root vector. The operators S** 
and Sy are obtained from the relations S^ = |[S'*. ± iSy]. Since /i^, Èa^, and 

E^oij satisfy the commutation relations defined in Sec. 13.2.21 the operators S'* , Sy, 
and Sl satisfy the su(2) commutation relations of Eqs. IT?l (i.e.. spin operators). 
Equation l5.37l can then be written as 

H = Ht + Ht, (5.51) 

where Ht = u^S^ + u^^Sy + z/^S* G su*(2), z/^ G M, and is an element of the 
orthogonal complement su*(2)-'- as defined by Eq. 15.481 

Therefore, Ut is defined as the operator in the group SU^{2) = e""**^^^ that 
diagonalizes Ht G su*(2). That is, 

Ut = exp [tif^iSi + f^lSl + filSl)] , (5.52) 

where /i^ G M, and 

UtHtU} oc 5* G i)D. (5.53) 
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This is a diagonalization problem in su(2) and the defining coefficients /i^ of 
Eq. l5.52l can be obtained through the diagonalization of a 2 x 2 matrix^. 
Remarkably, Ut leaves invariant the decomposition 



(5.54) 



where í)^"^ and f)^"*" are the sets of operators in í)/? and (í)4. © f)_), respectively, 
which are orthogonal to su*(2). To show this, I first notice that 



í,-L 



W(2),fè- 



0, 



(5.55) 



or otherwise a group operation in SU^(2) = e""**-^-' could transform operators 
in í)^"^ into operators in 5U*(2). This is not allowed due to the orthogonalization 
property of Eq. I5.48Í For the same reason, 



su*(2),f) 



± 



(5.56) 



Equations l5.55l and l5.56l are sufficient to guarantee the invariance of the decompo- 
sition in Eq. l5.54l under the action of Uf. Then, such Ut anihilates the í-component 
of Hmf and transforms it as in Eq. 15.491 

The same su(2) -diagonalization procedure is applied to and so on, to ob- 
tain H"^, ■ ■ ■ , RP. In each step, a certain component t is eliminated as described 
above. Therefore, gets closer to a diagonal form and 



lim CÍ)D- 

p— >oo 



(5.57) 



To show this, I first write as 



(5.58) 



k=l 



where = Hmf and 7^ 
defined as 



Lj. The square distance of io is 



dem = E 



,P|2 



(5.59) 



^In fact, different sets of coefficients can be used to diagonalize this problem allowing one 
to choose the sign of the proportionahty coefficient 



5.2 General Mean-Field Hamiltonians 



125 



which is equivalent to count for the off-diagonal elements as in the Jacobi's diago- 
nalization algorithm for symmetric matrices IIPTV92II . In particular, dc{H^) = 
when G í)z). If í denotes the index of the biggest 1^^^"^! in the p — 1 iteration, 
then 



Y,\ir'?<i\ir'?. 

Moreover, since Ut leaves invariant (Eq. l5.54l ). then 



(5.60) 



,P\2 



Therefore, 



j=i{^t) j=i 



p-l\2 



< 



l - 1 



(5.61) 



p-l|2. 



that is, 



l - 1 



p— 1^ 



Since (^) < 1, then 



lim dc{H'P) = 0. 



p— >oo 



; (5.62) 

(5.63) 
(5.64) 



The group operation U of Eq. l5.41l is then 



U= lim Tf Ui 



í(p), 



(5.65) 



p=i 



where t{p) is the index for the biggest at step p. Nevertheless, Eq. 15.631 as- 
sures a ràpid convergence and U can be well approximated by a finite product of 
operators (i.e., P < oo). If the desired accuracy is denoted by e, the number P of 
required iterations is defined by 



dcil·l"") < e, 



which is satisfied if (see Eq l5.63l) 



/ - 1 



dcilÍMp) < e- 



(5.66) 



(5.67) 
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Naturally, the larger the dimension M of f) is (i.e., the larger l is), the larger the 
number P of iterations required to obtain such an accuracy. 

To obtain P as a function of M, I start by rewriting Eq. 15 .671 as (M = 2/ + r) 

loge > Plog(l - l/l) + \og{dc{HMF)) ~ -P/l + logidciHMF))- (5.68) 
Since Hmf ^ the following approximation can be performed: 

dc[HMF) M(fc (M>1), (5.69) 

where rf^ G M referres to some characteristic reference distance. Therefore, 
Eq. l5.67l reads 



P > l 



log M - log ^ 



(5.70) 



that is, the desired accuracy e in the diagonalization is guaranteed if P satisfies Eq. 
15.701 In particular, if M 1, the integer P is bounded above by a polynomial of 
second order in M (i.e., M > log M for M > 1). 

Equation 15 .701 is necessary but not sufficient to assure the efficiency of this 
diagonalization method. It remains to be shown then that can be obtained ef- 
ficiently. For this purpose, I consider a simple classical algorithm based mainly on 
Standard matrix multiplication. The idea is to work in a certain ç x g dimensional 
and faithful representation of [), such as the adjoint representation (AppendixO. 
In brief, the algorithm is based on two main steps: The search for the biggest l/fl 
and the diagonalization in su*(2). These steps are repeated P times to diagonalize 
H. If (where p E [1 ■ ■ ■ P]) is the matrix representation of H^, the biggest 
|if I is found by projecting the corresponding matrix onto the matrices È^^ and 
È-ap with j G [1 ■ ■ ■ /] (i.e., the representations of the ladder operators). This can 
be done with a conventional computer in 0{q^ï) computational operations (i.e., 
addition and multiplication of complex numbers). Once ú^ has been found, the 
corresponding su*(2)-rotation has to be performed. This operation is also repre- 
sented by a g x ç matrix and the representation H^'^^ of H^^^ can be obtained 
with 0{q^) computational operations. 

In brief, the coefficients defining in Eq. I5.58l can be obtained with com- 
putational accuracy in 0{Plq^) computational operations. From Eq. 15.701 and if 
M < poly(A^) and q satisfies q < M, then can be efficiently obtained"^ in, at 
most, poly(A^) operations. 

"^This evaluation is done at the computational accuracy given by the CC. Such an accuracy 
decreases polynomially in the number of computational operations, so the method is efficient. 
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5.2.2 Example: Fermionic Systems 

Consider a feraiionic lattice system with Hamiltonian given by 



H 



N 



(5.71) 



where is the size of the lattice, Xjy = {Xj'j)*, and cj (cj) the fermionic creation 
(annihilation) operators for site j. Such a Hamiltonian is known to be exactly solv- 
able and can be diagonalized through a Bogolubov transformation IIBR86II . Nev- 
ertheless, it constitutes a simple example to apply the methods described above. 

For this purpose, I start by noticing that if G u(A^), where u(A^) = {cjCjr] C 
[1 ■ ■ ■ A^]} is the Lie àlgebra of dimension M = N"^ introduced in Sec. 13.3.31 
Therefore, it can be efficiently diagonalized by working in a low-dimensional, 
faithful, representation of u(A^) such as the one given by 



'-ni 



(5.72) 



where Tjji are the N x N matrices with +1 in the jth row and j'th column, and 
zeros otherwise. In this representation, H has associated the foUowing matrix: 



A 



11 



A 



12 



A 



IN 



H 



(5.73) 



, Atvi AAr2 • • • Aatat 

which can be easily diagonalized in a conventional computer, obtaining 





fe, 


• 


■■ \ 


Hd = UHU^ = 





£2 ■ 


■■ 




U 


■ 


■ ■ Bn / 



(5.74) 



where f7 is a iV x unitary matrix. In this representation, Tjj/ span the whole set 
of X matrices, and any possible Ü is directly associated with a group element 
induced by u(A^). For this reason, the diagonalization procedure described in 
the previous section to assure a group transformation needs not be performed. 
Therefore, Eq. l5.74l defines the diagonal form of H through 



N 



N 



Hd = ^ BkTkk ^ ^ ^kclck = Hd, 



(5.75) 



k=l 



k=l 
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where 5^, the eigenvalues of H, are now the single-fermion energies. 

In this case, the diagonalization method used tums out to be exactly equivalent 
to the one based on the Bogolubov transformation for fermionic quadratic Hamil- 
tonians. It is also directly applicable to diagonalize a bigger family defined by the 
femiionic operators of the àlgebra so (2 A^) (Sec. 15. 1.^ . Additionally, this method 
can also be used to diagonalize more complex problems in which a Bogolubov 
transformation is not straightforward to implement. 

5.3 Summary 

In this chapter, I have addressed two important problems in condensed matter the- 
ory through a quantum Information theory point of view: The characterization of 
quantum phase transitions and the efficient diagonalization of many-body Hamil- 
tonians. It has been shown that the concept of GE, which becomes very useful 
when there is a Lie algebraic structure associated with the problem, plays a sig- 
nificant role in these cases. 

First, since QE is related to the existence of quantum correlations and these 
dominate the different phases, I have shown that by choosing a relevant set of ob- 
servables, the relative purity contains Information about the critical exponents of 
the phase transitions in two models of interest. Second, as motivated by the results 
of Chap.m I have shown that whenever an interaction Hamiltonian is an element 
of a low-dimensional semi- simple Lie àlgebra, it can be diagonalized efficiently 
through algebraic methods. This constitutes a major result in condensed matter 
theory and, in principle, it could be extended for the case of infinite dimensional 
Lie algebras which are generated by a finite set of functions HOSDOSI . 



Chapter 6 
Conclusions 



When is a quantum computer useful?; Which problems can be solved more ef- 
ficiently with a quantum computer than with a conventional one? Although not 
known yet, finding the answer to these qüestions constitutes the main reason that 
makes the science of quantum information a prospering and rapid-growing field. 
In this thesis, I have addressed various subjects in quantum information theory, 
including the quantum simulations of physical systems, quantum entanglement, 
and quantum complexity, to prové some of the capabilities of quantum compu- 
tation. The main idea can be briefly stated as follows: To efficiently simulate a 
physical system with a quantum computer, the laws of quantum mechanics need to 
be exploited to our advantage. Quantum entanglement, a non-classical property, 
is at the core of many tasks in quantum information and, if entangled states are 
created in a quantum simulation, such phenomena cannot be easily reproduced 
with a conventional computer. Therefore, it is expected that a computer which 
imitates the system to be simulated, i.e., a quantum computer, will be the most 
efficient device for this purpose. In the foUowing, I present the principal results 
and conclusions from each chapter of this thesis. 

In Chap. 121 I addressed several broad issues associated with the simulation 
of physical phenomena by quantum networks. I first introduced the conventional 
model of quantum computation (Sec. 12.1.11) as the main model used to perform 
these simulations. I studied the implementation of deterministic quantum algo- 
rithms (Sec. 12.21) which allow one to obtain relevant physical properties, usually 
related with the evaluation of some correlation function of the system under study. 
In general, the physical system one is interested in simulating is expressed by 
some operator àlgebra that may differ from the operator àlgebra associated with 
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the quantum computer (i.e., Pauli operators). I pointed out that efficient mappings 
between these two sets of operators exist in many cases (Sec. 12.31) and are suffi- 
cient to establish the equivalence of the different physical models to a universal 
model of quantum computation such as the conventional model. In Sec. 12.31 I 
explained how these mappings can be used to perform quantum simulations of 
fermionic, anyonic, and bosonic systems, respectively, in a quantum computer 
made of qubits. 

I also explored various issues associated with efficient quantum simulations. A 
simulation is said to be efficient when the amount of resources required is bounded 
above by a polynomial in some variable A^, such as the size of the physical sys- 
tem to be simulated. The main tòpics I addressed were how to reduce the number 
of qubits and elementary quantum gates needed for the simulation and how to 
increase the amount of physical Information measurable. I showed that the eval- 
uation of some correlation functions in efficiently prepared initial quantum states, 
like fermionic product states (Sec. 12.3.11) . can be efficiently done on a quantum 
computer. In some cases, this presents an exponential speed-up with respect to 
the corresponding classical simulation (Sec. 12.51) . where no known efficient algo- 
rithms exist (i.e., it would require an exponentially large amount of resources). 
However, it remains to be shown wether other tasks related with physical simu- 
lations on quantum computers can be performed efficiently or not. For example, 
there is no known efficient quantum algorithm to obtain the ground state energy of 
the two-dimensional Hubbard Hamiltonian. This is due to an exponential decay 
of the overlap between the efficiently prepared initial state and the actual ground 
state of the model (Sec. 1211). 

As a proof of principies, in Sec. 12.61 1 presented the simulation of a quantum 
many-fermion system, the Fano-Anderson model, using a liquid-state NMR based 
quantum Information processor. Relevant correlation functions were obtained 
by executing the quantum algorithms described in Sec. 12.31 For this purpose, a 
pulse sequence consisting of rf pulses acting on an ensemble of trans-crotonic 
acid molècules was performed. Moreover, different approximation and refocus- 
ing schemes were used to optimize such pulse sequence and minimize the errors 
of the simulation. The results obtained were very accurate because the overall 
duration of the simulation was much smaller than the decoherence time of the 
system. This quantum simulation was performed efficiently, i.e., with polynomial 
complexity in the system size. 

Although the studies on efficiency were done by considering the conventional 
model of quantum computation, the results obtained are independent of the actual 
physical representation of the quantum computer. A generalization of these re- 
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sults can be obtained by means of the notion of generalized entanglement, which 
was presented in Chap. This generalization of entanglement goes beyond the 
Standard subsystem-based approach, and is a feature of quantum states relative to a 
preferred set of observables of the system under study: it is an observer-dependent 
concept. To understand the main differences with the usual notion, I described the 
main properties of quantum entanglement in Sec. 13.11 In Sec. l3.2.2l l tied together 
the theory of entanglement and the theory of coherent states whenever the pre- 
ferred set of observables is a semi-simple Lie àlgebra. Important results were 
obtained in this case, where generalized unentangled states were defined as the 
extremal states of the àlgebra or generalized coherent states. Furthermore, these 
states present least uncertainty and can be considered as the most classical states 
in some sense. Some examples were presented in Sec. l3.3[ 

The main conclusion of Chap. |3lis then that conventional entanglement is a 
special case of a much more general theory, and this should be deeply analyzed 
to take advantage of the quantum world in different quantum Information proto- 
cols like quantum teleportation, quantum computation, etc. For this purpose, in 
Chap. m I studied some of the capabilities of generalized entanglement. Since 
traditional entanglement is known to be a resource for several tasks in quantum 
Information (Sec. 14.11) . including quantum computation (Sec. 14.21) . one would ex- 
pect that a more general theory of entanglement would allow one to better un- 
derstand the reasons lying behind the power of quantum computers. Therefore, 
in Sec. 14.31 1 presented a wide class of problems that can be solved efficiently on 
both, a quantum computer and a conventional one. These problems are related 
with the evaluation of a particular type of correlation function. In these special 
cases, the corresponding quantum simulation involves only generalized unentan- 
gled states relative to a certain polynomially large (or polynomially bounded) set 
of observables (Sec. 14.41 ). I showed that if no generalized entangled states, with 
respect to these sets, is created at some step of the quantum simulation, this task 
can be efficiently reproduced with a conventional computer. 

This important result indicates that although entangled states (in the usual 
sense) could be involved in the quantum simulation, this is not a sufficient con- 
dition to State that a quantum computer is more powerful than a classical one 
for these simulations. Nevertheless, if generalized entangled states relative to 
all polynomially large sets of observables are created in the quantum simulation, 
such phenomena cannot be easily reproduced with a conventional computer, and 
no known efficient classical algorithms exist in this case. This represents, again, 
an exponential speed-up with respect to the classical simulation. 

In Chap. El I discussed some issues related with efficient initial state prepara- 
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tion when simulating a physical system with a quantum computer made of qubits. 
These results have been generalized in Sec. l4.4.Í1 from the point of view of Lie 
algebras. Again, if the initial state to be prepared is the generalized coherent state 
of a Lie àlgebra with polynomially large (or polynomially bounded) dimension, 
such a state can be prepared efficiently; that is, it can be prepared by applying a 
polynomially large number of elementary gates. The type of gates depend on the 
physical representation of the quantum computer but the existence of one-to-one 
mappings makes this result independent of such representation. 

In Chap.|51 1 addressed two important tòpics in condensed matter theory: The 
characterization of broken-symmetry quantum phase transitions and the exact di- 
agonalization of Hamiltonians, from the point of view of generalized entangle- 
ment. In Sec. 15.11 1 showed that the relative purity, which constitutes a measure 
of generalized entanglement in the Lie algebraic case, successfuUy distinguishes 
between the different phases present in the LMG model and in the anisotropic XY 
model in a transverse magnètic field. In these cases, the corresponding ground 
States can be exactly obtained, so choosing the preferred set of observables that 
contains the relevant correlations in the different phases becomes relatively easy. 
Nevertheless, applying these concepts to a more general case can be done, in prin- 
ciple, by foUowing the same strategy. However, determining in a systematic way 
the minimal set of observables whose relative purity is able to signal and charac- 
terize the quantum phase transition requires, in general, an elaborate analysis. 

Finally, in Sec. 15.21 1 introduced the general mean-field Hamiltonians as those 
operators that are elements of polynomially large (or polynomially bounded) di- 
mensional Lie algebras. This is a generalization of the known mean-field Hamil- 
tonians such as the one composed by quadratic fermionic or bosonic operators. I 
pointed out that the existence of low-dimensional faithful representations guaran- 
tees the existence of efficient classical algorithms for their diagonalization (Sec. 15. 2.11) . 
In particular, the Bogolubov transformation is a special case of this type (Sec. l5.2.2t . 

Much remains to be done to really understand the power of quantum comput- 
ers. As is the case for other investigators in the field, I believe that a complete 
understanding of quantum entanglement for pure and mixed states is the key that 
will unlock such power. I hope that this thesis has been an interesting approach 
to the subject. In the foUowing section, I list a set of problems that may deserve 
further investigation but are out of the scope of the present work. 
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6.1 Future Directions 

Most of the results about efficiency in quantum simulations of physical systems 
(Chap. El) were based on the implementation of a particular type of algorithms. 
Nowadays, adiabatic quantum computation HFGGOÜI has emerged as an important 
tòpic in which it is expected that certain quantum states, such as the ground state 
of the two-dimensional Hubbard model, could be efficiently prepared by slowly 
changing some Hamiltonian interactions. It is important then to investigate wether 
there is a connection between the quantum complexity associated to both types of 
algorithms. 

Regarding the theory of generalized entanglement, several issues related to 
mixed- State entanglement also need further investigation. First, the natural exten- 
sion makes any measure, such as the relative purity, very hard to compute and 
more efficient ways for its evaluation need to be studied. Second, it is important 
to extend the results about efficiency to the mixed-state case. For example, the 
simplest case would be to consider mixed states which are a finite convex combi- 
nation of generalized unentangled states relative to a polynomially large dimen- 
sional Lie àlgebra. Here, the computation of some correlation functions would 
still be tractable (efficient) on a conventional computer. 

The classical algorithms described in Chap. H] were used to show that certain 
tasks, like the evaluation of some correlation functions, can be done with polyno- 
mial complexity on a conventional computer. However, it may be possible to find 
even more efficient classical algorithms for this purpose, and it would be worth- 
while to compare their complexity with that of the corresponding quantum one. 
In fact, many quantum algorithms found in the literature, like Grover's algorithm, 
do not have an exponential speed-up with respect to their classical simulation. 

In Chap.|5l I analyzed some quantum phase transitions from the point of view 
of generalized entanglement. It yet remains to be understood how to choose the 
set of observables that captures the most relevant correlations, which distinguish 
different quantum phases, in a more general case. That is, by just performing a Lie 
algebraic analysis of the interaction in the Hamiltonian, can we always distinguish 
the preferred set of observables that characterizes the phase transition? 

Finding the solution to these and other problems, like characterizing a bigger 
set of exactly-solvable systems in a Lie algebraic framework, constitute, together 
with the results obtained throughout this thesis, an advance towards the unifica- 
tion of Physics and Information Processing Theory. At the end, Information is 
physical. 
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Appendix A 

Discrete Fourier Transforms 



In practice, to evaluate the discrete Fast Fourier Transform (DFFT) one uses dis- 
crete samples, therefore Eq. l2.33l must be modified accordingly. From Fig. l2.12l it 
is observed that instead of having 5-functions (Dirac's functions), there are finite 
peaks in some range of energies, close to the eigenvalues of the Hamiltonian. Ac- 
cordingly, one cannot determine the eigenvalues with the same accuracy as other 
numerical calculations. However, there are some methods that give the results 
more accurately than the DFFT. 

As a function of the frequency rji, the DFFT (S{rii)) is given by: 

AI 

S{vi) = AtY,S{t,)e''^''^, (A.l) 
i=i 

where tj = jAt are the different times at which the function S is sampled ('Eq. l2.32l) . 
rii = are the possible freqüències to evaluate the FFT of S{t) \ and M is the 
number of samples. 

Since one is interested in S{t) = J2n l7",P e"*^"* (Eq. 12.321) . then 

N-l 

S{r]i) = Aí 5^ |7„|2 ^ e*[^"-^"]*. , (A.2) 

n j=0 

and 

n 

'Only a discrete set of freqüències can be obtained from the evaluation of the DFT over a 
discrete sample. In this case, the Nyquist critical frequency is given by i^c = 
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If r]i is close to one of the eigenvalues A„, and these are sufficiently far apart 
to be well resolved, all terms in the sum of Eq. IA.31 other than n, can be ne- 
glected. Taking rji and rji+i 

\S{r]i)l\S{vi+i)\ »0,then 



r]i + ^oth close to A„ in such a way that 



sivi 



3Í(r);+i-A„)Aí 



(A.4) 



After simple algebraic manipulations (and approximating e^ini->'n)^i ^ 1 + i(^r]i — 
Xn)At and the same for the denominator in Eq. IA.4I ) the correction to the energy 
is 

Xn = Vi + AA„ (A.5) 

with 



27r 



MAt 



Re 



S{vi+i) 



S{r]i) - S{r]i+i) 



(A.6) 



Appendix B 



Discrete Fourier Transform and 
Propagation of Errors 



Theoretically, the function S{t) of Eq. I2.86l is a linear combination of two com- 
plex functions having different freqüències: S{t) = |7ipe~*'^^* + |72pe~*^2*, 
where Aj are the eigenvalues of the one-particle eigenstates, defined as |lPi), 
in the Fano-Anderson model with n = 1 site and the impurity (Sec. 12.71) . and 
\i = |(0|lPí)p, with 10) = 11112)- However, the liquid-state NMR setting used 
to experimentally measure S{t) adds a set of errors that cannot be controUed, and 
the function S{t) shown in Fig. l2.2ll is no longer a contribution of only two dif- 
ferent freqüències. 

As mentioned in Sec. 12.71 S{t) was obtained experimentally for a discrete set 
of vàlues tj = jAt, with j = [1, ■ ■ ■ , M = 128] and Aí = 0.1 s. Its DPT is 
given by Eq. lA.ll Since one is evaluating the spectrum of a physical (Hermitian) 
Hamiltonian, the imaginary part of S{'r]i) is zero^ In Fig. 12.221 I show S{r]i) 
obtained from the experimental points S{tj) of Fig. 12.211 Its error bars (i.e., the 
size of the line in the figure) were calculated by considering the experimental error 
bars of S{tj) in the foUowing way: First, I rewrite Eq. lA.ll as 



with Qij = M-'^[Re{S{tj))cos{riitj) - \m{S{tj)) sm{r]itj)] (real). Then, the ap- 
proximate Standard deviation ES'; of S{r]i) depends on the errors EQij of Qij as 



M 




(B.l) 



'Due to experimental errors, the imaginary part of S{rii) could be different from zero. How- 
ever, I only consider its real part because it contains all the desired information. 
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(considering a normal distribution ||Tay97||) 



M 



(B.2) 



On the other hand, EQij is calculated as [ |Tay97| | 



[Eg.. 



dQij 


2 

Er^ + 




dRe{S{t,)) 




d\m{S{t,)) 



(B.3) 



where Er and E| are the Standard deviations of the real and imaginary parts of 
S{tj) (see Fig. l2.21l) . respectively. For experimental reasons (Sec. l2.7.TTl these er- 
rors are almost constant, having Er ~ E| ~ Es independently of tj (see Fig. l2.21l) . 
where Es is taken as the largest Standard deviation. Combining Eqs. IB .21 and IB .31 
we obtain 



ES, 



M 



1/2 



M 



(B.4) 



In the experiment, M = 128 and E^ fa 0.04, obtaining ES*; ^ 0.0035, which 
determines the (constant) error bars (i.e., the size of the dots representing data 
points) shown in Fig. 12.221 

The Standard deviation E'qi in the frequency domain is due to the resolution 
of the sampling time At. This resolution is related to the error coming from the 
implementation of the z rotations in the refocusing procedure. A bound for this 
error is given by the resolution of the spectrum; that is. 



Em < 



2tt 
MAt 



0.5 . 



(B.5) 



Appendix C 

The Adjoint Representation 



The adjoint representation of an M-dimensional abstract Lie àlgebra í) = {Ói, ■ ■ ■ Òm} 
is the transformation that maps every operator in í) into a M x M dimensional 
complex matrix given by the structure factors of f) IIGeo99ll . If 

M 

= ^J2fíA^ (c.i) 

fc=i 

where [, ] is the antisymetric form (e.g., the commutator) and fjj, are the structure 
factors, the matrices Oj given by the elements 

\\OJ\\,'k = -^f^r J,f,k=[l---M] (C.2) 

define the adjoint representation of f). In other words 

M 

[o„o,'] =z^/J;.,o, (C.3) 

k=l 

where [À, B] = AB — BA is the usual commutator between matrices. Since the 
operators Oj are Hermitian, the factors f^-, are real and the adjoint representation 
is pure imaginary. 

The Killing form is the bilinear form : f) x í) ^ M given by the foUowing 
mapping: 

K{Òj,Òj,) =Jr[ÒjÒj,]. (C.4) 

This mapping defines a convenient inner product in [). In particular, if f) is compact 
one can always choose a linear transformation of the operators in í) such that 

K{Òj,Òj>) = 5jj.. (C.5) 
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In this case, the observables Òj are said to be Schmidt-orthogonal. Moreover, if 
the àlgebra I) is semi-simple and Eq. IC.5l is satisfied, the adjoint representation is 
a faithful representation of f); that is, every operator Oj has associated a unique, 
linearly independent, M x M matrix. 

A group operation induced by f) is a linear unitary transformation of the form 



U 



M 



where H = W = CjOj e \), and Q G 
defined by 



JH 



l + iH 



(C.6) 

The exponential mapping e*^ is 

{iHf 



+ 



2! 3! 

Then, its action over an element of the àlgebra is given by 

1 

2 L 



0] = e 



-iH, 



H, O, 



H,0, 



(C.7) 



(C.8) 



M 



which also belongs to the àlgebra, that is, O' = X] ^jj'Oj' G í)- The real 



coefficients Ujji define m M x M dimensional matrix u, whose properties are 
given by the nature of the adjoint representation of [). To see this, consider the 
foUowing decomposition: 

U = l[UA = l[e''^'' (C.9) 

where the infinitesimal group operation Ua can be approximated by 

Ua ^ l + iAH ; A ^ 0. (C.IO) 

Naturally, this infinitesimal transformation induces another infinitesimal transfor- 
mation over i) given by the matrix 



where /i is an M x M dimensional matrix with real coefficients fijj' . 
Keeping the first-order terms in Eq. IC.81 one obtains 



(C.ll) 



M 



Oj-iA. H,Oj 



+ ^·l]/^ii'^i'· 



(C.12) 



The matrix /i is then given by the adjoint representation H of the element H as 
defined by Eq. lC.2l 

H = ijj, (C.13) 
and the infinitesimal transformation is given by 

1 + A./i ^ e^^-^ . (C.14) 

Since U is obtained after the successive application of Ua, one obtains 




(C.15) 



where the M x M dimensional matrix e*^ is obtained by exponentiating the ad- 
joint representation of H. Equation lC.15l is usually referred to the adjoint action 
of the group. In many cases\ the matrix e'^ defines the adjoint representation of 
the group induced by (). 



'That is, only vàlid for simply connected Lie algebras. 
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Separability, Generalized 
Unentanglement, and Local Purities 

For a quantum system S whose pure states \ip) belong to a Hilbert space H of 
dimension dim(7Y) = d, the purity relative to the (real) Lie àlgebra of all traceless 
observables P) = su{d) spanned by an orthogonal, commonly normalized Hermi- 
tian basis {Ài, ■ ■ • , Am}, M = d'^ — 1, is, according to Eq. l3.32l given by: 

M 

n(IV^)) = K5^(A.)'- (D.l) 

i=l 

The normalization factor K depends on d and is determined so that the maximum 
purity value is 1. If Tr(AiAj/) = óui (as for the Standard spin-1 Gell-Mann matri- 
ces), then K = d/{d — 1), whereas in the case Tr(AiAj/) = dóu^ (as for ordinary 
spin-1/2 Pauli matrices), K = l/{d — 1). Recali that any quantum state {tp) G H 
can be obtained by applying a group operator U toa reference state | ref ) [a highest 
or lowest weight state of 3u{d)]; that is 

|V) = t/|ref), (D.2) 

with U = 6*^»'''^% and Q E M. Therefore, any quantum state \ip) is a GCS 
of su{d), thus generalized unentangled relative to the àlgebra of all observables: 
Pf,(|V')) = lforall|^). 

Let's now assume that S is composed of N distinguishable subsytems, cor- 
responding to a factorization H = <S)f=i'^j^ with dim{Hj) = dj, d = Ylj^j- 
Then the set of all local observables on S becomes f) = í)ioc = 0, su·' ((ij). An 
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orthonormal basis which is suitable for calculating the local purity Pf, may be ob- 
tainedby consideringa coUectionoforthonormal bases su-'((ij) = {A{, ■ ■ A-^^}, 
Lj = cP- — l, each acting on the jth subsystem; that is, 



N factors 



i^' = al ® l2 ^ . . . ® A, ® ■ ■ ■ ® ]1^ , (D.3) 

factor 

where F = 1/ y^j- Then for any pure state ^Tí one may write 

N 



j=i i=i 



(D.4) 



By letting f)j = span{A^} be the Lie àlgebra of traceless Hermitian operators 
acting on Hj alone, the above equation is also naturally rewritten as 

^ 1 

n(i^)) = K'$^^n,(i^)), K,-^- (D.s) 

j=l 3 3 

The [)j-purity P\^. may be simply related to the conventional subsystem purity. Let 
Pj = Trj/^j({|z/;) (í/'l}) be the reduced density operator describing the state of the 
jth subsystem. Because the latter can be represented as 

P3 = j + Y.{M)^. (0-6) 

^ i=l 

Eq. lD.4l can be equivalently expressed as 



N ^ M 

(D.7) 



that is, P(jj(|í/^)) = {djJrp^ — l)/{dj — 1). Clearly, the maximum value of either 
Eqs. ID.SI or iDJi will be attained if, and only if, each of the conventional purities 



Tr Pj = 1 ^ Pf,^ = 1 for all j, which allows one to determine the K'-normalization 
factor as 

11 1 
K' = ^ = ' T = ^ . (D.8) 
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Accordingly, 

Pi,,jm = max = 1 ^ 1^) = ® ■ ■ ■ ® 10^) , (D.9) 

and the equivalence with the Standard notions of separability and entanglement 
are recovered. Note that for the case of N qubits considered in Sec. 13.3.21 the 
above value simplifies to K' = 2/N which in turn gives the purity expression of 
Eq. l3.5ll once the Standard unnormalized Pauli matrices are used (Àl = aÍ/V2, 
thus removing the overall factor 2). 
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Appendix E 

Approximations of the exponential 
matrix 

To classically compute the correlation function 

(HW|e-*^W^|HW), (E.l) 

one works in a low dimensional representation of the àlgebra f), such as the adjoint 
representation. Here, |HW) is a highest weight state of f) in some representation 
associated with the Hilbert space 7Y, and H and G f) are operators that map 
States in ü. into states in the same space. As before, I consider such an àlgebra to 
be compact semi-simple. Then, if [;) = {Oi, ■ ■ ■ , Ojv/}, there is a Killing form (for 
the adjoint representation) such that 

ií(0„0,)=5,„ (E.2) 

where Ok is the matrix representation of the operator Ok- 

In order to approximate a matrix, it is necessary to define a suitable norm. 
Here, I use the second norm defined by 

|A|=max|Ax| (E.3) 

where A is a n x n matrix and x E C" is a unit vector. If A can be diagonalized, 
this norm is equivalent to the largest eigenvalue. Nevertheless, the results and 
proofs of this appendix apply to any definition of matrix norm. 

The first step is to obtain, with the best possible method, a good approximation 
to the exponential matrix e*^, where H is the representation of the operator H. In 



148 



Approximations of the exponential matrix 



Ref. IIML03I . the authors state that the best method to evaluate such an exponential 
is the so called scaling and squaring method, which uses Padé approximants. This 
method is the one used in software like MatLab, etc. (However, if one is interested 
in evaluating e**^ for different vàlues of í, the method of obtaining the eigenvalues 
and eigenvectors of the Hermitian matrix H might be more efficient.) 
A Padé approximation to {A E C^") is defined by 

~ Rp,{A) = [D„{A)]-'Np,{A) (E.4) 

with 

''^ ^ Uíp + qV-jKp-jV ^ ^ ^ 

and 

Interestingly, Padé approximants can be used if \A\ is not too large. In this 
case, choosing p = q gives the best approximation. To calculate the matrices 
Npg or Dpg takes the order of qn^ flop operations, defined as the computational 
operation a ^ ax + y. The idea is then to use the property 

Therefore, m must be chosen such that it is a power of two and for which e^^^ 
can be efficiently computed. Then, the final matrix is obtained by m matrix mul- 
tiplications. 

A common criteria for choosing m is given by 

\A\/m < 1. (E.8) 

This is a criteria that might be too restrictive, but I will use it here. Then, m = 2^*, 
where s will be given by Eq. IE.8I In this way, the matrix e"^/™ can be efficiently 
computed by using a Taylor expansion or a Padé approximant. 

In particular, if c^/^' is approximated by Rgg(A/2'^), then q must be chosen 
such that the approximation has a small error. In the foUowing, I present some 
proofs obtained in the Appendix of the mentioned paper and which I do not de- 
scribe in detall here. 

First, if A is a matrix with \A\ < 1, then 

|log(]l + A)|< (E.9) 
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Second, if \A\ < 1/2 andp > 0, then 

\Dp,iA)-'\ < iq + p)/p. (E.IO) 

Third,if |A| < 1/2, then 

Rp,iA) = e^+^, (E.ll) 

where 

\F\ < 8|AP+^+^ (E 12) 

{p + qy.{p + q + ly.' 

The nx n matrix F can be shown to commute with the matrix A. This is because 
F must be a function of A, since the Padé approximants are functions of A, too. 
Then, if p = q, \ A\/2' < 1/2,1 obtain 

[RUAmr =e^^^, (E.13) 

with 

W.smV' ^ , <AV'- \„ (E..4, 



\A\ - (2g)!(2g+ 1)! - V2y (2g)!(2g + l)! 

Naturally, a low ratio between norms would give a good approximation. In 
the following, I relate these results with the specific problem of evaluating the 
correlation functions of Eq. lE.ll 

As mentioned before, I am interested in the approximation of the matrix e*^ 
(i.e., A = íH), where H is some low-dimensional matrix representation of the 
operator 

M 

H = Y.^jÒ,- OeM, He\). (E.15) 

i=i 

Here, M denotes the dimension of the compact semi- simple Lie àlgebra í). If a 
Killing form exists for the representation, I obtain 



\A\ = \H\ < \ MH 



M 



\ ECI- (£-16) 



Then, if each < do, where is some bound for the coefficients that build the 
operator H, 1 obtain 

\A\ = \H\ < ^ < ^/Mdo. (E.17) 
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Assuming that |yl|/2* < 1/2, then 

yMdo2'<l/2, (E.18) 

yielding to 



s > log2(VMdo) + 1. (E.19) 

Equivalently, 



m = 2' > 2VMdo. (E.20) 

Equation lE.20l tells one that the exponential can be approximated by an amount of 
products of matrices Rqg{iH /m) ~ e*^/™ efficiently in M. 
The error of the approximation is given by 

a = \[Rg,{iH/2')fW[Rg,{-tH/2')f ~W\ = {e'^'ie^We^' - W)e-'''\. 

(E.21) 

Because of the properties of the norm (i.e., \AB\ < \A\\B\ and |y4 + i?| < \A\ + 
\B\), one obtains 

a < |E|el^l|íy|el^'l + e'^l |íy | |Et|el^'l, (E.22) 

and considering that 

M 

i=i 

with ç, < do, then 



\W\< 



M 



^ç|<v^rfo. (E.24) 
\ i=i 



Evenmore, because E can be diagonalized (it commutes with H) we obtain \E\ = 
\E^. (This property might be satisfied for every matrix). Then, 



a < \E\e^^^^\^do. (E.25) 

Equation lE.25l tells one that for small vàlues of the approximation can be 
performed with high accuracy. Then, if 

Éi, (E.26) 



\E\ 


\E\ 


Jaï 


~ W\ 



E.l Scaling of the Method 
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I obtain 

a < el\H\e'^''\"\^/Mdo < eiMdle^'''^'^" < eiMdge^^i*"^ < e, (E.27) 

where e is the maximum tolerable error in the approximation. This error deter- 
mines ei, which determines the integer q of Eq. lE.14l In fact, for a constant error 
e, the higher M is, the smaller ei. Then q increases and the approximation needs 
to be done by higher order Padé approximants. 

Since I am interested in the case when M < poly(A^), where N < log(c?) is 
an integer that depends on the dimension d of the Hilbert space H, it remains to 
be shown that for fixed e, the integer q scales at most polynomially with M or A^. 
First, Eq. lE.27l tells one that for fixed e, I can consider 

eiM = r,,do (E.28) 

where d„ > is the coefficient of proportionality. Then (see Eq. lE.14ll 

) (2g)!(2g+l)! " " l^" ^^"^^^ 

[The previous analysis did not consider roundoff errors and these might be taken 
into account if needed. Nevertheless, one can consider that every step was done at 
the accuracy given by the number of bits of our computer (Turing machine). Also, 
some bounds to the errors can be improved.] 



E.l Scaling of the Method 

In this section, I am interested in obtaining the number of operations required to 
obtain cr < e as a function of norms of operators, etc. First, Eq. IE.14I can be 
bounded as foUows: 



Since 



I obtain 



a <2\E\e^\^\\W\, (E.31) 



16\H\ 161^1 - 16\H\ - 
(T < j:^e^^\W\ < j^:;^e\W\ <e, (E.32) 



(2g)2'í ' ' - (2g) 
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where I have assumed that < ^ 1- Since 

2glog2(2g) > 2q, (E.33) 

I obtain 

2g > 4 + log2 e + \og^ \H\ + log^ \W\ + \og^{l/e) (E.34) 

for the desired accuracy. 

Because calculating Npg or Dpg takes qM^ flop operations in the adjoint rep- 
resentation, calculating the approximated exponential matrix takes (2g + m)M^ 
flops. Then, to obtain the approximated matrix e~^^We^^ one needs n ~ 0[{2q+ 
m + 2)M^] operations. That is, 

n ~ [(log2 1/e + log2 1^1 + log2 \W\) + \H\/2]M' (E.35) 

flops, where I have considered e ^ 1. 
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Efficient Classical Evaluation of 
High-Order Correlation Functions 



In this case, I am interested in the evaluation of correlation functions of the form 



{HW\W'P---W'^\H\N) 



(F.l) 



where \HW) is the highest weight state of a compact semi-simple Lie àlgebra 
i) = {Oi, ■ ■ ■ , Oai}, where each Hermitian operator Oj maps states of the Hilbert 
space H into states in the same space. 

In a Cartan-Weyl decomposition, each operator W'"" can be decomposed as 



(F.2) 



í;=l 



Again, the roots Uj are considered to be positive so that the state | HW) satisfies 



K,|HW) 
/ifclHW) 



OVj, 

efclHW). 



(F.3) 
(F.4) 



Assume that 5^^ = (si, ■ ■ ■ , Sq) is a vector whose q components can take an 
integer value in the set [1, ■ ■ ■ , Z] (where / is the number of positive roots). Then, I 
want to show that 



EfiE_a,^ ■ ■ ■ E_ 



|HW) 



V x'^''^ ■ È ■■■È + 



154 



Efficient Classical Evaluation of High-Order Correlation Functions 



E 



x^''^ ■ È 



■Jq — 1 



|HW), 



(F.5) 



where P is also a positive root. To show this I use an inductive method. For this 
purpose, some commutation relations are needed. These are 



E/3, E^as^ 



Ep, E^as-, 



SI -I 



k=l 



= Nf3,-as,Ep.as, i-í P ^ asi- 

Then, for q = 1, the desired result is satisfied: 

r 

è^Í_,,JHW) = J]aJ^efc|HW)if/? = a,„ 

k=l 

E^E_,,JHW) = Oif(3>a,„ 



(F.6) 
(F.7) 



(F.8) 

(F.9) 
(RIO) 



where Eq s . IF 3 1 and IF.4I have been used here. («j > aj if they differ in a positive 
root.) 

I now assume that Eq. lFSl is vàlid for some value of q. Then, I need to show 
that it remains vàlid for g + 1 . In the latter case, I need to obtain 



Í2 ^^^^■E. 



Jl>---J9 + 1 -«Jl 



Í1>--J<7 + 1=1 



E_ 



X'' 



7f 



ji,---j<ï=i 



|HW)(F11) 



where «j is also a positive root and = (z, si, ■ ■ ■ , s^). Again, one can show the 
validity of Eq. lFTTl using the commutation relations of Eqs. lF.7l and lF6[ 
Consider then that (3 = ai. Therefore, 



(F.12) 



k=l 



155 



where each operator hk behaves as the constant (e^ — X]j=i Q;^,) = when 
acting over the corre sponding weight states. The validi ty of Eq. lF.l ll is obtained 
by noticing that the operator E^a- acting on the state of Eq. lF31 increases the order 

g in 1. For this case, the coefficients x"^" can be obtained from the coefficients 
xj^'... in the foUowing way: 



X ' — Z^k=l ^k, 



xf'^ = x^'^, (F.13) 



while the other coefficients remain zero. 
Consider now that (3 > a^. Then, 



■■■i5-a.JHW), 
(F.14) 

where at = (3 — ai is a positive root. Because of Eq. lFSl the first term on the rhs 
is a linear combination of states with order g + 1 in -È^q, while the second term is 
of order q. Then, Eq. lF.l ll is satisfied for this case, too. The coefficients can now 
be obtained by the foUowing recursion: 



13. 

x'l 



r/3, 5^ 



X^ 



X[ 



X' 



Nf3,-a^x2:..j„ forn e [1 ■ ■ ■ q],ji ^ i, 



(F.15) 



while the other coefficients remain zero. 

Similar results are obtained for the case /3 < ai. Then, 



(F.16) 

where ak = P — ai is a positive root. Again, Eq. lF.SI tells one that the first term 
on the rhs is a linear combination of states of order ç + 1 in E^a while the second 
term is just a weight state of order g + 1 in E^a- Therefore, Eq. IF.5l is satisfied 
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and the coefficients can be obtained as 

xy = 



X' 



X'. 



X' 



X 



forn e [l,---,^], 



(F.17) 



13 -Oi, 



while the others remain zero. 



For a fixed value of q, the number / of coefficients Xj^^.j^ V/?, is given by 

- 1 , 



/ 



/ - 1 



(F.18) 



and considering that q E [1 ■ ■ ■ p — 1], where p is the order of the correlation in 
Eq. IFll the total number F of coefficients Xj^^... ^^ V/5, satisfies 



F <p- 



IP+i 



l - 1 



-F+i < poly(/). 



(F.19) 



These coefficients can be computed easily, in polynomial time with respect to /, 
by using Eqs. lF13l IF151 and IF 171 The calculation of every coefficient is needed 
for the foUowing results. 

Remarkably, Eq. lFSI vields to the same results for the action of the operators 
W'' as defined by Eq.lF2l 



W"i ■ ■ ■w'^\mj) 



E 



E 



y1 E 



s<j-i 



Sl,---Sg-l = l 



|HW). 



(F20) 



The idea is then to update the coefficients zl ... whenever one of the operators W" 
is multiplied by the left of Eq. lF20l until q = p (see Eq. IFll) . The result is then 



{mi\w'p ■ ■ ■w'^\mj) 



(F.21) 



For this purpose, one needs to present a recursive method to update these co- 
efficients. First, if the operator '~fl^^hk acts on the state W''^ ■ ■ ■ iy'^|HW), the 
coefficients can be easily updated as 



(¥22) 



where n E [1 ■ ■ ■ q]. For a fixed value of q, the number of computational operations 
to update all the coefficients z'^ due to the action of a single operator 'yl^^hk is 
given by 



F,^ = 2 + ^/''(g + 2) = 2 



n=l 



l'^+\q-l) /'?+^ - / 



(F.23) 



However, if the operator {X^^^yÈ^a acts on the state W"^ ■ ■ ■ W'^\H\N), the 



coefficients need to be updated as 

9+1 



(F.24) 



while the other coefficients need not be updated or remain zero. The number of 
operations to update all the coefficients in this case is 



1 + 



1"+^ - 1 
l-l ' 



(F.25) 



Finally, when the operator \j E^^ acts on the state W"^ ■ ■ ■ Í^'^|HW), the 
coefficients need to be updated, too. For this purpose, I first rewrite Eq. lF.20l as 



^tíi ■■■ <^mE—a , • • • -Z?_Q, „ ) -|- 

f ^ í—^ •-'ml l'-'m *m ' 

m=l -çf 



|HW), 



(F.26) 

where the vector ~§ m = {S^^, ■ ■ ■ Sl^) has m components which can take vàlues 
in the set [1 ■ ■ ■ /]. Therefore, 



É^W"^---W'^\H\N)^YL=iT.^ 

O n 



zl, ™^ ... . x"/:.. .""-E.^, ■■■E_ 



Then, the coefficients are updated as 



x9+l 



^9+1 



A 



9+1 



V^9 9 ™aj,5^„ 



\9+i y^9 



lí,7É r. 



(F27) 



(F28) 
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Given that the coefficients are known, the calculation of the coefficients 

in Eq. lF.28l take ^-^j-j^ + ''''^^7'" computational operations. Therefore, the number 
of operations to update all the coefficients in this case is given by 

In brief, to update the coefficients due to the action of the operator W"^~^^ over 
the State W"^ ■ ■ ■ W'^ \ H W) it takes 

pg = ^pQ + + ^9) (F.30) 

computational operations (assuming that the coefficients x"^^ ' " are known). There- 
fore, the total number of computational operations JF to obtain the coefficients 
4i for the State W'p ■ ■ ■ W'^\H\N) satisfies 

< pPP, (F.31) 



that is, it takes at most a polynomially large amount of operations, with respect 
to / and r, to classically compute the correlation of Eq. lRTI (Though, the same is 
not true with respect to p. In fact, the complexity of the method described here 
increases exponentially with p, the order of the correlation function.) 
If M = r = 21 (i.e., the dimension of [)) satisfies 

M < poly(A^), (F.32) 

and ~ \og{d), where d is the dimension of the associated Hilbert space, the 
method presented here allows one to compute Eq. lF.ll efficientlv on a conventional 
Computer. 



Appendix G 

Classical Limit in the LMG Model 



As I mentioned in Sec. 15.1.11 some critical properties of the LMG, such as the 
order parameter or the ground state energy per particle in the thermodynamic 
limit, may be obtained using a semi-classical approach. Here, I sketch a rough 
analysis of why such approximation is vàlid (for a more extensive analysis, see 
Ref. ||SOB04J). 

Defining the coUective operators 

N 
k=l 

where a, a' =t or | and the fermionic operators c|,^ (c^g.) have been defined in 
Sec. 12.3. li The coUective operators satisfy the u(2) commutation relations; that is 

[E{a,a-'), E(^a-",a"')\ = S^' a" E[a,a"') — Sa-a'" E[a" ,a') ■ (G.2) 

If the number of degenerate levels N is very large, it is useful to define the 
intensive coUective operators E(^^„i) = E(^^y)/N, with commutation relations 

E{a;a-'), -E'(cr",íT"') 

Therefore, the intensive coUective operators commute in the limit N ^ oo, they 
are effectively classical and can be simultaneously diagonalized. Similarly, the 
intensive angular momentum operators Jx/N = (-É'(-f j^) + _É'(| |))/2, Jy/N = 
(%,|) - %,|))/2z, and JJN = - ^Qa))/2 (with J„ defined in Eqs.l53l 

15.41 and l5.51) commute with each other in the thermodynamic limit, so they can be 
thought of as the angular momentum operators of a classical system. 



— : [6a'a"E(^a-,a"') 



5nnlllEl„ll 



(o-",cr') 



(G.3) 
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Since the intensive LMG Hamiltonian H/N, with H given in Eq. (15.81 1. can 
be written in terms of the intensive angular momentum operators, it can be re- 
garded as the Hamihonian describing a classical system. The ground state of the 
LMG model \g) is then an eigenstate of such intensive operators when ^ oo: 
{Ja/N) \g) = ja\g) , ja bcing the corresponding eigenvalue. In other words, when 
obtaining some expectation vàlues of intensive operators such as Ja/N or H/N 
the ground state l^f) can be pictured as a classical angular momentum with fixed 
coordinates in the three-dimensional space (see Fig. l5.ll) . 

This point of view makes it clear why such operators ought to be intensive. 
Otherwise, such a classical limit is not vàlid and terms of 0{1) (order 1) would 
be important for the calculations of the properties of the LMG model. Obviously, 
all these concepts can be extended to more complicated Hamiltonians such as 
the extended LMG model, or even Hamiltonians including interactions of higher 
orders as in HGilSlI . 
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